Cosmic Ray Backgrounds for 
Dark Matter Indirect Detection 



Balliol College 
University of Oxford 

A thesis submitted for the degree of 
Doctor of Philosophy 
Trinity 2010 



Philipp Mertsch 




Cosmic Ray Backgrounds for Dark Matter Indirect Detection. 

Philipp Mertsch, Balliol College. 
A thesis submitted for the degree of Doctor of Philosophy. Trinity 2010. 

Abstract 

The identification of the relic particles which presumably constitute cold dark matter is a key chal- 
lenge for astroparticle physics. Indirect methods for their detection using high energy astrophysical 
probes such as cosmic rays have been much discussed. In particular, recent 'excesses' in cosmic 
ray electron and positron fluxes, as well as in microwave sky maps, have been claimed to be due 
to the annihilation or decay of dark matter. In this thesis, we argue however that these signals are 
plagued by irreducible astrophysical backgrounds and show how plausible conventional physics can 
mimic the alleged dark matter signals. 

In chapter 1, we review evidence of, and possible particle candidates for, cold dark matter, as well 
as our current understanding of galactic cosmic rays and the state-of-the-art in indirect detection. 
All other chapters contain original work, mainly based on the author's journal publications [1, 
2, 3]. In particular, in chapter 2, we consider the possibility that the rise in the positron fraction 
observed by the PAMELA satellite is due to the production through (hadronic) cosmic ray spallation 
and subsequent acceleration of positrons, in the same sources as the primary cosmic rays. We 
present a new (unpublished) analytical estimate of the range of possible fluctuations in the high 
energy electron flux due to the discreteness of plausible cosmic ray sources such as supernova 
remnants. Fitting our result for the total electron-positron flux measured by the Fermi satellite 
allows us to fix the only free parameter of the model and make an independent prediction for the 
positron fraction. Our explanation relies on a large number of supernova remnants nearby which 
are accelerating hadronic cosmic rays. Turning the argument around, we find encouraging prospects 
for the observation of neutrinos from such sources in km^-scale detectors such as IceCube. 

Chapter 3 presents a test of this model by considering similar effects expected for nuclear 
secondary-to-primary ratios such as B/C. A rise predicted above O(IOO) GeV/n would be an unique 
confirmation of our explanation for a rising positron fraction and rule out the dark matter expla- 
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nation. 

In chapter 4, we review the assumptions made in the extraction of the 'WMAP haze' which has 
also been claimed to be due to electrons and positrons from dark matter annihilation in the Galactic 
centre region. We argue that the energy-dependence of their diffusion means that the extraction of 
the haze through fitting to templates of low frequency diffuse galactic radio emission is unreliable. 
The systematic effects introduced by this can, under specific circumstances, reproduce the residual, 
suggesting that the 'haze' may be just an artefact of the template subtraction. 

We present a summary and thoughts about further work in the epilogue. 
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1. Introduction 



1.1. Dark Matter 

One of the most astounding results of modern astrophysics and cosmology is the finding that about 
20 % of the energy density of the Universe and about 80 % of the matter in the Universe is in the 
form of non-baryonic dark matter (DM). It turns out that DM plays an important role in structures 
on a wide range of scales ranging from faint satellite galaxies to the largest known structures of the 
Universe. DM, however, also presents a veritable problem for particle physics. Since the discovery 
of the discrepancy between the amount of luminous matter and total matter, many different pieces 
of evidence have been gathered, the most important of which we review in Sec. 1.1.1. In Sec. 1.1.2 
we discuss some relic particle candidates for DM, emphasising their connection to new particle 
physics beyond the Standard Model. 

1.1.1. Evidence for dark matter 

The probably most commonly known piece of evidence for DM (evidence at least for matter with 
a much larger mass-to-light ratio M/L than usual) comes from galaxy rotation curves, that is 
the variation of rotational velocity of stars with distance r from the galactic centre. Most of the 
luminous mass of spiral galaxies is in fact contained in the central bulge, and assuming Newtonian 
gravity and dynamics, the rotational velocity of stars should increase up to a certain radius rium 
that contains most of the luminous matter, and decrease beyond that like Xj \/r. What Vera Rubin 
however discovered [4] in 1970 for the nearby spiral galaxy M31, Andromeda, was ?>, flattening in the 
rotation curve beyond the radius numi implying that the Galaxy contains mass with a distribution 
different from the distribution of luminous matter. Flat or even increasing rotation curves have 
been also found for many other spiral galaxies [5], including the Milky Way (for a recent analysis, 
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see [6]). 

Already almost 40 years earlier, the amount of luminous matter in the Coma cluster had been 
determined [7, 8] by Fritz Zwicky, applying the virial theorem to measurements of the velocity of 
individual galaxies. He discovered that the amount of matter to explain the large velocity dispersion 
of up to 1000 km s^^ would require M/L up to hundred times the solar value. 

Another possibility of inferring the mass of a distant system, e.g. a galaxy cluster, is by grav- 
itational lensing (for a review, see [9]). This makes use of the bending of light from a distant 
bright source, like a quasar, by the gravitational potential in between the source and the observer, 
as predicted by General Relativity. If the distant source, the observer and the mass in between 
are aligned, the picture of the source becomes strongly distorted to multiple images or arcs which 
allows the determination of the mass of the object in between. Such constellations that lead to 
this so-called strong lensing are however rare and the effect on the shape of individual sources is 
mostly too small to determine the gravitational potential. It is however possible to analyse the 
systematic alignment of a set of background galaxies around the interjacent mass. This technique, 
called weak gravitational lensing, is able to resolve the (DM dominated) mass distribution purely 
by statistical means. The strength of the lensing signal compared to the intrinsic, so-called "shape 
noise" depends on the nature of the interjacent mass: Lensing by galaxy clusters usually gives the 
strongest signals and can be used in conjunction with observations of baryonic matter to determine 
the mass-to- light ratio of these systems. Lensing by foreground galaxies has much smaller effects 
(due to the smaller mass of galaxies compared to clusters) and usually requires combining, called 
"stacking", the signals from many galaxy lenses. Finally, it is also possible to determine the align- 
ment of background galaxy shapes along the large scale structure of the Universe, called "cosmic 
shear" and resolve the three dimensional (again DM dominated) mass distribution using also red- 
shift information. For instance, this allows the determination of the matter power spectrum and 
to constrain important cosmological parameters. The general advantage of gravitational lensing is 
that the mass of an object can be estimated without any assumptions about its dynamical state 
or even without any detectable baryonic matter in the first place. In the cases where such infor- 
mation is available, the inferences on the mass, from weak or strong lensing and from kinematic 
measurements, usually agree very well. 

Of course, the above observational results do not require that the non-luminous form of missing 
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matter is necessarily non-baryonic. In fact, a certain fraction of the non-luminous matter may well 
be in the form of heavy, compact objects, like black holes, neutron stars, white dwarfs or brown 
dwarfs, collectively referred to as MACHOs (MAssive Compact Halo Objects), and diffuse, hot, 
interstellar and intergalactic gas. However, the total estimated mass (including DM) of astronomical 
objects (galaxies, cluster, large scale structure) is much higher than the total amount of baryonic 
matter which can be determined from cosmological measurements as follows. The abundance 
of light elements as predicted by big bang nucleosynthesis (BBN) [10, 11] is a function of the 
baryon abundance 0,b where fli denotes energy densities in terms of the critical density pc = 
3H'^ / (SirG) with the Hubble parameter H and Newton's gravitational constant G. Observations of 
the (primordial) Deuterium abundance constrain the baryon abundance to « 0.02, where h = 
if/(100kms^^ Mpc^^) ~ 0.7. Baryons therefore only account for about 20% of the matter in the 
Universe. Furthermore, the power spectrum of cosmic microwave background (CMB) anisotropies is 
also sensitive to the amount of baryonic and the total amount of matter through the acoustic peaks. 
More precisely, the recent release [12] of the WMAP 7-year data reports l^cdm/i^ = 0.1109 ± 0.0056 
and Qbh'^ = 0.02258j:[|:[][}[j^^ for F = /i lOOkms'^ Mpc^^ = 71.0 ± 2.5kms-i Mpc'^ 

According to the concordance model, the same quantum fluctuations that lead to the anisotropies 
in the CMB are also the seeds for the formation of structure on all scales in the Universe today 
which provides more evidence for DM. On the one hand, the amount of density perturbations at the 
time of recombination is rather small (the CMB is isotropic to one part in 10^). On the other hand, 
structure formation through gravitational collapse could not start from baryonic matter before the 
time of recombination because of electromagnetic interactions. It turns out that the time that 
baryonic matter alone therefore could have had to form structure is much too short. Non-baryonic 
DM, however, would have started collapsing before recombination and it turns out it can explain 
structure formation rather accurately as shown by A^-body simulations [13, 14, 15]. Furthermore, 
structure formation requires that DM is non-relativistic at and after the time of recombination, 
that is cold dark matter (CDM). 

Other ideas, like Modification of Newtonian Dynamics (MOND) [16] (for a review, see [17]) 
and particular theoretical representations, for example. Tensor Vector Scalar (TeVeS) theory [18] 
can reproduce some of the above mentioned results without the need for new, non-baryonic DM. 
For example, MOND reproduces galactic rotation curves better than DM, in particular needing 
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less free parameters. Furthermore, it gives an explanation of the "Tully-Fisher relation", that is 
the observed correlation between luminosity and rotation velocity, but fails on the scales of galaxy 
clusters [17]. 

In particular, the 'smoking gun' signature for DM comes from the so-called 'bullet cluster' [19]. 
Weak lensing and x-ray expose the distribution of dark matter and hot gas (which makes up the 
majority of baryonic matter), respectively, in two galaxy clusters that have collided ~ 100 Myr 
ago. Not only are the distributions clearly distinct, but they also prove that while the hot gas 
has interacted (electromagnetically) during the collision, the DM has just passed through without 
interacting. Such a behaviour is basically impossible to imitate by modifying gravity. 

Despite all the successes of the CDM paradigm in explaining the above mentioned astrophysical as 
well as cosmological data, there are a large number of questions/problems that CDM cannot answer 
(so far) : Why are the DM haloes predicted by (pure DM) A^-body simulations much "cuspier" than 
those actually observed? Will including baryonic matter in A^-body simulations flatten the profile 
in the inner kiloparsecs? Why do we observe less satellites (DM substructure) than predicted by A^- 
body simulations? Can galaxy dynamics, e.g. SN feedback or (local) reionisation [20] reconcile this 
discrepancy? Why are there less massive galaxies today than expected? Although the quantitative 
discrepancy between the predictions and observations in some of these cases might be rather large, 
the expectation is that the uncertainty due to astrophysics is also still quite large and more realistic 
modelling will resolve many, if not all contradictions. 

1.1.2. Particle candidates 

What all the different pieces of evidence presented in the last section have in common, is that 
they are based on the gravitational interaction of DM only. However, from a particle physicist's 
viewpoint it is desirable to describe DM by a fundamental particle, in a similar way that the 
standard model (SM) of particle physics describes the fundamental matter particles and force 
carriers. The properties that characterise a particle in a relativistic quantum field theory are its 
mass, spin, quantum numbers under (gauge) transformations and also its coupling to other (SM) 
particles. 

Even before considering any particular model, the astrophysical evidence of DM gives us some 
constraints on its particle physics nature. First of all, the DM particle needs to be stable on 
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cosmological timescales, that is its lifetime must be larger than the age of the Universe. Secondly, 
zero results for searches for exotic, heavy nuclei on Earth constrain DM to be at most weakly 
interacting since DM would form such states if it was interacting via strong or electromagnetic 
forces. Furthermore, DM must be cold to be able to explain structure formation. This rules out 
the only known neutral and purely weakly interacting particle, the neutrino, as at least one mass 
eigenstate is relativistic today. Hence, none of the SM particles can accommodate for CDM and this 
clearly hints at physics beyond the SM. Interestingly, many theories of new physics that address 
the short-comings of the SM predict new particles including some therefore well-motivated dark 
matter candidates. 

WIMPs 

The SM has been very successful in providing a fundamental theory of matter and its interactions 
that has been tested to a great accuracy in colliders and non-accelerator experiments. However, 
although the Higgs mechanism explains the breaking of electroweak symmetry, one of the short- 
comings considered more serious is the so-called hierarchy problem: How is the weak scale stabilised 
with respect to radiative corrections that would normally boost it to the Planck scale? One of the 
most elegant solutions to this problem is supersymmetry (SUSY) which entangles the usual Poincare 
algebra with a new set of generators transforming fermionic into bosonic degrees of freedom and 
vice versa. In its simplest implementation, this effectively amounts to roughly doubling the particle 
content of the SM by mirroring each SM particle in a supersymmetric partner (sparticle) with a spin 
different by 1/2. These sparticles cancel the quadratic divergencies arising in radiative corrections of 
the Higgs mass parameter and hence stabilise the weak scale. To keep the necessary cancellations 
natural, the new particles must have masses mcM close to the weak scale. Among these new 
particles there are a number of DM candidates. 

If such new particles at the weak scale were also weakly interacting, that is they are weakly 
interacting massive particles (WIMPs), then their relic density today as predicted from production 
by freeze-out, i.e. thermal decoupling [21, 22, 23, 24, 25, 26], is in excellent agreement with the 
constraints from astrophysics and cosmology (in particular ilcdm from the WMAP experiment, 
see above). More precisely, at early times, the DM particles are in thermal equilibrium. As the 
Universe cools, their density is being suppressed by the Boltzmann factor oc e~"^oM/T ^^^^ £qj, j^^^g 
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times and hence low temperatures T, the abundance would normally vanish. However, apart from 
cooling, the Universe also expands, such that the rate of production/annihilation becomes smaller 
than the Hubble rate and the WIMPs drop out of thermal equilibrium. Quantitatively, the WIMP 
density n is governed by the Boltzmann (continuity) equation (first given in [21]), 

^ = -3Hn - (cJannf) {n^ " ^gg) (1.1) 

where H is the Hubble rate, {aannv) the thermal average of the WIMP annihilation cross section 
and velocity and Ueq the WIMP equilibrium density. In general, this needs to be solved numerically 
but a simple analytical estimate gives 

S^WIMP ~ 7-;— (<7annf )~ • (1-2) 

with Xf = "Hi-wiMp/T/ ~ 20, Tf (Tq) the temperature at freeze-out (today) and Mp\ the Planck 
mass. The WIMP mass does not enter Owimp directly but in most theories it is the only mass scale 
that determines the annihilation cross section. In particular, for weak interactions and a mass at 
the weak scale, one finds an s-wave annihilation cross section, 

(TAV- / a ~ 3 X 10-26 cm^ (1.3) 

that reproduces the relic density Owimp = C(l)- This fact has been called the 'WIMP miracle'. 

One might wonder why such a new, weak-scale particle should not decay, seeing that all SM 
particles above a GeV are unstable. Usually, one postulates a discrete symmetry that forbids 
interactions leading to decay of these particles. In the case of the minimal supersymmetric standard 
model (MSSM) this symmetry is called i?-parity. Originally, it was introduced [27] to satisfy bounds 
on the proton lifetime. However, to protect the proton it is enough to forbid one of the interactions 
necessary for its decay which can be achieved by, e.g. baryon or lepton parity. Therefore, to make 
the lightest supersymmetric particle (LSP) stable, i?-parity has to be introduced by hand. It was 
also argued [28] that i?-parity violating interactions must be suppressed to prevent the washout of 
the cosmic baryon/lepton asymmetry but this argument does not hold when lepton-mass effects 
are included [29]. 
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Examples for WIMPs arise naturally in weak scale SUSY theories as motivated above. If SUSY 
breaking is gravitationally mediated, the lightest supersymmetric particle (LSP) is usually the 
neutralino [30, 31], a mixture of the supersymmetric partners of the hypercharge gauge boson, the 
neutral component of the W boson and the neutral higgs partners. 

Theories of universal extra dimension (UED) [32] , for example, do not try to address the hierarchy 
problem, but also lead to new particles at the weak scale and thus WIMP DM candidates [33, 34]. 
The general idea is that all SM particles propagate in a higher dimensional space. In the simplest 
version, a single additional dimension is compactified on an S1/Z2 orbifold of radius R ~ 10^^ m 
or smaller. In 4 dimensions this leads to an infinite spectrum for each SM particle, so-called 
Kaluza-Klein (KK) [35] particles, equally separated in mass by R~^. KK particles posses a discrete 
symmetry, ii'-symmetry, which makes the lightest KK particle, the LKP, stable. In most models, 
the LKP is the first KK state of the hypercharge gauge boson, B^. 

superWIMPs 

On the one hand, the big advantage of the WIMP "miracle" is that it predicts the right DM 
abundance for stable, weak particles. On the other hand, constraining ourselves to stable and 
weakly interacting particles turns out to be too rigid a presumption. In fact, a number of the new 
particles suggested by beyond standard model (BSM) theories turn out to be either unstable or 
interact much more weakly. 

However, if every WIMP produced by thermal decoupling would decay to a "superWIMP" , that 
is a particle which is also at the weak scale but is only "super weakly" interacting, the virtue of 
the right relic density can be saved [36, 37], 

"superWIMP — "WIMP • (,l-4j 

"T-WIMP 

Another production mechanism is reheating [38, 39, 40, 41]. It can be shown that today's density is 
proportional to the reheating temperature Tr and one finds Tr ~ 10^° GeV for a 100 GeV gravitino 
superWIMP. 

Weak scale gravitinos are in fact a typical example of superWIMPS [41, 36, 37, 42, 43, 44, 45], 
realised in SUSY theories with gravity mediated SUSY breaking. As there is no reason to believe 
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that gravitinos are systematically lighter or heavier than the other superpartners in these scenarios, 
the gravitino is the LSP in about half the parameter space. The role of the WIMP is then played by 
the next-to- lightest SUSY particle (NLSP), for example the stau, which can decay to the gravitino 
with lifetimes naturally of the order of hours to months. 

Another possible example are axinos, the superpartner of axions (see below), and both could be 
possibly contributing in a multi-component DM scenario [46]. 

Hidden dark matter 

The fact that DM is not interacting through electromagnetic or strong forces has led us to the 
conclusion that it could at most be weakly interacting. The alternative, i.e. that it is only grav- 
itationally interacting is usually disfavoured because of lack of predictivity, missing connections 
to new physics and the loss of an automatic prediction for the relic density. Generally speaking, 
constraints on hidden sectors can only be obtained from their gravitational interactions, e.g. from 
constraints on the expansion rate at BBN [47]. There are however counter-examples resolving some 
or all of these problems and in the following we mention one particular class of examples that 
provide a DM candidate with the correct relic density from thermal freeze-out, too. 

In fact, going back to Eq. 1.3, one realises that every extension of the SM that predicts similar 
ratios of jm} will reproduce the observed relic density, no matter whether <^ is a SM coupling. 
Let's consider for example a setup of gauge mediated SUSY breaking which one could imagine 
is also mediated to a hidden sector whose particles are not charged under any SM gauge group; 
string theories, for example, predict many such hidden sectors. The masses of the hidden sector 
superpartners, m/i, are however set in the same, generation-independent way as the SM superpartner 
slepton masses, m, by the factor F/Mm from the SUSY breaking sector, 

"^~lfci£ and similarly rnr^-r^,^. (1-5) 

with Qh the hidden sector gauge coupling. Therefore, the ratio rrih/ gj^ that determines the relic 
density from freeze-out, see Eqs. 1.2 and 1.3, is universal and we expect that the hidden sector 
contains (some) DM candidates with just the right relic abundance. This relation is called the 
"WIMPless miracle" [48]. 
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Asymmetric dark matter 

It is questionable whether the freeze-out paradigm for DM production is necessarily the right one. 
Applied to baryon production, for example, it fails spectacularly - not only, because the predicted 
baryon abundance, ub/u^ ~ 10~^^, is about ten orders of magnitude smaller than the observed 
one, but also because it cannot explain the observed matter-antimatter asymmetry. Instead, the 
baryon abundance can be explained from an asymmetry of baryons over anti-baryons generated, 
for example, by non-perturbative sphaleron processes. If the annihilation cross section of baryons 
and antibaryons is large enough that all antibaryons will have annihilated away, we are left with 
the right density of baryons only. 

The fact that the abundance of DM and baryons is only different by a factor 5, a relation which 
in a freeze-out scenario of DM would be merely coincidental, may hint at a common physical 
origin [49, 50]. If one wants to explain the DM abundance today from an asymmetry generated 
in a similar way, DM must have a quantum number B\ similar to the baryon number B. At 
temperatures above the electroweak phase transition, T > T^,, electroweak anomalous processes are 
in thermal equilibrium and equilibrate the lepton, baryon and DM asymmetries. Depending on 
wether the DM mass m-oM is smaller or larger than T=k, the B' asymmetry gets frozen at Yab or 
Boltzmann suppressed as Ias'/^ab ~ e""^™/-^*, and the DM abundance, 0^/, can be related to 
the baryon abundance, Qb, 



with c an order one parameter and ms the nucleon mass. The observed DM abundance ^Ib' = 5 il^ 
can therefore be recovered by m-DM ~ 5 GeV for ttidm ^ [50] or ttt-dm ~ 2 TeV for moM ^ 
T* [49]. Technicolour, a strongly interacting theory suggested for electro- weak symmetry breaking, 
for example, rather naturally predicts iec/ini-baryons with TeV masses. 

Axions 

The axion was originally postulated to solve the so-called 'strong CP problem' [51, 52, 53]. The 
QCD Lagrangian contains a four-divergence proportional to the angle 9 that would lead to P and 




(1.6) 
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CP violation through non-perturbative effects. The non-observation of, e.g. a neutron dipole 
moment, however, constrains |^| to < 10^^*^. The most elegant way to guarantee a small value is 
the Peccei-Quinn mechanism that promotes ^ to a dynamical field with a classical potential that is 
minimised at ^ = 0. This is achieved by adding an additional, chiral symmetry f7(l)pQ, the Peccei- 
Quinn symmetry, which is spontaneously broken at the scale fa- The pseudo Nambu-Goldstone 
boson is the axion a; "pseudo" , because the global symmetry is not exact at the quantum level and 
therefore the axion is not massless but has a mass of order AqcD/Za, 

rUa^OSeY ^^'^f''^ . (1.7) 

Ja 

Also, this anomaly leads to a potential for the axion, fixing it to around (a) = —^/a/const, which 
cancels the above four-divergence, thus solving the strong CP-problem. Expanding the axion field 
around its minimum (a) one obtains the axion Lagrangian which is still to be complemented by 
the axion couplings to, e.g. photons and fermions, 

>Ca77 = -57-T-E • B and C^7r = igj—af-f^f , (1.8) 

where a is the fine-structure constant, E and B the (colour) electric and magnetic field, respectively, 
and 57, 9/ are model-dependent coefficients of order one. 

There is a variety of constraints on the axion mass. Collider searches for rare decays vr"*" — t- a{e~^e~)e 
rule out very short-lived (lifetime < 10^^^ s) axions with masses above 1 GeV. For long-lived ax- 
ions (lifetime > 10^^^ s) production {p + N^a + XoTe + N^a + X) as well as interaction 
{a + N ^ X) cross sections are constraint from beam dumps and rule out axions heavier than 50 
keV [54]. Astrophysical constraints [55] like bounds on the lifetime of red giants limit the cooling 
due to axions and give 200 keV > nia > 0.5 eV (> 10^^ eV even, if the coupling to electrons is 
large). The duration of the neutrino burst observed from SN 1987a finally also limits the fraction 
of axion cooling and excludes 2 eV > nia ^ 3 x 10^^ eV. Axions lighter than 10~^ eV are ruled out 
because they would over-close the Universe {Q,a oc /a oc m^^). The axion is thus extremely light 
and weakly interacting, making it a potential DM candidate [56, 57, 58]. 

Depending on whether the reheating temperature of inflation is smaller or larger than the tem- 
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perature Tpg of the spontaneous breaking of the PQ symmetry, the axion field is homogenised over 
enormous distances or carries strings and domain walls as topological defects. Axions therefore get 
produced by vacuum-realignment only [59] or also by string and domain wall decay, and the critical 
density today is 

with ai the initial misalignment angle. 

1.2. Cosmic Rays 

1.2.1. Galactic and extra-galactic cosmic rays 

The Earth's atmosphere is constantly bombarded by a flux of elementary particles: cosmic rays 
(CRs). The spectrum of these particles measured on Earth covers 12 orders of magnitude in 
energy ~ corresponding to 40 octaves in frequency! On the other hand, the steeply falling power 
law spectrum between E~^ '^ and implies quickly declining flux rates of 1 particle m^^ s^^ sr^^ 
above 100 GeV, 1 particle m~^ yr~^ sr~^ above 10^^ eV and 1 particle km~^ yr~^ sr~^ above 10^^ eV. 

Spectrum 

The cosmic rays with energies between ~ 1 GeV and (at least) ~ 3 x 10^^ eV are considered to be of 
galactic origin and are hence called Galactic Cosmic Rays (CGRs). They exhibit a rather featureless 
power law with spectral index a ~ —2.75. The low energy cut-off is due to solar modulation, that 
is cosmic rays of energy lower than a few hundred MeV lose all their energy by running up against 
the electric potential generated by the solar wind [60]. At ~ 3 x 10^^ eV, a softening to a ~ —3.1 
is observed, a feature called the "knee". At around 5 x 10^^ eV the spectrum further softens to 
a ~ —3.3 ("second knee"). As the Larmor radius at these energies starts exceeding the spatial 
dimensions of the Galaxy of kiloparsecs, cosmic rays at least beyond this energy must be of extra- 
galactic origin. At 3 x 10^^ eV, the spectrum hardens again to —2.7 ("ankle"), before it gets cut 
off at ~ 5 X 10^^ GeV [61], probably due to the so-called GZK cut-off, the suppression of the flux 
by photo-pion production on the CMB [62, 63]. The flux of cosmic rays measured on Earth, scaled 
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Figure 1.1.: Flux of galactic and extra-galactic cosmic rays, scaled by E'^''' (from [11] and with 
additional data [64]). The energy for the helium flux is in GeV/n. 

by E'^ '^ to amplify the features, is shown in Fig. 1.1. 

Below a GeV cosmic rays must be of local, i.e. solar system origin. In fact, the solar wind is a 
prime laboratory for testing models of cosmic ray acceleration (e.g. by shock waves) observed by 
interplanetary probes, like the Voyager spacecrafts. As these low-energy particles are however not 
important for dark matter indirect searches and cosmic ray backgrounds we do not consider solar 
cosmic rays any further. 

Composition 

Most chemical elements observed in OCRs occur in ratios that are consistent with their relative 
abundance in the solar system, hinting at a global validity of these abundances elsewhere in the 
Galaxy, potentially even the local Universe. Some elements however, in particular those which are 
rare in the solar system, are overabundant in cosmic rays. The interpretation is that although the 
source abundances are similar to solar system values, these so-called 'secondaries' get produced 
from the more abundant 'primaries' by spallation on the interstellar medium. This explanation 
gives in turn estimates of the average matter or 'grammage' of a few tens gcm~^ that primaries 
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must traverse to produce the observed number of secondaries. The energy dependence of the ratios 
of secondaries-to-primaries is an important test of the propagation model (see Sec. 1.2.7). It is 
worth noting that the electronic component has a softer spectrum (see Sec. 1.3.2) and is smaller 
than the nuclear one by ~ 10~^ at 10 GeV. 

The change of composition at spectral features is used to infer the nature of the sources of ultra 
high energy cosmic rays (UHECRs). Besides the observation of an increasingly heavy composition 
around the knee (see above), particular interest has been generated by recent contradictory obser- 
vations by the HiRes [65] and Auger [66] experiments at the ankle. However, it seems likely that 
this discrepancy could be resolved once hadronic showers at the highest energies are understood 
better. 

Anisotropies 

Another important information is of course encoded in the arrival directions of cosmic rays. Up 
to ~ 10^^ eV these directions get randomised by scattering on magnetic turbulence in the galactic 
interstellar (but probably also in the intergalactic) medium. The Auger collaboration has, for ex- 
ample, published evidence of a correlation of their highest energy events with a catalogue of nearby 
AGNs [67, 68]. The correlation has however decreased ever since, as have claimed correlations of 
other observations and source candidates. The question of the origin of UHECRs therefore remains 
open. 

1.2.2. The standard model of galactic cosmic rays 

Since their discovery by Victor Hess in 1912, a standard picture for OCRs has emerged that can 
explain a large number of observations. Usually, one assumes that the problem of GCRs factorises 
in two parts that can be treated separately: acceleration in confined sources and transport through 
the interstellar medium. This assumption seems to be supported by the observed decrease of 
secondary-to-primary ratios (see Sec. 1.2.7). In fact, if cosmic rays were mostly accelerated in the 
interstellar medium by some stochastic process (first order Fermi acceleration by interstellar shock 
waves or second order Fermi acceleration by interstellar turbulence) then these ratios would be 
logarithmically rising [69, 70]. 

The acceleration of both the hadronic (proton and nuclei) and the electronic component of 
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primary GCRs is assumed to take place in the shocked shells of supernova remnants (SNRs) by 
first order Fermi acceleration. Cosmic rays get accelerated as long as they are confined to the shock 
region by up-stream turbulence and only some high-energy particles might escape if no Alfven 
waves of low enough wave number are present for their scattering [71, 72]. In any case, once the 
SNR has entered its radiative phase all high-energy particles are released and diffuse through the 
interstellar medium by scattering on magnetic turbulence in the form of Alfven waves. It is expected 
that convection by a possibly CR generated wind and reacceleration due to stochastic, second- type 
Fermi acceleration also plays a role. 

In the following sections, we describe the basic building blocks in more detail. We start with 
a brief review of our understanding of supernova remnants, in particular their dynamics and the 
properties of the so-called Sedov- Taylor phase in which the bulk of particle acceleration happens. 
We then derive the transport equation and apply it to diffusive shock acceleration, claimed to be 
responsible for the acceleration of cosmic ray electrons and nuclei up to TeV energies. We also 
specialise the transport equation to cosmic ray propagation in the Galaxy, outlining the different 
physical processes at work and briefly review two simple analytical solutions which we extend 
on in later chapters. We conclude this section by reviewing the main predictions for fluxes and 
secondary-to-primary ratios and comparing them to observations. 

1.2.3. Supernova remnants as sources 

The assumption that the population of GCRs is powered by SNRs is based on three indications. 
The first one is the presence of non-thermal populations of electrons as observed in radio and x- 
rays, for example, for SN 1006 [73]. This is usually explained as synchrotron radiation of relativistic 
electrons on the ambient magnetic fields, amplified by compression and possibly by the Bell-Lucek 
mechanism [74, 75]. Estimates of the highest energies reach up to hundreds of TeV (for a modelling 
of SN 1006, see [71]). 

Furthermore, there is a well-established theory that can explain the generation of such a non- 
thermal population of particles by first order Fermi acceleration in shocked shells, see Sec. 1.2.5. As 
SNRs exhibit strong shocks with very high Mach numbers, applying this theory to the parameters 
derived from, e.g. kinetic observations, it can be shown that efficient particle acceleration is possible. 

Finally, even as early as 1953, it was pointed out [76] that SNRs with their benchmark total bulk 
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kinetic energy of 10^^ erg could easily provide the right order of magnitude energy needed to power 
the galactic population of nuclear cosmic rays. More precisely, with a SN rate of about 0.03 yr~^, 
a volume of the extended cosmic ray halo of 7r(15kpc)2 x 3kpc = 5.7 x 10^^ cm and an average 
residence time of cosmic rays in the Galaxy (that is the time until escape form the galactic cosmic 
ray halo) of 20Myr, O(10)% of the kinetic energy must be transferred into GCRs to maintain the 
(local) energy density of 0.3 GeVcm~^ - a reasonable efficiency. 

The spectroscopic classification of supernovae (SNe) into type I (without H lines) and type II 
(with H lines) does not coincide with the nature of the progenitor system. While only type la 
SNe originate in the thermonuclear burning of carbon-oxygen in a white dwarf, all other types (lb, 
Ic, IIP, IIL and Iln) are core collapse (CC) supernovae. Interestingly, the kinetic energy of the 
ejecta of the subsequent supernova remnants (SNRs) are quite similar, typically of the order of 
]^q51 g],g_ QnQQ the mass of the interstellar medium swept up by the shock front exceeds the mass 
of the ejecta, that is from the Sedov- Taylor phase (see below) onwards, there should not be any 
phenomenological distinction between the remnants of thermonuclear and CC SNe. 

Conservation of mass, momentum and energy across a planar, adiabatic shock front leads to the 
Rankine-Hugoniot relations [77]. In particular, for the compression ratio r = n\/n2 of the shock 
front, one finds, 

ti2 1 7 - 1 2 1 , , 

— = - = \ fl 10) 

and for strong shocks with a ratio of the specific heats 7 = 5/3 and a Mach number 7W S> 1, 
r = 4. Furthermore, this analysis is based on the assumption that radiative processes and loss of 
non-thermal particles can be neglected. Also, back reaction of the energetically important fraction 
of accelerated cosmic rays will change the shock structure, modifying 7 towards the fully-relativistic 
value of 4/3. 

Typical explosion velocities y^2£'sN/Mcj of the ejecta in the first phase of the SNR are 10"^ kms^^ 
for type la and 5000 km s~^ for CC SNe. This is much higher than the estimated sound speed in 
the surrounding medium, and therefore a blast wave with M > 10'^ forms, that is a shock front 
followed by self-similar {v{r) oc r), quickly cooling ejecta. After only a few days, when the shock 
front has decelerated sufficiently, the SNR enters the so-called "ejecta-driven" phase and a reverse- 
shock forms, that reheats the ejecta of the SNR. A contact discontinuity (with constant pressure) 
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develops between the heated ejecta and the shock front, hke a piston pushing into the ISM. This 
discontinuity is expected to be unstable to Raleigh- Taylor instability and therefore produces strong 
turbulence. 

The reverse shock will reach the centre of the SNR and disappear after hundreds of years which 
marks the transition to the Sedov- Taylor phase [78]. The time tch at which this happens can be 
estimated by dimensional analysis, assuming that the only scales available to the problem are the 
energy E, the mass of the ejecta Mej and the ambient density po, 

t,^ = E-'/^Mfp,'/-\ (1.11) 



The transition usually occurs when the mass swept up by the shock front is a few times the mass 
of the ejecta and therefore the value of M cannot enter into the subsequent evolution. Only using 
E and po, the dynamics of the Sedov- Taylor phase can be again be determined dimensionally [79], 

u=-- = 1.2xl0^ U^) - kms-S (1.13) 



5 t ylO^^ergy Vl.4/ Vcm^/ , 

where /U is the mean mass per particle in units of the proton mass. Extensive analytical and nu- 
merical simulations of the ejecta-driven and Sedov- Taylor phases have been presented, for example, 
in [79]. 

The Sedov- Taylor phase eventually comes to an end when the shock front has slowed down 
enough, so that radiative processes can become dominant. If the ejecta are hot enough at this 
point, the piston can still be powered by the pressure although the expansion is not adiabatic 
any more. If the ejecta have cooled too much, the former shock front just continues outwards, 
conserving momentum. Hydrodynamical simulations [80] show that under realistic assumptions 
this happens at 

^tr =2.9 X lOM -5, yr, (1.14) 



10^^ erg y Vcm^ 



E \ ^/ / \ — 7/17 
^tr=9l(T7^^) PC, (1-15) 
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M...^10M^V^'"'(^)-*'"Me. (1.10) 



After a few times ttr the velocity is as low as 100 — 300 km~^ and the Mach number is of order 
3 — 6 such that particle acceleration comes to a halt. 

1.2.4. The transport equation 

In the following we sketch the derivation of the transport equation governing the dynamics of a 
test-particle under magnetohydrodynamical turbulence in a moving plasma. We closely follow the 
treatment given in [81] which we refer the reader to for a more detailed calculation. 

The fundamental relation is the Vlasov equation [82, 83] for the Lorentz invariant phase space 
density / = /(t,x, p) at time t, 

We incorporate the processes of momentum, pitch angle and spatial diffusion by converting the 
microscopic Lorentz force term, encoding the electromagnetic interactions of the test-particle with 
the plasma, to macroscopic effective collision operators, thereby transforming the Vlasov into a 
Boltzmann equation. 

Momentum difFusion 

Scattering on magnetohydrodynamical turbulence leads to diffusion, both in momentum p and in 
pitch angle 9, fx = cos 6 = cos (p • B/(pi3)) where B is the large scale magnetic field. We can 
treat the diffusion in momentum space in a Fokker-Planck framework (see, e.g. [82, 84]), that is 
considering the scattering of charged particles on magnetic inhomogeneities or Alfven waves as a 
Markov process. The Fokker-Planck equation. 



|,,..V)/.|.D„.| wi. D„.l(^), a.lB) 



describes diffusion in momentum space and leads to second order Fermi acceleration as originally 
envisaged by Fermi [85]. Anticipating isotropy of the phase space density with respect to pitch 
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angle (see below), the RHS of Eq. 1.18 simplifies to 

y^PP^ • (1-19) 



p2 Qp PP Qp 



The momentum diffusion coefficient, Dpp, is calculated by considering that the momentum gained 
on bouncing off scattering centres in the plasma moving with velocity V, Ap = — (p • V)/f , occurs 
every At = L/v where L is the collision mean free path, 

J, _i / (P:V)!W_ip^ .,20) 

Pitch angle scattering 

The next necessary ingredient is pitch angle diffusion. Particles of momentum p with Larmor radius 
rL = p/{ZeB) interact resonantly with Alfven waves of similar wave-length, tl ~ 1/A;, changing 
their pitch angle 0. Interacting with waves with random phases, the pitch angle performs a random 
walk and after ~ (B/dB)'^ interactions it has changed by vr and the particle can be considered 
as scattered, that is having lost all information about its initial direction. Here, 5B denotes the 
turbulent component of the magnetic field. The Fokker-Planck equation for pitch angle scattering 
reads, 

(f).4|((-^>|). 

with the diffusion coefficient in pitch angle (i.e. the scattering frequency). 



\ / k£, 



At I 4 V^VSvr 



(1.22) 



and where E}^ is the energy in the mode k while B"^ /St: is the energy of the ambient magnetic 
field and Vt = eB/{'ymp) the Larmor frequency with nip the proton mass. Having assumed spatial 
uniformity, df/dx = 0, pitch angle diffusion can be implemented into the Vlasov equation 1.17 by 
identifying the RHS of Eq. 1.21 as an effective collision operator. 
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Spatial difTusion 

If the mean free path length is short, spatial transport needs to be treated in the diffusion equation. 
To prepare for the discussion of acceleration across a shock front, we also want to consider the 
possible motion of the scattering centres. The idea is to transform the Vlasov equation into the 
plasma frame and determine the isotropic part of the phase space density, /. It turns out that after 
averaging over pitch angle, / satisfies the transport equation, 

|[ + (u . V) / - V • V/) = ^ (V • u) , (1.23) 

where the coefficient for diffusion parallel to the magnetic field is 

1-2 f k8k , , 

^11 ^ — « ^ , , (1.24 

and the coefficient for diffusion perpendicular to the magnetic field, 

O.^i^^(^), (1.25) 



\B^/87T 
has been ignored. 

The second and third term on the LHS of Eq. 1.23 are the usual convection and spatial diffusion 
term and the RHS describes the energy loss (gain) of particles in a diverging (converging) flow. 
Adding Eq. 1.19 to the RHS completes the transport equation. 

1.2.5. Diffusive shock acceleration 

Fermi's original idea [85] for the acceleration of cosmic rays was stochastic scattering on magnetic 
turbulence in the form of "magnetic clouds". Unfortunately, the acceleration rate is only second 
order in the speed u of the scattering centres (hence the name second order Fermi acceleration): 
Although the average energy gained in head-on collisions is oc u/c, the average energy lost in 
an overtaking collision is equal but of opposite sign. In fact, only an asymmetry between the 
probability for either type of interaction, also of order u/c, leads to an overall energy gain, but only 
of second order, i.e. {u/cf'. It turns out that with the low speed of magnetic interstellar turbulence 



25 



downstream upstream 



n 



/o(p) 





a; 



a; 



Figure 1.2.: Left: DSA setup in the rest frame of the shock front, ui {U2) and n\ (712) denote 
upstream (downstream) plasma velocity and density, respectively. Right: Profile of the phase 
space density at a fixed momentum p. Particles within a momentum dependent distance D/u oi 
the shock front participate in the acceleration process. 

this mechanism would take much too long. Even for the much faster turbulence in the shocked 
shells of SNRs, second order Fermi acceleration is still too inefficient. 

The standard theory for cosmic ray acceleration is therefore a first order Fermi process in the con- 
verging flow of upstream and downstream plasma across a shock front: diffusive shock acceleration 
(DSA). For a contemporary review, see [86]. 

Macroscopic picture 

We consider steady diffusive shock acceleration (DSA) at a one-dimensional, parallel shock (that 
is the ambient magnetic field is parallel to the shock normal) in the test-particle approximation, 
that is ignoring the back-reaction of a potentially energetically important population of hadronic 
cosmic rays on the shock front. This was worked out almost simultaneously by three different 
groups [87, 88, 89]. (For the microscopic treatment of DSA [90], also published at the same time, 
see 1.2.5.) Here, we follow the treatment of [91]. 

The setup is as follows: The shock front is in its rest-frame in the y-z-plane at x = and 
upstream plasma is fiowing in from x < and downstream plasma fiowing out to x > (see left 
panel of Fig. 1.2). The density (velocity) of the background plasma are ni and n2 {ui and U2) in 
the upstream and downstream region, respectively. 

We now consider the transport equation 1.23 in this one-dimensional setup, however, ignoring 
the momentum diffusion term. 




(1.26) 
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where we have set D = Dy. The general steady state solution for x 7^ is, 



f{x,p)=gi{p)exp(^J dx^ ^^J^ +g2{p) ■ 



(1.27) 



^ /i (p) , where /i (p) is the phase space 



Imposing the usual boundary conditions, i.e. f{x,p) — 
density far upstream, that is the injection spectrum, and \f{x,p)\ < 00 for x — t- 00, we find 



fix,p) 



fi{p) + 91 (p) exp {Jq dx'u/D) for x < , 
92ip) = f2{p) forx>0. 



(1.28) 



The exponential term reflects the fact that only upstream diffusion can counterbalance the convec- 
tion by the background plasma in the upstream half-plane. The spatial dependence is illustrated for 
a homogeneous diffusion coefficient in the right panel of Fig. 1.2. Furthermore, the scale height of 
the exponential, di{p) = Di{p)/ui, is the distance that particles can on average diffuse away from 
the shock before being lost from the acceleration process. As the diffusion coefficient is growing 
with momentum, lower energy particles are confined more closely to the shock than high energy 
particles. 

To relate /i {p) , gi (p) and /2 (p) , we need to find matching conditions at the shock. These can be 
found most easily by multiplying Eq. 1.26 by a test function and integrating over an infinitesimal 
interval around x = 0. For the test-function dx' / D{x'), we find. 



[fte = 0{e) 



(1.29) 



that is / is continuous across the shock. For the test-function 1 on the other hand, 



D 



df 
dx 



df 



+ 3 ^^J- ^-dp 



0(e), 



(1.30) 



which for e — )• is the matching condition for df/dx. 

Substituting Eq. 1.28 into Eqs. 1.29 and 1.30, and eliminating gi, we find 



df2 



(1.31) 



27 



which is integrated to give 

/2(p) = ap-- r dpV»- Vi(p') + bp-^ , (1.32) 
Jo 

where a = 3r/(r — 1) is 4 for a strong shock, r = 4 (see Sec. 1.2.3). Ignoring the second term, which 
can be interpreted as acceleration of particles from the thermal background, it turns out that, as 
long as fi{p) is softer than the resulting spectrum /2(p) is a power law with spectral index a. 
The coordinate space density of accelerated particles is then n{p)dp = 47rp^/2(p) oc E~'^ which is 
the well-known result of DSA. 

Microscopic picture 

Although the macroscopic picture nicely proves that DSA leads to a power law spectrum oc E~'^, 
there is no physical intuition involved that could explain how the power law form comes about. 
Therefore, we re-derive the result of the last section, following a more physical microscopic ap- 
proach [90] (see also [92, 93]). 

As mentioned above, in second order Fermi acceleration, the test-particles can gain or lose en- 
ergy in approaching or following scatterings, respectively. In DSA, a first order Fermi process, 
however, the test-particles only encounter approaching scatterings: Seen from either plasma frame 
(downstream or upstream), the other side is always approaching. After crossing the shock, the 
test-particle quickly isotropises in the new frame experiencing only approaching scatterings which 
lead to energy gains. The energy gain rate is therefore proportional to the energy gained in every 
scattering, u/c. 

To calculate the spectrum we need in fact both, the momentum gain Ap of a test-particle in a 
cycle across the shock-front, i.e. coming from the upstream plasma, crossing to the downstream 
side and crossing back to the upstream side; and the probability that a particle gets advected 
downstream so that it is lost from the acceleration process. 

Consider a particle with momentum p, velocity v and pitch angle fx in the upstream frame. In 
the shock frame its momentum is p(l + fiui/v). Crossing the shock front, its momentum does not 
get changed, but in the downstream frame it is measured as p{l + fi{ui — U2)/v). To calculate 
its average momentum change on crossing the shock front once, we average over the pitch angle. 
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including the weighting factor 2/i, 

{Ap)=pl dnn{ui-U2)/v2ij,= '^p— — (1.33) 
Jo 3 V 

Crossing back, ui and U2 get swapped but the integration also runs from to —1, so the results is 
the same. 

The probability that a particle gets advected downstream can be calculated from the flux of 
particles to downstream infinity, n2U2, divided by the flux of particles crossing the shock front 
from upstream to downstream, djj, /if 71,2/2 = n2f /4. Hence, the probability that a particle gets 
advected to downstream infinity is 4:U2/v. 

We now consider the cumulative spectrum Ni{x,p) = dp' Anp'"^ fi[x , p) and assume that at 
upstream infinity all particles have the same momentum pQ, that is A'^i(— oo,p) = Nq6{p — pq). 
Since all the particles injected will eventually end up in the downstream region and can only have 
gained energy, N2{(x>,p) = vNq for p < pQ. 

The momentum pn and velocity u„ of a particle that has performed n cycles and has therefore 
crossed the shock 2n times, are of course random variables but for large re, they are expected to 
peak sharply around a mean (deterministic) value. Then, 

~ TT ( 1 + Q (""1 ~ ^2) /wj ) po log ( — ) ~ - (ni - 'U2) V — . (1.34) 

The probability Pn, that the particle has not been lost from the acceleration after re cycles is, 
^"~n(l- — ) ^ log (P„)~ -4^25^-. (1.35) 



1=1 2=1 

Combining Eqs. 1.34 and 1.35 gives 

log (P„) ~ -3^^ log (^) ^ Pn= (^-\ "'"'^^"'""'^ (1.36) 
The cumulative spectrum of particles downstream is then just the total density of particles 
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A''2 = tNq times the probability P„, 

-3u2/{ui—U2) 



N2(pn) = PnN2 = r ( ^ ) No ioT Pn > Po (1-37) 

Po, 



and 



f2[P) = --. = ~A ~ (1-^^ 

47rp^ op 47r — U2 \po J 

Although this result only reproduces the result of the macroscopic approach, Eq. 1.32, it becomes 
clear that as in every Fermi process, the power law spectrum is the result of the form of the energy 
gained on every crossing which is proportional to the energy of the incoming particle, and the 
probability to be lost from the acceleration process. In the case of acceleration at a shock front, 
both these number are fixed by the kinematics of the shock, as expressed by the compression ratio 
(and the Mach number). 

1.2.6. Galactic propagation 

The transport equation can, of course, not only be applied to the acceleration of cosmic rays 
but also to cosmic ray transport in the Galaxy. Expressing the original Eq. 1.23 in differential 
particle density n(t,x,i?) instead of phase space density /(tx,j)), n{t,j(., E)dE = Airp^ f{t-K,p)dp, 
and adding the injection of cosmic rays by SNRs as well as their production and losses due to 
fragmentation and decay, we find. 



dui d ^ d f dp p 

— - V . {D.. . Vn. - un.) - -p D„--^n. - - ^-n. - - 



g + (c /3 rigas (Tj-,i + 7 T-^-) rij - [c (3 rigs^s CTi + 1 T- ^) rii , (1.39) 



i<j 

where (3 is the speed in units of the speed of light c and 7 = (1 + /3^)^^/^. In the following, we 
discuss the different terms and their physical meaning in some more detail. 

• diffusion: —V • D^x • Vnj 

The main transport mode of cosmic rays is resonant pitch angle scattering on magnetic 
irregularities. Diffusion is strongly anisotropic {6B <^ B) locally, but the particle density gets 
isotropised by the fluctuations of the magnetic field at larger scales of 0(100) pc. The parallel 
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diffusion coefficient D^x ~ {^Bj-cs/ B)~'^vrg/3, with the Larmor radius and 6B^cs the field 
strength at the resonant wave number k = l/vg, is inversely proportional to the energy 
density in the turbulent component (see also Eq. 1.24), w{k)dk ~ k^^^^dk, which can be 
of Kolmogorov {S = 1/3) or Kraichnan {6 = 1/2) type. Expressing the diffusion coefficient in 
energy E, one finds Dxx oc I3E^ . Furthermore, with (5i?rcs ~ -B — few /xG at the principal scale 
of ~ (100 pc)"^, Dxxo = Dxx{lGeV) ~ lO^^cm^s"^. The effective values for 6 and Dxxo are 
usually determined from local observations of secondary-to-primary ratios (see Sec. 1.2.7). 

• convection: V • u nj 

There is evidence of a (CR driven [94]) wind in other galaxies and one might wonder whether 
SNRs also power a similar wind in the Milky Way [95, 96, 97, 98, 99, 100]. Direct observational 
evidence, however, only comes form x-rays close to the Galactic centre. A common way to 
achieve symmetry with respect to the galactic plane, u{—z) = u(z), is to consider a velocity 
linearly increasing with distance from the plane, u oc {du/dz)z. As the convection rate 
is decreasing with energy, convection can only play a role at lower, 0(1) GeV, energies and 
convection flattens out the steepening of source spectra (due to energy-dependent diffusion), 
observable for example in secondary-to-primary ratios. These also allow to constrain the 
convection velocity to n = 0(10) kms^^ (for one-zone models [101, 102]). 

• adiabatic energy losses/gains: (V • u)ni 

Diverging flows lead to adiabatic energy losses. This effect is for example important in 
propagation models with non-uniform galactic winds. 

• reacceleration: -jz-^P^Dpp-^^rn 

Substituting for the spatial diffusion coefficient, Eq. 1.24, in Eq. 1.20 gives the relation 
DxxDpp = p^v^/d where the Alfven velocity va is somewhere around 30kms^^. (In fact, 
the energy dependence of the diffusion coefficients can be modified if one accounts for the 
energy lost from the ISM turbulence into CRs which leads to damping and consequently a 
steep rise in the spatial diffusion coefficient at lower energies [103].) As already mentioned, 
distributed acceleration in the ISM cannot be the main source of acceleration because of 
the observed energy dependence of secondary-to-primary ratios. However, as the diffusion- 
loss time, z^g^^/Dxx, is decreasing with energy and the time-scale of distributed acceleration. 
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p^/Dpp oc Dxxlv\^ is increasing, second order Fermi acceleration by diffusion in momentum 
space, so-called reacceleration, can play a role at lower energies. For example, reaccelera- 
tion can show up in secondary-to-primary ratios and can potentially explain a bump in the 
boron-to-carbon ratio around 1 GeV [104, 105]. 

• continuous energy losses: —-§^^^1 

All CRs lose energy due to ionisation and Coulomb interactions which are however only 
important at energies below a few GeV. In addition, electrons and positrons interact with the 
ISM emitting bremsstrahlung (again, only important at ~ GeV energies), but also synchrotron 
radiation on the galactic magnetic fields and inverse Compton scattering (ICS) on interstellar 
radiation fields (ISRFs). In the Thomson approximation, the energy loss rate is proportional 
to the energy squared [106], 

|._1^^,,,,,,,,(_|)\ 

where cjt is the Thomson cross section. The energy densities for a 3 ^G magnetic field 
is 0.22 eVcm^'^; for the CMB, IR and stellar radiation they are 0.26eVcm~^, 0.2eVcm~^ 
and 0.45eVcm~^, respectively [107]. (See [108] for a sophisticated 3D modelling of ISRFs.) 
At higher centre-of-mass energies, the Thomson approximation is not valid any more and at 
even higher energies the cross section is in the Klein-Nishina regime and the loss rate becomes 
suppressed [107, 109] (see also [110]). For ICS, the critical energy depends on the average 
energy of the radiation background considered and is about 1.1 x 10^ GeV for the CMB, 
7.6 X 10^ GeV for IR and 8.7 x 10^ GeV for starlight [111]. 

• primary injection: q 

As discussed in Sec. 1.2.2, acceleration of primary CRs is supposed to take place in SNRs and 
the spectrum is a power law in energy (or rigidity) with spectral index close to —2 (in the 
test-particle approximation) . 

• nuclear spallation: Y^-^j /3 Jigas cij^i + 7 ^jllj) - {c l3 Ugs^s cri + j t~^) rn 
Spallation or decay of primary (secondary) cosmic rays during their propagation lead to their 
depletion and to the injection of secondary (tertiary) CRs. Furthermore, catastrophic energy 
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Figure 1.3.: Propagation setup. The thin disk of sources and interstellar gas is contained within a 
cylindrical cosmic ray halo of half-height ^max and radius R. Within the halo cosmic rays diffuse in 
coordinate and momentum space, get convected, spallate and lose energy, depending on the details 
of the propagation model. 

losses can in principle be accounted for by these terms. Note that spallation is dominantly on 
interstellar hydrogen and helium, i.e. ngas c'j-^i = iT'H + '^He '^f^i and similar for Tigas 
Parametrisations for meson production are discussed in Refs. [112, 113] and an overview of 
nuclear spallation cross section models is given in Ref. [114] and references therein. 

Usually, GCR propagation is considered to be limited to a cylindrical volume of radius R and 
half-height Zmax that the galactic plane is contained in, see Fig. 1.3. The transport equation is thus 
required to satisfy the boundary conditions, 

n(r|| ,z,t) =0 for |r|| | = R and for z = z-^aax , (1-41) 

where r|| = e^x + e^y . Alternatively, solutions without boundary conditions, but with a quickly 
increasing diffusion coefficient at the halo boundaries have been considered in [115] (see also [116]). 

Solar modulation is usually modelled by a simple effective (electric) potential <I> ~ 0(100) MV 
[60] which is assumed to be the effect of dynamical effects in the solar wind which are still not 
fully understood. The effect on charged CRs is that the flux at the top of the atmosphere (TOA), 
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Jtoa{E), is suppressed with respect to the interstellar flux Jis{E), 

Jtoa{E) = 2 "^isl^ + , (1.42) 

{E + \Z\e9)'^ — 

where {Z e) is the charge and m the mass of the CR particle. 
Leaky box model 

One of the earlier attempts to solve the problem of cosmic ray transport is the so-called Leaky Box 
model [117]. In this approach, all spacial dependencies are ignored and the Galaxy is described 
as a system that particles can escape from with a certain - in its simplest formulation, energy- 
independent - probability. In the steady state, this results in an exponential distribution of the 
column depth or grammage that cosmic rays particles have traversed. The average value of the 
order of a few gcm~^ must be a function of energy to explain the falling secondary-to-primary 
spectra. The simple Boltzmann-type equation describing this model, 

dnj rii Ui s—^ ( _i \ , _is 

-T7 = + 2^ (c/3ngasCTj^i + 7Tj_^i ) '^i - (c/3ngasO-i -F7Ti jni, (1.43) 

fl'' ''"esc ''"cool • , • ^ ' 

can be derived from the full transport equation by integrating over the cosmic ray halo and the 
contributions to the escape time Tesc from diffusion and convection are ~ -D^zZ-^max ~ "u/zmax, 
respectively, and the cooling time TcooI ~ p/(dp/dt). 

Applying the simple leaky box model to a network of stable nuclei, one can compute the rel- 
ative abundances of different isotopes assuming ratios at source that are similar to solar system 
abundances. For stable nuclei, it can be shown that the leaky box model reproduces the results of 
a diffusion model in certain limiting cases [118]. This can be understood considering that nuclei 
have a rather long residence time in the halo and therefore basically achieve the spatial averaging 
thereby justifying the leaky box approach. A generalisation of the leaky box model is the so-called 
weighted-slab technique [119, 120] that allows for more flexible path-length distributions. 



34 



Green's function approach 

A more sophisticated, but still analytic approach is to determine the solution to the transport 
Eq. 1.39 for a 5-like injection 6{t)6{r)S{E — Eq), that is, to find the Green's function G{t,r,E). 
The differential spectral density for a continuous source distribution is then given by its convolution 
with the Green's function. A particularly nice and useful example has been worked out for the 
propagation of electrons and positrons^ [121] (see also [122, 123, 124]), which we will also employ 
in Chapter 2. 

As we are mainly interested in energies above ~ 10 GeV, we can ignore convection (u = 0) and 
reacceleration {Dpp = 0), i.e. the transport equation 1.39 for the differential spectral density n of 
electrons and positrons reads, 

_ _ V . (Z) • V) n - — {b{E)n) = q{r, t) , (1.44) 

where the source term q now includes primary and secondary sources and the diffusion coefficient 
D = Z)|| is assumed to be homogeneous in the following. If one could ignore the energy loss term 
b{E) the energy of individual electrons and positrons would not change during their propagation, 
and the energy would simply be a parameter of a diffusion problem of the heat equation type and 
the corresponding Green's function is 



e-^'/^^* . (1.45) 



(47r£'t)3/2 



However, electrons and positron do lose energy and as the energy losses are continuous, the energy 
of a particle is monotonously decreasing from the energy at injection, Eq = E(t = 0), obeying, 

^ = b{E) = -boE' ^ l.-^ = -bQt. (1.46) 

where we limit ourselves to the Thomson approximation. Given the time t since injection and 
the energy E at observation, the energy at injection Eq is therefore unambiguously defined. As 
the diffusion coefficient is energy-dependent, we need to average it over time or equivalently over 



In the following, "electrons" is meant to denote both, electrons and positrons. 
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intermediate energies, E' . The diffusion length squared £^ = ADt is a function of both, E and Eq, 



\E,Eo)= / dt'iDiE{t')) =4 dE^ 
Jo JEo 



^ ,D{E') 



b{E') 



(1.47) 



With Eq. 1.46 and D{E) = DqE^ we find. 



4-Do f^s-i _ ^s-i 







(1.48) 



or with Eq = E/{1 - boEt), 



i\E,t) 



4Do 

b{l-6) 



E' 



1 - boEt J 



s-i 



(1.49) 



The contribution from an injection of electrons of energy £'o at t = 0, r = to the electrons of 
energy E at t, r is therefore. 



g{t,r,E) 



-r'^/f'(E,t) 



[lTl^{E,t)] 



3/2 



5 \ E 



En 



1 + boEotJ ' 



and a burst-like injection of a source spectrum Q{Eq) hence contributes. 



(1.50) 



-r2/e'2{E,t) / F \ 



To get the density for a steady source, this Green's function is integrated over time. 



G"*'^y(r,^) = / dtGit,r,E) 



bnE^ 



di^O . 3/2 ^(^0) ■ 



[7ri^{E,Eo)Y 



(1.52) 



This result can also be obtained from the time-independent transport equation, i.e. Eq. 1.44 with 
dn/dt = 0, see [125]. 

As electrons and positrons of tens and hundreds of GeV will lose most of their energy before 
travelling more than a few kiloparsecs and reaching the radial boundary of the cosmic ray halo, we 
ignore the radial boundary condition n{t, r = R, z) = 0. The boundary condition at z = Zmax can 
be implemented by the method of "mirror charges" [125], 
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X 



Gdisk{t,r,E)= G{t,rn,E) where r„ = y . (1.53) 

n=-oo \ {-!)"■ Z + 2Zjnaxn ) 

It is useful to factorise the spatial dependence into a dependence on (x, y) and z, r = ry + e^z, 

oo ^ 

Gdisk(t,r,S)= 7;^^^"'"^''^(^o)(l-&o^t)-' (1.54) 

n=— oo ^ ' 

= i^^"-"/'^Q(^o)(l - hoEtr-^x (f , ^) (1.55) 

where the sum has been expressed in terms of the ehiptic theta function, i?3, 

x(i,i) = ^[^3(^,e-^"')-i?3(i + |,e-^"')] , (1.56) 

and Zcr = 42;inax/7r. As most primary and secondary sources as well as our position are basically in 
the thin galactic disk, z and x{z^£) %(£) = x(0,^). 

The Green's function with boundary condition for a steady source is 



O^^E) = ^ I dE„^.-V^ Q(£„)_, ( ^ I . (1.57) 



(Semi-)numerical codes 

Both these propagation models oversimplify the problem but their assumptions are justified in 
certain limits: the leaky box model, for example, is accurate when considering stable nuclei, and 
the Green's function approach can be used for the propagation of electrons and positrons. 

A more realistic approach must take into account the morphology of the galactic magnetic field 
(GMF) and the interstellar radiation fields (ISRFs) as well as the resulting spatial dependence of 
the energy loss rate and the production of secondary radiation like diffuse radio and gamma-ray. It 
should further use the information on the gas density in the ISM, which does not only enter through 
the energy loss rate but also in the production of gamma-rays through 7r° decay. The resulting 
transport equation cannot be solved analytically any more and therefore numerical solutions are 
the only resort. 
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There are a number of numerical codes for GCR propagation in the hterature, some of them 
are pubhcly available. GALPROP [126, 127] numerically integrates the transport equation 1.39 on a 
spatial lattice by a Crank-Nicholson scheme. It can be run in a time-dependent mode or iteratively 
until a steady state, drii/dt ~ is reached. GALPROP uses input from 21-cm and CO (tracer of H2) 
maps for the ISM gas density and a sophisticated modelling of the interstellar radiation fields. It 
therefore allows not only calculations of nuclear and electron-positron fluxes but also more realistic 
predictions for gamma-rays and synchrotron radiation. The broad range of the predictions possible 
allows for important cross-checks, however, at the moment the full code is still not fast enough to 
allow for automated scans of the full multi-dimensional parameter space of the diffusion model. 

A similar, though slightly more general code is DRAGON [116] which allows, for example, for 
position dependent, anisotropic diffusion and separate injection spectra for different species. As an 
example for a semi-analytical approach we mention USINE [128, 129, 130] which employs analytical 
solutions found for a simplified setup with, e.g, simplified gas maps. This allows for sophisticated 
studies of the parameter space using MCMC techniques [131]. 

1.2.7. Observational results on the nuclear component 

In the following, we summarise the observational results in nuclear OCRs which are, to a large 
extent, in agreement with the predictions from the standard model of GCRs, as described above. 
We postpone the discussion of the leptonic component as well as synchrotron and gamma-rays, all 
of which have recently been claimed to show anomalies possibly connected to DM, to Sec. 1.3 where 
DM indirect detection is discussed. 

Primary nuclei 

As mentioned before, it has been shown [118] that for stable primary and secondary nuclei the 
diffusive-convective transport equation 1.39 is in some limits equivalent to a leaky-box model. In 
particular, further neglecting secondary production, energy losses and reacceleration, the ambient 
spectrum Ij of a primary species i is [118], 

li oc QiTTT? , 1— with — = , (1.58 
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where the column depth or grammage Xesc,i ^ ^/D is a power law in energy, oc E~ , above a few 
GeV/nucleon and 1/Xcsc,i dominates over the loss term ai/m. With a source spectrum Q oc E~~^ , 
7 ~ 2, the theoretical prediction for absolute primary fluxes is therefore li oc E"' with a = —7 — 5. 
From the local observation of primary nuclei (protons, Carbon etc.), we have a ~ —2.75 which 
would imply 5 ~ 0.75. However, from secondary-to-primary ratios (see below), one can determine 
6 independent of the source spectrum and finds 6 « 0.6 for plain diffusion, 6 ~ 0.4 when including 
reacceleration (see Sec. 1.2.7). This implies a source spectrum of 7 ~ 2.15 and 7 ~ 2.35, for plain 
diffusion and diffusive reacceleration, respectively, which is only marginally in agreement with the 
expectation from diffusive shock acceleration, cf. Sec. 1.2.5. At energies below a few GeV, ionisation 
energy losses need to be taken into account. 

In Fig. 1.4 we show a number of observations of the total proton and carbon flux, together with 
GALPROP results [103], for illustration. 




Kinetic energy, GeV Kinetic energy, GeV/nucleon 



Figure 1.4.: Left: Absolute proton flux in a diffusive reacceleration (DR) model (from [103]). 
The upper line shows the local interstellar (LIS) flux, the lower line the solar modulated one with 
$ = 550 MV and the thin dotted line the LIS spectrum best fitted to the data above 20 GeV [132]. 
Data: AMS [133], BESS 98 [134], CAPRICE 94 [135], IMAX 92 [136], LEAP 87 [137], Sokol [138], 
JACEE [139]. Right: Absolute carbon flux in a diffusive reacceleration model (from [103]). Data 
from ACE [140, 141], HEAO-3 [142], for other references see [143] (symbols are changed). 
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Stable secondary nuclei 

Secondary-to-primary ratios are used as test of the propagation model or to constrain its param- 
eters, like the spectral index 6 and normalisation Dxxo of the diffusion coefficient as well as the 
height of the diffusion zone -Sniax* However, as Dxx and -s^ax only enter into the fluxes of stable 
nuclei as the ratio D^x/ Zmax, there is a degeneracy between these quantities. 

The steady-state ambient spectrum of species i of secondary cosmic rays can also be understood 
in the framework of the leaky box model, 

Ij oc — ^ — -, ^— . (1.59) 

The numerator is dominated by primary species 1^ and above a few GeV/nucleon, Ik{E) oc E"'^^ . 
The denominator behaves similarly as for primaries, X^scj oc . Therefore, the secondary flux 
is softer than the primary one by E~^ and the secondary-to-primary ratio is falling, I2/I1 oc E~^ . 




Figure 1.5.: Left: Absolute antiproton flux from PAMELA and some other contemporary mea- 
surements [144, 145, 146, 147, 148] (from [149]). The dotted and dashed lines show the upper and 
lower limits for a range of diffusion models, taking into account the uncertainties on diffusion model 
parameters and cross sections [129] and the solid line shows the prediction from a different, plain 
diffusion model [103]. Right: Antiproton-to-proton ratio as measured by PAMELA and some 
other contemporary experiments [144, 145, 146, 147, 150] (from [149]). The dashed lines show the 
upper and lower limits for a leaky box model [151], the dotted lines for a diffusive reacceleration 
convection model [152] and the solid line the prediction for a plain diffusion model [103]. 
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Measurements of the antiproton-to-proton ratio and the absolute antiproton flux are shown in 
Fig. 1.5 together with some typical theoretical predictions. Figure 1.6 shows the B/C ratio, with 
predictions from the same GALPRDP model as in Fig. 1.4 and for a diffusive reacceleration model [103]. 
The range of spectral indices allowed by the data ranges from 0.3 for diffusive reacceleration to 
~ 0.6 for plain diffusion models. 




10"^ 10"^ 10° 10^ 10^ 10^ 10"^ 10"^ 10° 10"' 10^ 10^ 

Kinetic energy, GeV/nucleon Kinetic energy, GeV/nucleon 



Figure 1.6.: Boron-to-carbon ratio as a function of kinetic energy per nucleon in a plain diffusion 
(PD) model (left panel) and a diffusive reacceleration (DR) model (right panel) (from [103]). The 
lower lines show the local interstellar (LIS) flux, the upper lines the solar modulated one with 
^ = 450 MV. The experimental data below 200 MeV/nucleon are from ACE [140], Ulysses [153] 
and Voyager [154] and some high energy data are from HEAO-3 [142]. For other references see [143]. 



Unstable secondary nuclei 

Stable nuclei propagate over kiloparsecs before escaping from the cosmic ray halo, and hence the 
leaky box approximation is a good approximation for, if not even equivalent to, a diffusion model. 
For unstable nuclei the situation is however different. If their lifetime is smaller than or of the order 
of the average residence time in the Galaxy, their decays must be taken into account and escape is 
now competing with nuclear decay. As the survival probability, however, only depends on the time 
spent since production and not on the distance (the amount of matter traversed), secondary-to- 
primary ratios, like ^"^Be/^Be, can in principle break the degeneracy between diffusion coefficient 
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Do and halo height Zmax- At the moment however, the quahty of the data, in particular the 
statistics, is not good enough to really allow determining both parameters independently. Future 
cosmic ray experiments will provide better data, that will hopefully allow for the determination of 
either parameters, as both are in fact important for other observables, for example diffuse radio or 
gamma-ray backgrounds. 



1.3. Dark Matter Indirect Detection 

If the DM particle is in fact a weakly interacting particle produced by thermal freeze-out, see 
Sec. 1.1.2, then DM particles and antiparticles exist in equal amounts. By virtue of the weak anni- 
hilation cross section, the rate of annihilation is appreciable. A simple estimate of the annihilation 
rate per volume of a solar system DM density gives, 

2 "^DM - ^ ^ ' U.3GeVcm-3j llOOGevJ Vs x 10-^6 cm^ s'^ ' 

This amounts to a total rate of 1.6 x lO^^s-^ in a 1 kpc sphere around the Sun or 4.1 x 10'^^ s^^ 
for an NFW DM density profile, p{r) = 2pQ{r/rQy^{l + (r/ro))-^ [155]. 

Another possibility is DM decay. Although many DM models invoke a discrete symmetry to 
make the DM particle stable what is needed for an effectively stable DM candidate is actually only 
that its lifetime is much longer than the age of the Universe. In order to produce observable fluxes, 
the lifetime should however not be too large. A similar estimate as above gives, 

moMTDM^ lo.3GeVcm-3j llOOGevJ llO^^sJ ' ^^'^^^ 

Lifetimes of 10^^ s are naturally expected from, e.g. dimension-6 operators in grand unifying theories 
(GUTs) which are suppressed by two powers of the GUT scale Mqut = 2 x 10-*^^ GeV, 

^^^^^ 1026s^ (1.62) 

for mDM = 1 TeV. 

The particles produced this way, i.e. cosmic ray protons, antiprotons, electrons, positrons. 
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gamma-rays and neutrino, could in principle be observed on or around the Earth. Furthermore, 
annihilation products can produce secondary radiation like radio/microwaves from synchrotron or 
gamma-rays from ICS in addition to those from astrophysical CRs. The idea of indirect DM de- 
tection is then to extract the (possible) signal of DM annihilation or decay from the astrophysical 
backgrounds and use their spectral or spatial information to constrain the particle physics model 
of DM. 

Of course, a crucial question and the main purpose of this work is the determination of the 
contribution from backgrounds, i.e. the fluxes from purely astrophysical sources, to the new signals 
from the annihilation or decay of DM. On first sight, the constraint on these fluxes seem to be 
rather robust: According to the standard paradigm (cf. Sec. 1.2.2), the primary sources of OCRs 
are supernova remnants (SNRs) which are expected to produce power law spectra (see Sec. 1.2.5) 
and even after propagation, these fluxes should have rather featureless spectra. DM annihilation 
or decay is however expected to have very different injection spectra: if DM annihilated in a 2 — )• 2 
process or decayed to two particles, 1—7-2, the energy of the annihilation/decay products would 
be Mdm or Mdm/2, respectively. Hadronisation, for example of qq pairs directly produced or from 
gauge bosons, W'^ — t- qq' , Z ^ qq etc., would considerably broaden the energy spectrum. In any 
case, the injection would be much more concentrated in energy and even after propagation, the 
spectra should be different from the generic power law type of astrophysical origin. 

In the following sections, we will briefly review the ideas, prospects and current status of DM 
indirect detection in different channels assuming the standard astrophysical backgrounds as defined 
by the standard picture of GCRs. We present the claims of DM signature made in the charged 
lepton channels as well as the microwave and gamma-ray sky maps. 

1.3.1. Antimatter 

Using antimatter produced by annihilation or decay of DM for indirect detection harnesses the 
standard predictions for antimatter from astrophysical sources: a generally much lower abundance 
than matter and a particular energy dependence of the ratio of antimatter to matter. Due to their 
secondary nature, the spectrum of antimatter cosmic rays, like antiprotons or positrons, is softer 
than (primary) matter spectra and the secondary-to-primary ratios are falling above a few GeV, 
DC E-\ see Sec. 1.2.7. 
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Figure 1.7.: The positron fraction as predicted in different cosmic ray propagation models for a 
ratlier soft (left panel) and harder (right panel) electron spectrum (from [157]). The yellow band 
is spanned by the fluxes from a range of sets of propagation model parameters. Three particu- 
lar models, MIN, MED and MAX (named according to the effect on secondary antiprotons) are 
shown by the red solid, long dashed and short dashed lines. The data are from CAPRICE [158], 
HEAT [159], AMS [160, 161], MASS [162] and PAMELA [163]. 

In the following, we will briefly review the prospects for detection in two possible secondary- 
to-primary ratios, the positron fraction and the antiproton-to-proton ratio and quickly mention 
antideuterons. 

Positron fraction 

The positron fraction is defined as the ratio of the flux of positrons to the total flux of electrons 
and positrons (Je+ -|- Je-)i 

positron fraction = — — . (1.63) 

As mentioned above, in the standard picture of GCRs, positrons are only produced as secondaries 
by spallation of primary protons and nuclei which leads to lower abundances and softer spectra 
for the positrons and hence a falling positron fraction above a few GeV [156]. The benchmark 
prediction in [126] was one of the first applications of the GALPRDP code. 

The various sources of uncertainties and their impact on both the absolute positron flux and 
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the positron fraction have been assessed in [157]. The authors calculate the possible range in the 
absolute positron flux and the positron fraction from a broad range of models which differ in the 
primary proton fluxes adopted, the cross sections used and the primary electron spectrum assumed. 
The range induced by the uncertainty in the propagation model covers a factor of six between the 
lowest and the highest possible flux at 1 GeV but decreases to a factor 2.9 at 100 GeV. In particular, 
two different spectral indices for the locally measured electron flux are used, a = —3.35 (hard) and 
a = —3.53 (soft), and the positron fraction is always found to be falling for the hard electron 
spectrum and only very slightly rising for the soft electron spectrum, cf. Fig. 1.7. In fact, the 
(still preliminary) absolute electron flux measured by PAMELA [164], can be fitted with a power 
law with spectral index a = —3.226 it 0.020 between 10 and 60 GeV and is even harder below. 
Therefore, the soft electron spectrum is ruled out and even within the rather large uncertainties, 
the positron fraction from secondary production is always falling. 

This is however in disagreement with the recent findings [163] from the PAMELA experiment 
which shows a rising positron fraction above ~ 5 GeV, see Fig. 1.7. This trend was already apparent 
in the combined (1994, 1998, 2000) HEAT data [165], although with much larger uncertainties. Even 
more than 20 years ago, the rise in the positron fraction was realised [166] and tried to explain with 
some exotic contribution, possibly from pulsars. 

Of course, in particular in connection with claims of excesses in the absolute electron plus positron 
flux by PPB-BETS [167], ATIC [168] and Fermi-LAT [169] (see also Sec. 1.3.2 below), this result 
has generated much attention, mainly prematurely interpreted as evidence for WIMP annihilation. 
In fact, within a year of the presentation and publication of the PAMELA data, a multitude of 
models was put forward (see, for example, [170] for a comprehensive, though non-exhaustive list of 
references) that could explain the upturn in terms of the onset of a new, harder primary positron 
component from DM annihilation or decay. 

The contribution from DM annihilation can be taken into account both in (semi) analytic and 
in fully numerical computations by an additional source term on the RHS of the transport equa- 
tion 1.39, 



Here, (o'ann'y) denotes the averaged annihilation cross section times velocity, p^{r,z) the DM halo 




(1.64) 
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profile and g{E^±) is the differential production spectrum. The latter encodes the particle physics 
(obviously dependent on the DM model considered). Adding this flux to the one from astrophysi- 
cally produced secondary positrons, predictions for the positron fraction can be made. 

However, with a thermal annihilation cross section, (cianni') ~ 3 x 10~^^cm^s~^, and standard 
astrophysical assumptions about the DM halo, e.g. NFW [155] or isothermal [171] profile, local DM 
density of ~ 0.3GeVcm~^, it turns out that the additional fluxes fall short by between a factor of 
10 to 1000, depending on the DM mass (see, e.g. [172]). It would be very difficult to detect such 
small deviations from the astrophysical backgrounds, in particular considering the uncertainties 
involved. 

In fact under certain assumptions one can expect an additional "boost factor" to appear in 
Eq. 1.64, either from astrophysics or particle physics. Astrophysical boost factors are induced by 
overdensities in the distribution of DM as for example predicted by A^-body simulations. This does 
not only lead to a global amplification compared to a smooth DM density, J dV > {J dV p)'^, 
but in particular to potentially "bright" nearby clumps of DM. However, the probability of such 
a DM overdensity close enough to the Earth to explain the excess in the positron fraction is 
very small [173]. Particle physics boost factors can come from a low- velocity enhancement of 
the annihilation cross section, for example by a resonance just below 2Mdm [174] or so-called 
Sommerfeld enhancement [175, 176, 177]. It is also possible to overcome the helicity suppression 
of the s-wave annihilation into two fermions by emitting a photon. Alternatives include giving up 
the WIMP paradigm by, e.g. assuming non-thermal DM production (see [178] for the example of 
Wino-like neutralino DM). 

In terms of annihilation channels, the PAMELA positron fraction alone allow annihilation to 
charged leptons for a wide range of DM masses, rriDM ^ 100 GeV [172]. For light DM, moM ~ 
100 GeV, annihilation into W bosons is also possible. 

Absolute antiproton flux and a ntiproton-to- proton ratio 

For a generic WIMP model, there is a priori no reason why DM should only annihilate or decay 
into electrons and positrons or other charged leptons. In fact, heavy quarks possibly produced by 
DM annihilation or decay will hadronise to all sorts of baryons, ultimately also yielding protons 
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and antiprotons. The antiproton-to-proton ratio, 



antiproton-to-proton ratio ratio 



Jp 
Jp 



(1.65) 



and the absolute antiproton flux can therefore be used as a cross check for DM explanations for 
the positron excess. 

First of all, one should note that predictions for the background of antiprotons from the spallation 
of (mostly) GCR protons on the interstellar H and He are able to reproduce the measurements by 
BESS [147], CAPRICE [145] and AMS-01 [148] on the absolute antiproton flux as well as the 
most recent PAMELA data (see Sec. 1.2.7 and Fig. 1.5). A detailed study of the background 
investigating the uncertainties from the cross sections used and the diffusion parameters adopted, 
has been presented in [129]. 

The potential contribution from the annihilation of DM into antiprotons can be calculated within 
the theoretical frameworks presented in 1.2.6 or using fully numerical codes like GALPROP, starting 
from a spatially varying injection term, 



Here, g{Ep) is the differential production spectrum, assembled from the branching ratios Bh into 
quarks or gluons h in different channels F, as well as the fragmentation and hadronisation functions 



where Ep is the antiproton's kinetic energy. 

Considering again the uncertainties introduced by the GCR propagation parameters but also by 
the DM halo model adopted, the semi-analytic analysis of [179] finds a large proportion of a scan 
over a particular MSSM parameter space to be consistent with the antiproton fluxes. We note that 
in this study no boost factors were introduced. 

However, antiproton measurements can give stringent constraints [172, 152] on a number of 
models invoking boost factors to explain the anomalies in the lepton channels, e.g. the positron 
excess. Lepton channels are still available at all masses, but the antiproton data basically exclude 




(1.66) 



diV^/dSj 




(1.67) 
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hadronic annihilation channels with m-DM ^ a few TeV. Above 10 TeV, both leptonic and hadronic 
channels (excluding, perhaps, direct annihilation into quarks) give both good fits to positron and 
antiproton data. For light and intermediate masses, however, annihilation into gauge and Higgs 
bosons must somehow be suppressed. For an example of a model building way out of this, see [180]. 

Anti deuterons 

Following the idea of DM indirect detection in rare nuclei it was suggested [181] to look for anti 
deuterons D from the annihilation of DM. Secondary antideuterons get produced by GCR p, He and 
p impinging on interstellar H and He, with pH and pile interactions dominating but contributions 
from p becoming comparable below 1 GeV/n. Uncertainties from the cross section and again from 
the cosmic ray model amount to an order of magnitude at ~ O.lGeV/n, decreasing to about a 
factor four at ~ lOOGeV/n [182], similar to the behaviour for the antiproton flux [129]. 

Although currently only upper limits on the D flux exist [183] , it turns out that with reasonable 
assumptions for the propagation model the ratio of DM signal to astrophysical secondary back- 
ground is usually larger than 0.5 below IGeV/n for DM masses up to hundreds of GeV. Anti 
deuterons are therefore one of the most promising detection channels for light and intermediate 
WIMP masses, and it has been shown [182] that with the sensitivity of the forthcoming GAPS 
long-duration balloon flight experiment [184] a large fraction of a low-energy MSSM parameter 
space is accessible. 

1.3.2. The total electron-positron flux 

In contrast to stable nuclei, electrons and positrons lose their energy quickly through synchrotron 
radiation and inverse Compton scattering, see Sec. 1.2.6. They can therefore only travel finite, 
energy-dependent distances before losing all their energy. The fiux on Earth at the highest, 
0(1) TeV energies, is therefore dominated by the closest and youngest sources. In particular, 
considering that due to the discreteness of the sources in space and time there is/are necessarily 
a (few) source(s) of minimum distance and age, one predicts a propagation cut-off in energy. The 
older and further sources however add up to a smooth spectrum below ~ 100 GeV. For a more 
detailed and quantitative discussion see Sec. 2.3.1. 

The combined differential flux of GCR electrons and positrons, (Jg- + Je+), has been measured 
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by a number of experiments over the last decades, however, with considerable scatter, see Fig. 1.8. 
Here, we focus on two more recent, but somewhat contradictory results which both hint at some 
sort of excess with respect to pre- Fermi, standard GCR propagation models [126, 185]. 
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Figure 1.8.: Absolute electron plus positron flux J{E), scaled by E"^ (from [169]). The red data 
points with statistical error bars are from the Fermi-LAT measurement [169], the grey band de- 
notes the systematic uncertainty and the double headed arrow the energy scale uncertainty. Other 
measurements [148, 168, 167, 186, 187, 188, 189, 190] are referenced in the legend. The dashed line 
is from a GALPROP conventional diffusive model [191]. 

On the one hand, a sharp feature, i.e. a hardening from an spectrum at ~ 100 GeV and 
a rather strong shoulder at ~ 800 GeV, has been observed independently by PPB-BETS [167] and 
ATIC [168]. On the other hand, the Fermi collaboration has performed an analysis [169] with data 
from their Large Area Telescope (LAT) which shows a smooth continuation of the E~^ spectrum up 
to hundreds of GeV where the spectrum starts to soften. This is in agreement with a determination 
of the flux by HESS [186, 192] which also shows a softening to E~'^ around a TeV and is basically 
also inconsistent with the PPB-BETS/ATIC findings. Both signals, i.e. either a sharp feature or 
a broader excess, can in principle be interpreted as a contribution from DM annihilation or decay 
into leptonic channels in the galactic DM halo, cf. Sec. 1.3.1. 



49 



We note that an total electron-positron flux per se does not constitute an excess with respect 
to the standard picture of GCRs. For example a steady, homogeneous E~'^ '^ injection of electrons 
together with a diffusion coefficient oc E~^'^ could reproduce such a spectrum. The pre-Fermi 
models were adjusted to earlier measurements which indicated a softer spectrum. On the other 
hand, the slight dip and bump in the Fermi-LAT data at ~ 100 GeV and ~ 300 GeV, respectively, 
as well as a rather soft primary injection spectrum derived from, e.g. gamma-rays from SNRs (see 
2.4.1) can however be interpreted as indications for the presence of an additional component in the 
total electron-positron flux. 

Again presuming that the excesses seen in the total electron-positron flux are DM induced, we can 
further constrain possible annihilation channels and DM masses. Whereas the PPB-BETS/ATIC 
data favour annihilation into /U^^^ with moM ~ 1 Tev, the smoother Fermi-LAT spectrum is better 
fit by annihilation to r^r^ with moM ~ 2TeV [193]. For a harder background the Fermi-LAT 
data also allow the z^"*"^" channel [194]. Heavier DM can however not reproduce either excess and 
therefore, DM models leading to hadronic excesses (even beyond the current reach of antiproton 
measurements) are ruled out. 

1.3.3. Gamma-rays 

The major appeal of DM indirect detection through gamma-rays is that in contrast to charged 
particles, gamma-rays free-stream through the Galaxy and allow to be traced back to their sources. 
This allows, for example, to investigate directions in the Galaxy that are particularly apt for DM 
detection, e.g. due to low astrophysical backgrounds. Furthermore, gamma-rays reach the Earth 
without attenuation or energy losses, apart from redshift losses. (This is true at least up to energies 
of hundreds of TeV; above, the mean free path for gamma-rays sharply drops due to e"^ e~ pair 
production on the CMB down to lOkpc at a few PeV.) 

The physical processes possibly contributing to gamma-rays from DM annihilation or decay are 
prompt emission, i.e. internal bremsstrahlung (final state radiation and virtual internal bremsstrahlung) 
and neutral pion decay (from hadronic decay modes), or ICS of electrons and positrons, all of which 
have continuous spectra. The annihilation or decay to 77, 7Z or is, although loop-suppressed 
(DM is electrically neutral) and with a branching fraction of 10"^ or smaller, important because 
the resulting mono-energetic line emission is a very clean signature. 



50 



Again, we distinguish between DM annihilation and decay because of the different dependence 
on the DM density. As the annihilation rate depends on the DM density square, the differential 
flux of gamma-rays (p^^^ in a particular solid angle AO around a direction {i, b) in the sky contains 
the line of sight (l.o.s.) integral of the density square, p'^{r), 

Jr(^,b) = P^-^Y.^Bf f ^^f dsp\r{sj,b)), (1.68) 
^ 47r 2 mfj^ ^ dE^ Jao Ji.o.s. 

where (cTann^) is the thermal average annihilation rate, dNj/dEj is the differential spectrum for 
final state / and Bj is the branching ratio for gamma-rays. 

The decay rate on the other hand is proportional to the DM density, so the differential gamma-ray 
flux (p'^'^^ from DM decay reads, 

J^'^%i,b) = ^—Y,^Bj [ dn [ dsp{r{s,l,b)), (1.69) 
' Att moM ^ dE^ Ji.o.s. 

with r the DM decay rate. 

The recent data from the Large Area Telescope (LAT) on board the Fermi satellite [195] allow 
for various searches for possible DM signatures with a so far unprecedented accuracy. The LAT 
is a pair-conversion gamma-ray detector consisting of 4 x 4 towers of tungsten trackers on top of 
a electromagnetic calorimeters with a thickness of 7 radiation lengths and is surrounded by anti- 
coincidence detectors. It has a broad energy range of 20 MeV to 300 GeV and with its large field 
of view of 2.4 sr it covers the whole sky in two orbits, i.e. three hours. Its energy-resolution is 
around 10% and its point spread function decreases from a few degree at 20 MeV to less than 0.1° 
at hundreds of GeV. 

In the following, we briefly review the different targets for DM indirect detection in gamma-rays. 
Galactic centre 

The galactic centre is a prime target for DM searches since it contains the highest DM density and is 
therefore expected to be the brightest DM source in the sky. Of course, astrophysical backgrounds 
are also very bright close to the galactic centre, and so the prospects strongly depend on the 
form of the DM density profile in the inner kiloparsecs. On the one hand, density profiles fitted 
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Figure 1.9.: Significance for detection in hb (left panel) and r^r^ (right panel) in the ttidm — (cann?^) 
plane (from [197]). The region above the black line is excluded by EGRET, the region below the 
grey line is not detectable by Fermi-LAT and the region in between is detectable by Fermi-LAT if 
the 'conventional' [198] or 'optimised' [191] background model is realised. The shaded region can 
only be detected under assumption of the 'conventional' background [198]. 

to "observations" in A^-body simulations hint at a "cuspy" NFW profile [155, 196] with a rather 
steep, decline in the central kpc. On the other hand, kinematical observations show rather flat, 
iso-thermal profiles [171]. Again, the amplification due to the (P'ir) factor leads to much larger 
signals from DM annihilation than from decay. 

A sensitivity study [197], considering two different GALPROP diffuse backgrounds [198, 191], but 
no point sources, has been performed pre-flight. For generic DM annihilation into W~^W~ , bb, 
tt and r+r^, an NFW density profile and accounting for detector response, the sensitivity after 
5 years of data is good enough to dig deeper into the m-DM — {o'a.nnv) plane than possible with 
EGRET [199, 200]. In Fig. 1.9, we reproduce the sensitivity plots for gamma-rays for annihilation 
to bb and r^r~ from [197]. 

Shortly after the public release of the first 11 month data, an excess above astrophysical back- 
grounds has been claimed. A bump-like feature around ~ 2 GeV above power-law background 
has been identified [201] within the inner 3°. The excess has a steeper radial profile and is more 
spherically symmetric than astrophysical backgrounds. It could be explained by 25 — 30 GeV dark 
matter particle, annihilating to bb with ((Jann^) ~ 9 x lO^^^cm'^s^^ in a cusped density profile. 

The Fermi collaboration is of course performing their own studies of the emission around the 



52 



Galactic centre but has, so far, only presented very preliminary results. A binned likelihood analy- 
sis [202] with a simultaneous spatial and spectral fit has been performed within a region of interest 
(ROI) of 7° X 7° around the galactic centre, using only high-quality reconstructed events gathered 
during the first 11 month of operation. Various contributions have been considered, in particular, 
point sources have been modelled and subtracted, the diffuse galactic emission is accounted for 
by a GALPROP model and an isotropic, extra-galactic background is allowed for. The residuals, 
see Fig. 1.10, also show an unmodelled excess in the 2 — 5GeV range. The Fermi collaboration 
however concludes, that better modelling of the galactic diffuse emission and possibly unresolved 
point sources is necessary before any excess can be confirmed. 
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Figure 1.10.: Residual spectrum from the likelihood analysis of [202], integrated over the 7° x 7° 
ROI (from [202]). The blue band shows the systematic error from the uncertainty on the effective 



area. 



Galactic halo 

Although DM contribution to gamma-ray fluxes from the galactic halo are comparatively lower than 
from the galactic centre, the astrophysical backgrounds, in particular from (unresolved) sources, 
are also less strong. Still, diffuse backgrounds remain a formidable challenge for these studies. 
Furthermore, the spectrum from DM annihilation or decay is the same anywhere in the halo and 
is expected to have a particular shape with sharp cut-off at the DM mass (half the DM mass 
for decay) which is difficult to explain astrophysically. Usually, one excludes a region around the 
galactic centre or the galactic plane, thereby also limiting the influence of the uncertainty of the 
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halo profile in the inner kiloparsecs. 

The pre-flight study [197] has estimated the necessary (aannv) for a detection with 3 a sensitivity 
after only one year of data as a function of the DM mass. The uncertainty in the background 
model turns out to be large and the necessary ((Tann'f) increases from 7 x 10~^^ cm^ s~^ (for moM = 
50 GeV) to 6 X 10^^^ cm^ s^^ (for moM = 250 GeV) for the 'conventional' background model [198]. 
It was concluded that Fermi-LAT could probe a large region of MSSM or mSUGRA parameter 
spaces. 

Constraining the analysis to intermediate galactic latitudes 10° < |6| < 20°, the Fermi collab- 
oration has investigated [203] the diffusive background, maximising the fraction of background 
produced within a few kiloparsecs from the solar system and thereby minimising the dependence 
on the uncertainty of ISM and CR densities elsewhere in the Galaxy. The shape of the spectrum 
is consistent with an a priori diffusive emission model, an update of the 'conventional' GALPROP 
model [198]. The excess measured by EGRET [204] is also not reproduced and the deviation is 
therefore likely to be an instrumental effect of the EGRET experiment [205]. 

A conservative analysis that does not attempt to subtract or fit any astrophysical backgrounds 
is presented in [206]. Model- independent two-body annihilation or decay to leptonic and hadronic 
channels are considered and the resulting gamma-ray flux from final state radiation and vr'^-decay 
as well as from ICS of are required not too exceed the Fermi-LAT measurements by more than 
3(T in different, selected regions of the sky. This gives already quite considerable constraints in the 
"iDM ~ {o'ann'v) plane. Interestingly, a large fraction of the parameter space needed to explain the 
electron-positron excesses from DM, see Sees. 1.3.1 and 1.3.2, are ruled out. In particular, for NEW 
or Einasto profiles all fits to PAMELA, Fermi-LAT and HESS e^-data from DM annihilation are 
excluded, and even for a cored isothermal profile, only annihilations into muons remains marginally 
consistent. 

A different analysis [207] is based on template subtraction of foregrounds that are assumed to 
be traced by spatial templates. An excess, called the 'Fermi haze', is found up to |6| ~ 40° above 
the galactic centre which is satisfactorily fit by a bivariate Gaussian. This is argued to be most 
likely due to inverse-Compton scattering (ICS) by relativistic electrons, and that the underlying 
electron distribution is compatible with the 'WMAP haze' [207], see Sec. 1.3.4. While such a 
signature in ICS is naturally expected if there is indeed an additional population of electrons with 
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a hard spectrum, it was pointed out [208] that the template maps apphed in Ref. [207] are in 
fact inappropriate and underestimate both the vr*^ decay and ICS contributions to the gamma-ray 
emission, in particular in the galactic centre region. The 'Fermi haze' may therefore be an artefact 
due to incorrect foreground removal. Furthermore, with 1.6 years of data from Fermi-LAT, a recent 
analysis [209] finds this excess to be distributed above the galactic centre in an hourglass-shaped 
morphology. Both because of this elongation and the peculiar angular profile, the authors conclude 
that a DM explanation of this signal seems to be disfavoured. 

The Fermi collaboration has so far not confirmed this excess but is performing analyses with more 
sophisticated foreground models and is investigating possible correlations with (local) structures 
seen in radio maps, like the North polar spur, see [210]. 

Milky Way satellites 

The CDM paradigm and in particular A^-body simulations [13, 14, 15] predict a large num- 
ber of bound substructures. The masses of these so-called Milky Way satellites go down to 
(10~^ . . . 10~^^) Mq, depending on the free-streaming length of the particular DM model consid- 
ered. However, the minimum of the satellite masses that will be observable with Fermi is rather 
IO^Mq. 

A simple estimate of the significance expected from a satellite with truncated NFW profile and 
typical WIMP annihilation cross section sets the number of Milky Way satellites detectable by 
Fermi to 5fT within 5 years to ~ 12 [197]. Particular care must be taken not too misidentify 
statistical fluctuations in the gamma-ray flux as DM annihilation or decay from substructure but 
it was shown in [197] that a log-likelihood analysis can distinguish between both cases for a '5(t' 
satellite. Furthermore, the DM substructure needs to be distinguished from astrophysical sources. 

It is worth stressing that in the case of DM annihilation, substructures are not only targets in 
themselves but also increase the diffuse flux by the dependence (see Eq. 1.68). Estimates relying 
on the frequency of such structures seen in A^-body simulations however limits the boost factor to 
0(10) [211]. 
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Dwarf spheroidal galaxies 

Dwarf spheroidal galaxies (dSplis) are particularly faint companion galaxies of the Milky Way or 
Andromeda. They may well be the most abundant type of galaxies in the Universe but difficult to 
detect due to their faintness. Before 2005, there were only 9 known dSphs [212, 213, 214, 215, 216] 
but the Sloan Digital Sky Survey (SDSS) [217] has increased their number by 11, also improving 
our understanding of this type of galaxies. From stellar kinematics it can be inferred that some 
have mass-to-light ratios of up to 0(1000), that is many times more than in conventional types 
of galaxies. Furthermore, dSphs contain only little neutral or ionised gas which could otherwise 
contribute to its gamma-ray emission. Therefore, dSphs are the most extremely DM dominated 
environments known which makes them an interesting target for DM searches. 

A recent study [218] of 14 local group dSphs with data from the first eleven month of Fermi-LAT 
operation does not find any significant gamma-ray emission above 100 MeV and sets upper limit 
on their gamma-ray fluxes, both for power-law injection and spectra motivated by different WIMP 
models. Using stellar velocity data to model the DM content of a subset of 8 of the above dSphs, 
the Fermi team can constrain the cross section for annihilation into gamma-rays by bb, t'^t~' and 
a mixture of both, as motivated from neutralino DM. More precisely, the limits are starting to be 
competitive with cross sections from a scan over mSUGRA parameter space and already exclude a 
fraction of the MSSM parameter space. A model of UED with S^^^ KK DM is not constrained but 
SUSY models with Wino-like DM in the context of anomaly mediated SUSY breaking (AMSB), 
see e.g. [178], can be ruled out for ttidm ^ 300 GeV. Interestingly, models invoked to explain the 
cosmic ray lepton excesses (see Sees. 1.3.1 and 1.3.2), can be constrained in so far as masses above 
~ 1 TeV can be ruled out by considering the ICS emission produced. 

Line emission 

Although DM is found to be electrically neutral and annihilation or decay to photons can therefore 
only occur at loop order, the preferred smoking-gun signature of DM in indirect detection is the 
gamma-ray line. The internal width of the line is for most annihilation (excluding, perhaps anni- 
hilation via a Z boson resonance) and decay (due to the necessarily long lifetime) processes rather 
small. The width of the line, ~ 10^^, is therefore mostly due to Doppler broadening although this 
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effect is also small because of the non-relativistic WIMP velocities today. The lack of astrophysi- 
cal sources with lines in the GeV - TeV regime finally makes this signature virtually background 
free. Theoretical frameworks include SUSY WIMP and gravitino DM models. Unfortunately, the 
branching fraction of the necessary processes is found to be 10"'^ or smaller. 

A recent Fermi study [219] has found no significant excess and has set limits on gamma-ray lines 
between 30 and 200 GeV. By considering three different halo profiles, NFW [196], Einasto [220, 221] 
and isothermal [171], upper limits on the annihilation and decay cross section have been derived. 
Although for the annihilation cross sections the bound is about an order of magnitude weaker than 
the cross section for thermally produced WIMPs, some non-thermal production models, e.g. [178], 
can again be ruled out. The upper limits on the life-times are also constraining for some models 
with monoenergetic lines from decay of gravitinos. 

Other targets 

Other targets include so-called cosmological DM [222], that is extra-galactic haloes and large- 
scale structure which contribute to the isotropic gamma-ray flux, and DM signals from clusters of 
galaxies [223]. 

1.3.4. Radio and microwaves 

The measurement of anisotropics in the cosmic microwave background (CMB) by COBE and 
WMAP has ushered in an exciting new era in cosmology. The study of the cosmic signal requires 
careful subtraction of galactic foreground emissions and this will become even more crucial for stud- 
ies of the 'B-mode' polarisation signal by PLANCK and the proposed CMBPol satellites [224]. A 
by-product of this foreground subtraction is the study of diffuse galactic microwave emission in its 
own right. (In fact most of the diffuse foreground emission is of galactic origin, a counter-example 
being unresolved extra- galactic radio sources.) Usually, one distinguishes three different physical 
components which have different underlying processes and frequency behaviour: 

1. Free- free emission 

Thermal electrons produce bremsstrahlung on the gas in the ISM. The intensity produced 
by free- free emission, Ii, oc z^"*^'^^, is proportional to the line of sight integral of the ISM gas 
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density squared which is why it is beheved to be traced by recombination Ha hne maps [225]. 
Dust can absorb some of this emission and usually regions close to the galactic plane with 
dust optical depth too high are excised. 

2. Dust-correlated emission 

Dust grains vibrating in equilibrium with the surrounding radiation fields are producing 
emission in the microwave range. This 'thermal dust' emission has been mapped [226] and 
evaluated at 94 GHz [227]. In addition, larger dust grains can have an electric dipole moment 
and be excited by collisions with ions [228]. While individual dust grains radiate with a bump 
like spectrum, the superposition of many dust grains with different frequencies can generate 
a soft power law spectrum. To a first approximation, this 'spinning dust' is expected to 
spatially correlated with thermal dust and in fact an anomalous, soft component was found 
in Ref. [229]. Although the presence of an anomalous component is widely accepted, its 
interpretation as spinning dust is still being debated, for example by the WMAP team itself 
who do not conclusively confirm the spinning dust hypothesis. 

3. Galactic synchrotron 

GCR electrons and positrons produce synchrotron radiation on the galactic magnetic field 
(GMF). For a power law spectrum of relativistic electrons, ngi oc E~", the synchrotron 
emissivity is also a power law, 1^ oc z^^^^")/^. For an electron spectral index a ^ 3, 1^ oc 
which is close to the observed i>~^'^ [230]. Tracer maps include radio maps or difference maps 
at microwave frequencies, used, for example by the WMAP collaboration [231]. 

For dark matter indirect searches more important is, however, the claimed presence [232] of a new 
foreground that does not correlate spatially with any of the above processes. Using the template 
fitting technique (see 4.2), a residual remains after subtraction of the known foregrounds - the 
so-called 'WAMP haze'. This residual has a roughly spherical morphology localised around the 
centre of the Galaxy, see Fig. 1.3.4, and a harder spectrum [233] than synchrotron radiation by 
relativistic GCR electrons from standard astrophysical sources like SNRs. An independent analysis 
has confirmed the existence of the haze [234], but others [235, 236] do not find the evidence to be 
significant, including the analysis [237] by the WMAP collaboration itself. 
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Figure 1.11.: Sky map of the 'WMAP haze' (from [233]). 



Initially, it was believed that the haze is free-free emission from ionised gas too hot to be traced 
by recombination line maps but too cold to be visible in X-rays [232]. However it was suggested 
later that it is in fact synchrotron emission from a new population of relativistic electrons and 
positrons, produced by dark matter annihilation [238]. It is indeed possible to explain the haze [239] 
by the synchrotron radiation of an additional population of relativistic electrons and positrons, 
produced by WIMP annihilation with a thermal cross section. Other authors however argue that 
the annihilation cross-section needs to be significantly boosted over the usual value [236] or that 
the parameters of the diffusion model used are somewhat non-standard [240]. There have also been 
attempts to fit both the morphology and spectrum of the haze by ascribing it to electrons emitted 
by pulsars with a hard spectrum [241, 242]; however the expected haze is then less spherical since 
most pulsars are in the galactic disk. This is also true of SNRs which have in fact recently been 
invoked [243, 2] as sources of positrons with a hard spectrum to explain the rise in the cosmic ray 
positron fraction at high energies measured by PAMELA [163]. 

More conservative constraints on DM annihilation or decay can be derived by demanding that 
the possible contribution from DM does not exceed the observed fluxes in radio and microwaves 
(including astrophysical backgrounds). For examples of such studies, see [244, 245]. 
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1.3.5. Neutrinos 

Dark matter indirect detection with neutrinos benefits from the same advantages as indirect detec- 
tion with gamma rays: Neutrinos propagate through the Universe truly unimpeded and also point 
back to their sources. Their small interactions however require detector volumes of cubic kilometre 
size and provide event rates of a few per year, both for conventional astrophysical application as 
well as for DM searches. The weak interaction cross section are however also an advantage as they 
could possibly make neutrinos the only messengers that can escape very dense environments. 

One such situation can lead to a particularly interesting idea for DM indirect detection. The 
Sun is constantly moving through the DM halo. Although the scattering cross section of WIMPs 
off the nucleons in the Sun is small, occasionally, a WIMP can scatter elastically and become 
gravitationally trapped inside the sun. Over the lifetime of the Sun, a large number of WIMPs 
can become trapped and an equilibrium between WIMP capture and annihilation will be reached. 
Neutrinos detected from the Sun with energies of tens to hundreds of GeV can be unambiguously 
assigned to these annihilation processes and used to constrain WIMP annihilation cross sections 
and DM models. 

Such neutrinos could be detected by the SUPER-Kamiokande experiment which has the best 
sensitivity for WIMP masses of about a GeV, whereas for larger masses the neutrinos could be 
detected in IceCube [246] (in particular in its low-energy extension, "DeepCore"), currently under 
construction under the South pole. IceCube consists of a cubic kilometre of clean ice, about 2000 km 
under ground and instrumented with thousands of photo multiplier tubes. A muon converted from 
a muon neutrino in the surrounding material or the detector itself leaves a track of Cerenkov light 
in the detector which allows the reconstruction of the neutrino direction and energy. 

The beauty of detection of neutrinos from WIMP annihilation inside the Sun is its relative robust- 
ness against astrophysical uncertainties; the flux depends only the local DM density and velocity 
distribution which also affect DM direct searches since in equilibrium the neutrino production only 
depends on the capture rate. The rate has contributions from both spin-independent (SI) and 
spin-dependent (SD) cross sections and it turns out that for the SI cross section, upper limits from 
direct detection experiments basically rule out any prospects for event rates 0(1) yr~^. 

SD cross sections, on the other hand, can become much larger than SI ones, e.g. for certain 
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regions of the MSSM parameter space that are experimentany not too weh constrained and could 
lead to event rates of up to 1000 yr~^. Large SD cross sections for neutralino WIMPs are usuaUy 
associated with large Z couplings and hence a large Higgsino mixing. Therefore, the focus point 
region of the mo - ?Tiy2 plane is particularly promising. 

IceCube has already provided limits on the elastic scattering cross section better by 2-3 orders 
of magnitude for the spin-dependent case [247]. For spin-independent cross sections, however, the 
event rates allowed by direct detection limits are much too low to be discovered by IceCube. 

Of course, similar to gamma-ray searches, neutrinos can be searched for from other targets, for 
example the galactic centre [248, 193]. For preliminary result from IceCube, see [249]. 

1.4. Conclusion 

Indirect detection of DM, that is the search for deviations in the fluxes of cosmic ray nuclei and 
charged leptons as well as radio/microwaves, gamma-rays and neutrinos, is a very promising idea. 
Not only capitalises it on the large amount of dark matter contained in the Milky Way halo and 
possibly also in extra-galactic haloes, but it also uses a huge variety of astrophysical environments 
or targets, ranging from substructures of the halo of a few solar masses to large galaxy clusters. 

Astrophysical backgrounds, however, play a crucial role for the prospect of detecting DM in these 
channels. In particular, with most DM explanations of the positron and electron-positron excesses 
starting to be ruled out by gamma-ray constraints, for example, it is clear that these signatures 
must be of astrophysical origin. Therefore, a better understanding of the astrophysical backgrounds 
is crucial both to resolve these problems in our understanding of current excesses but also for future 
indirect searches. 

Even when DM indirect detections provide a (first) evidence of DM annihilation or decay, pro- 
vided it had some non-gravitational interactions, it is unlikely that we could identify the interesting 
DM parameters, mass, couplings etc., let alone the full DM model. In most cases, the uncertain- 
ties inherent to the DM modelling itself as well as uncertainties induced by the astrophysics of 
its detection are considerable and all we can hope for are some rough indications. DM indirect 
detections can therefore only unfold its true potential in combination with other approaches to DM 
detection. Direct detection (for a recent review, see, e.g. [250]) and accelerator searches (e.g. [251]) 
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are two other, equally important approaches that, unfortunately, could not be discussed in this 
introduction due to page limitations. Other possible constraints can come from other effects that 
DM would have, for example, on cosmological observables: DM with large annihilation rates or 
short lifetimes could also affect big bang nucleosynthesis [252, 253, 254, 255, 256]. One of the most 
stringent constraints however comes from the effect of DM annihilation or decay on the CMB and 
during the epoch of reionisation [257, 258, 259, 260]. 

For the phenomenology of indirect searches, what is needed at the moment are a better un- 
derstanding of the backgrounds that go beyond the (often) oversimplified assumptions and hence 
predictions of the standard picture of GCRs. In the rest of this work we present our humble 
contribution to this ongoing challenge. 
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2. Additional Electrons/Positrons from 
Astrophysical Sources 



2.1. Introduction 

The excess in the positron fraction recently measured by the PAMELA collaboration clearly hints 
at some additional positron component. The astrophysical background is expected to be falling at 
energies above ~ 1 GeV if positrons are indeed mainly produced as secondaries by decay of charged 
pions from the spallation of cosmic ray protons and nuclei on the interstellar medium (ISM). In 
addition, recent measurements of the total electron-positron flux by the Fermi collaboration have 
shown it to be harder than previously expected. 

The evaluation of the astrophysical predictions has gained in importance due to the possibility of 
these excesses being due to exotic physics, in particular due to additional electrons and positrons 
produced in annihilation or decay of dark matter (DM) in the galactic halo. A large number of 
these models is starting to be ruled out by constraints from other messenger, e.g. radio waves 
from synchrotron radiation or gamma-rays from inverse Compton scattering (ICS). It is therefore 
of utmost importance to investigate how these excesses can be understood in terms of astrophysical 
sources. 

One of the most minimal approaches is to consider whether the additional electrons/positrons 
could possibly be produced in the standard setup without making any additional assumptions 
about new sources, new interactions or the like. Such minimality is a basic principle of scientific 
heuristics as most famously expressed in 'Occam's razor': 'Numquam ponenda est pluralitas sine 
necessitate' [261]. (Plurality ought never be posited without necessity.) or 'Frustra fit per plura 
quod potest fieri per pauciora' [262]. (It is futile to do with more things that which can be done 
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with fewer.) Instead, the assumptions of the current model of galactic cosmic rays (GCRs) are to 
be reconsidered, possibly premature arguments about the importance of different effects are to be 
revisited and the available parameter space is to be reassessed. 

In particular, we are going to question the assumption that positrons get exclusively produced 
by spallation of CR protons and nuclei on the ISM and, instead, we consider the production of 
secondary electrons and positrons by protons and nuclei inside the paradigm sources: supernova 
remnants (SNRs). Not only can the total fluxes of these secondaries reach levels comparable to the 
primaries but the secondary spectra will also differ quite drastically from the primary ones. Hereby, 
we follow a recent idea to explain the rise in the positron fraction by such a harder injection of 
secondaries from SNRs [243], however improving on this model in several respects. 

The effect of secondaries produced inside the GCR sources had so far always been neglected due 
to a crude column depth argument, recently reiterated [263]: The average residence time of primary 
nuclei in a SNR is necessarily shorter than the SNR lifetime, 0(10^) yr, and therefore happens to 
be much smaller than the residence time in the ISM, ~ 0(10) Myr. Even if the gas density inside 
the SNR is higher than in the ISM, due to compression by the shock, the grammage experienced 
inside the SNR, 0.2 (ngas,SNR/(10 cm-^)) 

gem is still much smaller than the average grammage 
of the ISM, a fewgcm"^. Therefore, the production of secondaries inside SNRs was believed to be 
negligible. 

What this argument ignores is that the secondaries produced inside the source get also accelerated 
and their spectrum becomes considerably harder than that of the secondaries produced during 
propagation of the primaries in the ISM. As we will see below, the accelerated secondary positrons 
in sources have a power law index of —7 + 1 where —7 is the index of the primaries at source. After 
propagation this gives a spectral index —7 at energies above 10 GeV where cooling is the dominant 
loss process. The secondary positrons from propagation, however, have a spectral index of —7— (5—1. 
We remind ourselves that the value predicted for 7 by diffusive shock acceleration (DSA) is 7 ~ 2 
and 6 is determined from nuclear secondary to primary ratios as 5 ~ 0.6. We now compare the fluxes 
«^rc = Jsrci^ / ^o)~'^ of Secondaries from the source and Jism = Jisui^ / -^0)^''^^^^ °f secondaries 
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from propagation. Although the ratio of the total number of particles is indeed small, 

^'^^ = -^^0 ~ ^'^^ ' ~ 2 6 X 10^ (2 1) 

iViSM di?JisM(i?) 4m7-1 

if we chose, for example Jsvc/ ^isM = 2 x 10~^ and E'max ^ Eq, the differential fluxes become 
comparable at {Ji^mI '^srcY^''^^^^ — SOE'o already. 

The appeal of the acceleration of secondaries, besides its minimality, is that its effect is guaran- 
teed, in the sense that secondaries will definitely be produced in cosmic ray source and are subject 
to subsequent acceleration. However, the normalisation of the flux of these secondaries, for example, 
with respect to the secondaries from propagation is the crucial question. 

In this chapter, we will study the acceleration of secondaries in sources in some detail. In Sec. 2.2 
we will calculate the spectra of the secondaries in the simple test particle approximation of DSA 
and explain why the secondary spectrum is harder than the primary one. We will then consider the 
effect of discreteness of sources and suggest a realistic distributions of sources in the spiral arms 
of the galactic disk and how to determine it from cosmic ray spectra alone in Sec. 2.3. Sec. 2.4 
contains a summary of the different populations of electrons and positrons and some comments on 
how to normalise them by fitting to data. Previously, the flux of secondary e~ and e"*" in the sources 
has been normalised with respect to the primary electrons in an ad hoc fashion [243]. Instead, we 
exploit the hadronic origin of these secondaries and normalise using the gamma-ray fluxes (assumed 
to be from vr'^ decay) detected from known SNRs by Imaging Air Cerenkov Telescopes (lACTs) like 
HESS. We can thus fix the only free model parameter by fitting the total e~ + fiux to Fermi- 
LAT and HESS data. The fraction is then predicted up to TeV energies and provides a good 
match to PAMELA data (Sec. 2.5). Having constrained the distribution of the closest SNRs via 
the measured e~ and e"*" spectra, we present in Sec. 2.6 an example of a likely source distribution in 
order to illustrate that there are good prospects for IceCube to detect neutrinos from nearby SNR. 
In Sec. 2.7, we counter some criticism on this model, in particular the claim that the acceleration 
of mechanism does not hold in a time-dependent treatment. A consistent picture thus emerges 
for all presently available data in the framework of the standard DSA/SNR origin model of GCR. 
Some alternative, astrophysical explanations of the excesses in GCR leptons are being mentioned 
in Sec. 2.8. We conclude and comment on some open issues and grounds for concern in Sec. 2.9. 
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2.2. Acceleration of Secondaries in the Sources 



Similar to the diffusive shock acceleration of primary cosmic rays, we describe the acceleration of 
secondaries produced and accelerated in the source in the simple test particle approximation, that 
is ignoring the back reaction of the energetically important nuclear component of the cosmic rays 
on the structure of the shock front. 

The setup is similar to the one discussed in Sec. 1.2.5. We consider the phase space density, f±, 
of secondary e~ and e"^ produced by the primary GCR, both undergoing DSA, which is described 
by the steady state transport equation, 

- — (d — / ^ +-—p^ +q (2 2) 

dx dx \ dx J 3 dx dp ' 

where the source term q± is determined by solving an analogous equation for the primary GCR 
protons (see Sec. 1.2.5). (Ideally we should solve the time-dependent equation, however we do 
not know the time-dependence of the parameters and can extract only their effective values from 
observations. This ought to be a good approximation for calculating ratios of secondaries to pri- 
maries from a large number of sources which are in different stages of evolution.) We consider the 
usual setup in the rest-frame of the shock front (at x = 0) where ui (U2) and ni (77-2) denote the 
upstream (downstream) plasma velocity and density, respectively (cf. left panel of Fig. 1.2). The 
compression ratio of the shock r = ui/u2 = n2/ni determines the spectral index, a = 3r/(r — 1), 
of the GCR primaries in momentum space. 

For X ^ 0, Eq. 2.2 reduces to an ordinary differential equation in x that is easily solved taking 
into account the spatial dependence of the source term, cf. Eq. 1.28, 



( )gX«i/Z3{pp) forx<0, 

(2.3) 

Q±,2{P) forx>0. 



where the proton momentum pp should be distinguished from the (smaller) momentum p of the 
produced secondaries, the two being related through the inelasticity of production: ^ ~ 1/20. 
Assuming D oc p (Bohm diffusion) in the SNR, the solution to the transport equation 2.2 across 
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the shock, satisfying the boundary conditions, 



Um f±=0, hm 

x—^—oo x—^—oo ox 



and 



hm /± 

a;— >oo 



< OO , 



(2.4) 



can then be written: 



f±={ 



1±^l , f g=c/'li_g^x/di 



/,0 



for X < 0, 
for x > 0, 



(2.5) 



where di = D/ui is the effective size of the region where e~ and e"*" participate in DSA (see right 
panel of Fig. 1.2). 

Continuity at the shock front x = requires: 



D 



dfi 



dx 



D 



dfi 



x=0- 



dx 



1 



{U2 - Ui)p 



x=0+ 



djl 

dp ' 



(2.6) 



yielding the differential equation. 



P- 



dp 



(2.7) 



This is readily integrated with boundary condition /^(O) = and gives 



p 



(2i 



Assuming Feynman scaling for the pp interaction, i.e. pdcjpp/dp oc T,± we can express the momen- 
tum dependence of the source term as 



q±Ap) = dpNcRip 



~ — ^A^cr(p): 



dp 



4-Kp^ 



7-2 



(2.9) 



We can easily interpret the solution Eq. 2.5 in terms of power laws in momentum. The second 
term downstream, (^2 7^2)2;, follows the spectrum of the primary GCRs (oc p~") and describes the 
production of secondary e~ and e"*" that are then advected away from the shock front. However, 
secondaries that are produced within a distance ^ D/u from the shock front are subject to DSA 
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(see Eq. 2.3 and right panel of Fig. 1.2). The fraction of secondaries that enters the acceleration 
process is thus given by the ratio of the relevant volumes, i.e. {D / ui) / {u2 tst<}B.) , and the number 
density injected into the acceleration process is (l/e + r2)Dg±,i/uf. This rises with energy because 
of the momentum dependence of the diffusion coefficient {D{p) oc p) so the first term downstream 
in Eq. 2.5 gets harder: f±{p) oc p~^'^^. 

The total injection spectrum R± from one SNR is obtained by integrating the steady state 
solution over the volume of the SNR: 



R^=4ttp'^Att / dx x'^ f±{x , p) . 

Jo 

The resulting source spectrum, R±, is thus the sum of two power laws, 
where the "cross-over" nioinentuni, Pcross? 

satisfies 



(2.10) 



P 



(2.11) 



4 a{l/t, + r^) 

As has been noted [243], this mechanism is most efficient for old SNRs where field amplification 
by the shock wave is not very effective any more. We therefore introduce a fudge factor K-q 
that parametrises the effect of the smaller field amplification on the otherwise Bohm-like diffusion 
coefficient in the SNR, 

3.3 X 10- A-B(^)"'(^)cmV'. (2.13) 

The number of particles entering the acceleration process can of course not exceed the total 
number of secondaries produced inside the SNR. This effectively caps the growth of the term 
D{p')q± i{p')/ui once {D / ui) / {u2 Tst^R) becomes larger than unity, a relation that defines a char- 
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acteristic momentum scale Pbreak- We therefore substitute in Eq. 2.8 



D{p)qIi{p) 
ul 



D{p)qI,i(p) 



for p < Pbreak, 



— ^ < 



(2.14) 



for p > Pbreak- 



The source spectrum R± thus returns to a p ^ dependence around p = Pbrcak- At even higher 
energies the secondary spectrum cuts off at the same E'cut as for primary electrons (see Sec. 2.4.1). 

2.3. The Discreteness of Sources 

In the standard picture, see Sec. 1.2, the explanation for the observed galactic cosmic rays (GCRs) 
is factorised into two parts which can be treated separately: the acceleration in supernova remnants 
(SNRs) and the diffusive-convective transport through the Galaxy. 

Although SNRs can become quite large (up to 0(100) pc towards the end of their lifetime, see 
Sec. 1.2.3 and Eq. 1.15), they are under most circumstances still much smaller than the scale 
on which the propagation of GCRs take place, i.e. kiloparsecs, and their extent can usually be 
neglected for all practical purposes. 

We are therefore not dealing with a continuous distribution of sources within the galactic disk 
but with a large number of discrete sources. The usual assumption is that diffusion will wash out 
this source distribution. This is however only true as long as the diffusion length is long compared 
to the average distance of sources (which is of the same order as the average distance of an observer 
to the closest source). For protons and (stable) nuclei the diffusion length is indeed quite long 
- of the order of tens of kpc - since energy losses are negligible above a few GeV and therefore 
escape from the cosmic ray halo is the dominant loss process. For charged leptons on the other 
hand, in particular for electrons and positrons, the diffusion length is strongly energy dependent. 
The dominating energy losses are through synchrotron radiation and ICS and the diffusion-loss 
length, X{E) = y^P{E, 2E) oc (see Eq. 1.48), is quickly declining. Therefore we expect to 

start seeing contributions from individual, nearby sources above the energy at which the diffusion 
length drops below the average source distance. Depending on the particular values of the diffusion 
parameters, one predicts "bumps" , that is Green's functions of individual sources, to appear above 
(50 . . . 200) GeV. The superposition of a few such bumps leads to features which we will investigate 
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in the following section. 

In addition to the discreteness in space, GCR sources like SNRs also have a finite lifetime. A 
common assumption is that the bulk of cosmic rays will only be released after the SNR has entered 
the radiative stage and DSA comes to an end.^ This makes the distribution of sources also finite 
in the time-domain and will lead to even stronger features. 

2.3.1. Fluctuations in the electron flux 

Apart from some inferences about (only the observed) near SNRs we do not know much about the 
position of sources of GCRs. Even in the limit of a continuous source distribution, which is justified 
when considering old, far away sources and hence fluxes at lower energies, our knowledge of the 
global source distribution is limited, see also Chapter 4. We have however argued that at higher 
energies, the contribution to the electron-positron flux from particular sources depend strongly 
on their distances. Our ignorance of the source positions therefore introduces an uncertainty into 
our prediction for the electron-positron fluxes at Earth. Put differently, the electron-positron flux 
stochastically fluctuates between different realisations of a statistically smooth source distribution. 

In the following we will use central limit theorems and the analytical form of the Green's function 
for the diffusion-energy loss problem described in sec. 1.2.6 to infer the statistical properties of the 
flux, i.e. its average, its quantiles and the energy dependence of both. For numerical studies of this 
effect based on a Monte Carlo approach, see [264, 265, 266]. 

The flux from a source that injected a spectrum Q{E) = QqE^^ of electrons or positrons a time 
t ago at a distance L from the observer is given by the Green's function Gdisk (cf- Eq. 1.51) of the 
diffusion equation, 

J,{E) = -^Gdisk = f K)"' e-^'/^'goi?-^(l - boEty-'-x (-) , (2.15) 

47r 47r Zcr \ ^cr / 

where c/(47r) denotes the "flux factor" for relativistic particles and £'^{E,t), xi^/zcr) and Zcr are 
defined in Sec. 1.2.6. 



In more sophisticated modelling, one expects that a certain fraction of the particles accelerated will be able to 
escape the SNR during the Sedov- Taylor phase already [71, 72]. For example, due to a lack of turbulence at low 
wave numbers ahead of the shock, high energy particles cannot be confined to the shock region and are instead 
injected into the interstellar medium. 
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The flux J of identical^ sources at distances {Li} and times {tj} is the sum of the individual 
fluxes, 



J = J2 ME) = ^Y. K)"' ^-"-^""QoE-^I - boEur-'-x - . (2.16) 

To investigate the fluctuations introduced by the stochasticity of source distances and times, we 
calculate the expectation value and standard deviation of the total flux J exploiting the central 
limit theorem [267]. 

If the central limit theorem was applicable, at a fixed energy E the fluxes J{E) for different 
realisations of the same source density would follow a normal distribution with mean and 
standard deviation aj, 

f.j =^N^z = ^N{Z) , (2.17) 
aj =^VNaz = -^VNy^{Z^)-{Z)\ (2.18) 

were (Z™) denotes the moments of the Green's function Z = Gdisk{E, L,t), 

(Z-)= J dzfz{z)z"'. (2.19) 

As a function of the random variables L and t, Z itself is a random variable with the probability 
density fz- In the case under consideration, L and t are assumed to be independent random 
variables with probability densities /l and ft, respectively, and thus the joint probability density 
fL,t factorises, fL,t{L,t) = fL{L)ft{t). 

Let's assume that the sources are, although discrete, homogeneously distributed in a ring around 
the observer with an inner and outer radius ri and r2, respectively (see Fig. 2.1). The distance L 



^Here and in the following we assume that all the sources considered have a common injection spectrum Q{E). Of 
course, in reality the sources differ in total power output and spectrum which will introduce additional fluctuations 
into the spectrum observed at Earth. Here, we however constrain ourselves to investigate the effects of the 
discreteness of sources. 
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Figure 2.1.: Propagation setup. The sources are discrete, but homogeneously distributed around 
the observer in a ring of inner radius ri and outer radius r2- The disk is contained in a diffusion 
volume of half-height ^max- We neglect the boundary condition in radial direction. 

to some source is a random variable with the probability density, 

2L/(r2 — rf) for ri < L < r2 , 



ML) = { 



(2.20) 

otherwise . 



Also, we assume that the source rate, e.g. the supernova rate, is constant over the time scales 
considered, 0(100) Myr (electrons and positrons of GeV energies can diffuse over 0(100) Myr before 
losing their energy but are, of course, subject to escape losses). Therefore, we take the sources to 
be equally distributed up to a maximum time tmax = 1/ (^o-E'min), set by the minimum energy -Bmin 
which is to be considered, 



Mt) = { 



l/imax for < t < tmax , 

(2.21) 

otherwise . 



Rewriting the m-th moment of the Green's function as an integral over L and t (see Appendix A), 
we have, 

= J dzfz{z)z^ (2.22) 

r2 /"imax 

dL dtfL,t{L,t)G"'{L,t) (2.23) 
n Jo 

fL{L)dL dtft{t)G"\L,t) (2.24) 



ri JO 
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dt / dL- ^_ 

Jri "^2 '1 



X 



^cr \ '^cr 



max JQ Jrf '2 '1 \ cr \' 

1 1 /"^max 

max '"l -^0 



(2.25) 
(2.26) 



{rrf) QoE-^l-boEt) 



We now substitute, 



E, 



E 



dt 



dE, 











1 - boEt (1 - boEt)'^ boE"^ ' 

to obtain an integral in energy Eq at source, 

dEo 



= ^— ^ / 

Tmax ^1 J E 

















m 













with 



60(1 -<^) 



For m = 1, this is. 



(2.27) 



(2.28) 



, (2.29) 



(2.30) 



imax t\ - r\ Je boE^ 



-X 



''cr \ ^cr 



(2.31) 



which can, for example, be evaluated numerically. The expectation value (^J = c/{4:Tt)N{Z) for the 
flux is identical to the flux from a continuous, homogeneous and steady source distribution. 



^{^2 ''l) ^max 



(2.32) 



in a ring with inner and outer radius ri and r2, respectively. The first term is the surface density 
of sources and the second one the power spectrum per source (although we are, strictly speaking, 
comparing to a continuous source distribution). 
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Expanding the term xi^/^cr) in Eq. 2.29 (associated with the diffusion in the z direction) for 
^ *C Zcr, which is only justified for energies above ~ 50GeV for Zmax — 4kpc and amounts to 
ignoring the boundary (condition) in the z direction, 



^cr V Zcr J V 7r£*^ 



(2.33) 



we can get analytical estimates for the moments (Z™'). To this end we further substitute, 



4Z)o 

Eo = E{l- A^)^ 



s-i 



di?o = ^(l-A2)^-^dA2 



(2.34) 



similarly define. 



4.D0 



4Z)o 



(2.35) 



and find. 



1 1 



(4Do 



X rdA2(l-A2)=^(AY-f- 
Jo 



^-2+5+|m(l-(5)-m7 



Pi 



(2.36) 



Eventually, we want to send ri to zero since the distance to the nearest source is not physically 
limited. At the energies for which the above expansion holds, we also have i'^ *C r^, and therefore 
the exponential term in the integrand is w 1. For m = 1, the integral converges. 



dA2(i-A2)^^-i(A2)-V2 



^ ' 2(1-5) 



(2.37) 



where F is the gamma function and the average flux is. 



c c 1 N Qo 

IJ-J = fJ-z = -; . — o 

47r 47r y^4Do6o(l - 5) vrr| i 

max 



^-7-l+(l-5)/2_ 



7-l\ 
1-5) 



^ ' 2(1-5) 



(2.38) 
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With the parameters chosen as shown in Table 2.1 this gives, 

E^Hj ~ 150 GeV^i cm-2 g-^ sr^^ , (2.39) 

close to the featureless E~'^ spectrum measured with Fermi-LAT [169]. 

For m = 2, the integral in Eq. 2.36 however diverges; the second moment, {Z'^), and hence the 
variance cr^ is infinite. This could be cured by reinstating exp (— 2p^/A^) with a finite pi as a cut- 
off in the integral for small A'^, however, for the price of an additional and, even worse, physically 
unmotivated parameter. For ri = 0, the standard deviation of the flux therefore does not exist and 
cannot be used as an indicator of the amplitude of the stochastic fluctuations. Furthermore, the 
usual central limit theorem cannot be applied any more. 

The reason that the expectation value is finite but not the variance is that the probability density 
fz{z) has a broad power-law tail. For such cases, a generalised central limit theorem [268] is 
applicable: The centred and normalised sum of N independent and identically distributed (iid) 
random variables Zj converges against a stable distribution 5(q,/3, 1,0, 1) [269] if the probability 
density fz{z) behaves like |2;|-°-i for z — )• oo. In general, the distribution function for S is not 
known analytically but can be calculated as the inverse Fourier transform of its characteristic 
function. To determine the parameters a and /3, we need to find the asymptotic behaviour of the 
probability density fz{z) for large z. 



Table 2.1.: Summary of parameters used in the Monte Carlo simulation. 



Diffusion Model 


Dq lO^^cm^s-i 




5 0.6 


from GCR nuclear secondary-to-primary ratios 


2^max 3 kpc 




ho 10-16 GeV-i s-i 


CMB, IBL and B energy densities 


Source Distribution 


tmax 3 X 10« yr 


from £^min ~ 1 GeV 


N 1 X 10^ 


from supernova rate and tmax 


Source Model 


go 8.4 X lO*'^ GeV-i 
7 2.2 


1 fit to absolute + e flux 


r2 15 kpc 
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The distribution function Fz{z) satisfies, 



Fz{z) = dtdLmfUL) 



(2.40) 



with the region T>z given by the condition Z < z. With Z = G{t, L, E) this can be transformed 
into a condition on L: 



Z<z ^ L^> L^i„ = -f log z + f log ((7r^2)-3/2g^^-7(i _ 6o^t)T-2) , (2.41) 



and hence 



1 1 rimax rr2 

dt / dL2L 

J max[0;Lmin] 



max ' 2 JO 

max ?'2 \ Jo J L'^ . Ju JO 



(2.42) 

(2.43) 



where Lmin > for < t < and Lmin < for t > t^. The probability density can then be 
obtained by differentiating, 



dz tmax Jo V 



(2.44) 



With the substitutions, Eqs. 2.28 and 2.34, this reads. 



fz{z) 



1 1 ADo 



E^-^ / * dA2 (1 - A2) 



-<5/(5-l) 



z 



(2.45) 



with A^ defined by \\2^x2 = where 



^min = -^^\ log log 



47rL>o 
60(1 - S) 



3 

'2 3 



+ ^logA2 + ^log(l-A2) 



(2.46) 



As we are looking for the asymptotic behaviour of Fz{z) for large z which we anticipate corresponds 
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to ^ 1, we can neglect the (1 — A^) term in A^j^^ and solve A^j^^ for A^, 



Xl :^ . (2.47) 

AttDq 



With the same approximation, the integral in Eq. 2.45 gives (A^) /2 and 



M^) - t^^^E-^-'^-'QT^-'^' ■ (2.48) 



The asymptotic behaviour of the distribution function Fz{z) is therefore 

1 - ~ c+z-/^ with = l^^AryE-'-^'Qt ^ (2-49) 

^ '■max ' 2 " 

and = for z < 0. 

The generalised central limit theorem [270] then states that the centred and normalised sum Xjv, 

= ^ [Y^Xi-ai^ , (2.50) 

weakly converges to the stable distribution 5(a,/3, 1,0, 1) with a = 4/3 and (3 = 1. The normali- 
sation constants are 

aN = Nnz and = i — tt^^^— y-^l N''I\ (2.51) 

1^2r(|)sin(f|)^ 

The expectation value //j for the flux on Earth is therefore still c/ {A'k)N but instead of the 
standard deviation, we use quantiles of the stable distribution to quantify the level of fluctuations. 
For 5(4/3, 1, 1,0, 1) the 5%, 16%, 84% and 95% quantiles are approximately -3.3, -2.6, 1.0 and 
4.7, respectively. The energy dependence of the quantiles of J is contained in 67V {c^)^^'^ oc 
^-7-3(5/4^ which is always harder than /xj oc E"^^^'^^^'^/'^ , so the fluctuations are growing with 
energy. For the particular values from Table 2.1, oc Fig. 2.2 shows the uncertainty bands 

that 68% and 90% of the total fluxes should be contained in, that is [/xj + c/(47r)6ArXig%, /xj + 
c/(47r)6jva;g4%] and [//j + c/{ATT)b]\fX^%, fij + c/(47r)67va^95%] j together with the expectation value 
for the flux /nj. 
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1 10 

Energy [GeV] 

Figure 2.2.: Fluxes of Cosmic Ray electrons from ensembles of sources uniformly distributed in a 
disk around the observer. The solid line denotes the expectation value for the sum of fluxes from 
N discrete, transient sources. The dashed lines show the expectation values if the sources are 
limited to a ring with inner radius ri (normalised to the expectation value for ri = at 1 GeV. 
The coloured bands quantify the fluctuations and contain the fluxes 68 % and 90 % of the time, 
respectively. The fluxes from 50 realisations of an ensemble of N individual sources are shown by 
the grey lines. 

One notices that a spectral softening occurs in the spectrum around 10 GeV. Below ~ 50 GeV, 
the expansion for x{^/ ^cr) made above is not valid any more and the spectrum is not dominated by 
energy losses but by losses due to the escape from the halo. We also show the the total fluxes from 
a Monte Carlo calculation of 50 realisations of a uniform source distribution with the same model 
parameters, see Table 2.1. Although the bands do indeed contain the right fraction of fluxes, the 
majority appears to be not only lower but also more curved downwards than the expectation value 
for ri = 0. 

The deficit in the majority of realisations is in fact a consequence of the stable distribution which 
is highly asymmetric with a long tail. On the one hand, it is because of the asymmetry that a 
majority of the fluxes from different realisations of the source density is below the expectation 
value. The long tail of the probability density on the other hand makes the average from a large 
number of realisations finally and slowly converge to the expectation value. This convergence is 
slowing down with increasing energy as because of the decreasing diffusion-loss length less sources 
effectively contribute at higher energies. 
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A more physical way of explaining the deficit at higher energies is that every realisation of a finite 
number of sources necessarily contains one closest source at L^i = miuj [Li] which contributes the 
most at the highest energies. Once the diffusion-loss length \{E) becomes however shorter than 
the distance L^x even this contribution gets cut-off by the exponential term in the Green's function, 
Eq. 2.15, exp (— L^j/^^). The curving of the fluxes, called propagation cut-off, is therefore nothing 
but the shoulder of (a few) Green's function(s) from the closest source(s). We note that such a 
cut-off cannot simply be modelled by reinstating the minimum radius ri in Eq. 2.36 as shown by 
the dashed lines in Fig. 2. 2 for ri = 0.1 and 0.5 kpc. 

2.3.2. A realistic source distribution 

Of course, the sources of cosmic rays are not evenly distributed in a disk around the solar system. 
Above, this was only assumed to simplify the calculation of the effect of the discreteness of sources. 
In fact, SNRs are expected to spatially correlate with star formation activity and to be mainly 
based in the thin galactic disk, tracing the spiral structure. As discussed above, we have seen that 
effects of their actual distribution will play a role at energies above ~ 100 GeV. 

Some authors [107, 271] have assumed a continuous distribution of sources for distances beyond 
a few hundred parsecs, supplemented by a set of SNRs, known from x-ray or radio surveys, for 
smaller distances. This approach is however biased by the choice of young, nearby sources which 
have been detected in radio and/or X-rays. Older sources may not be visible in photons any longer 
but still be contributing to the GCR electron flux. This relation is illustrated in Fig. 2.3. We note 
that the effect of this incomplete assumed source distribution is a dip in the electron flux seen 
in both analyses [107, 271], although at different energies because of the different diffusion model 
parameters chosen. 

Determining the complete distribution of sources in our vicinity (i.e. up to a few kpc) from 
observations seems challenging. However it turns out that we do not need to know the exact 
distribution in order to make a prediction for the e"^ flux and fraction but require only a limited 
amount of information, all of which is encoded already in the total e~ + e"*" flux. By including 
the recent measurements by Fermi-LAT [169] and HESS [186, 192] of the total e~ + flux in 
the energy region of interest, we have sufficient information at hand to make a prediction for the 
positron fraction under the assumption that the additional positrons originate in the same sources. 
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Figure 2.3.: Distance-time diagram for nearby SNRs (after Ref. [266]). Left: The open circles mark 
supernova events and the world-hnes of the discovered remnants are indicated. The thick yeUow hne 
is our past hght-cone; all events lying on it, e.g. the SNR world-lines touching it, can be observed 
presently. The blue, purple and red shadings (top to bottom) show the relative contribution of 
sources to the diffuse e~ and flux observed at Earth at 10, 100 and 1000 GeV, respectively. 
The open red circle is an example of a hypothetical supernova whose remnant is too old to be 
visible any longer but which might still be contributing to the diffuse e~ and e"*" flux. Right: A 
distance-time diagram for hypothetical nearby SNRs. The open black circles are an example of a 
possible realisation of the supernova density (the world lines have been suppressed) as simulated 
by our Monte Carlo calculation (see Section 2.3.2). 

We perform a Monte Carlo calculation by considering a large number of random distributions of 
sources drawn from a probability density function that reflects our astronomical knowledge about 
the distribution of SNRs in the Galaxy. The better the flux of e~ and from such a realisation 
of the source density reproduces the measured fluxes, the closer is the underlying distribution of 
sources likely to be to the actual one. We do not consider any scatter in the parameters of the 
SNRs but assume a prototypical set of source parameters that we determine from a compilation 
of gamma-ray SNRs, see Sec. 2.4.1. Of course all SNRs are not the same, however variations of 
the source parameters would only introduce additional fluctuations into the fluxes without altering 
their average. We can choose the best "fit" to the data and thus determine the flux. 

The smoothed radial distribution of SNRs in the Galaxy is well modelled by [272] 




(2.52) 
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Figure 2.4.: The assumed distribution of SNRs Figure 2.5.: The probabihty density for the dis- 
in the Galaxy; the cross denotes the position of tance of a SNR from the Sun. 
the Sun in between two spiral arms. 

where A = 1.96 kpc~^, tq = 17.2 kpc, 9o = 0.08 and (3 = 0.13. To obtain a realistic probability 
density for the distance between the Earth and a SNR we have to also take into account the spiral 
structure of the Galaxy. We adopt a logarithmic spiral with four arms of pitch angle 12.6° and a 
central bar of 6 kpc length inclined by 30° with respect to the direction Sun - galactic centre [273]. 
The density of SNRs is modelled by a Gaussian with 500 pc dispersion for each arm [264]. The 
resulting distribution g{r,(j)) (see Fig. 2.5) has been normalized with respect to azimuth in such 
a way that the above radial distribution Eq. 2.52 is recovered. To obtain the probability density 
for the distances we transform to the coordinates {r',(f)') centred on the Sun. As the e~ and 
fluxes are assumed to be isotropic, we can average over the polar angle 0', such that the probability 
density f^' depends only on the distance r' to the source. 



This function is shown in Fig. 2.5. 

We assume that the sources are uniformly distributed in time, i.e. their probability density ft{t) 




(2.53) 
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IS 



ftit) = { 



1/t 

max 

for < t < fmax, 

otherwise, 



(2.54) 



with tmax denoting the earhest time considered, which is again related to the minimum energy for 
which our calculation is valid through 



(2.55) 



The total number N of sources that are needed in the Monte Carlo simulation to reproduce the 
(observed) number M ~ 300 of SNRs active in the galaxy at any given time depends on the 
average lifetime of a SNR, tsnr, which is suggested to be ~ 10^ yr [274], hence 



iV = 3 X 10' 



V300y VlO^yr 



Tsnr 
lO^yr 



(2.56) 



We take a Green's function approach to calculate the contribution from these N discrete, burst- 
like sources to the total electron-positron flux in the solar system today for a halo of extent itzmax 
in z direction, again neglecting the boundaries in the radial direction. We remind ourselves that 
the Green's function for the flux of electrons from a source at r that went off a time t ago with a 
spectrum Q{E), is 



E 



n=— oo ^ u / 
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(2.57) 



where 



(2.58) 



and the diffusion length £ is defined by 
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with Zcr = 42jnax/7r. If we neglect the spatial extent of the disk and set z = 0, the function 
x{i) = x(0, i) is approximately: 



X{1) ^ { 



for / > 5, 

(2-59) 



In practice both limits can be connected at / ~ 0.66 such that the approximated x(^) has a relative 
error of at most 0.5%. We motivate the choice of the parameters of our diffusion model from an 
analysis of nuclear secondary-to-primary ratios [275]: Dq = 1028cm2s-^ 5 = 0.6, L = 3 kpc, and 
from the galactic magnetic field and interstellar radiation fields [107]: 6o = 10^^^ GeV^^ s^^. 

2.4. Fitting the Total Electron-Positron spectrum 

A schematic description of the present framework is shown in Fig. 2.6. Cosmic rays are shock 
accelerated in SNRs and then diffuse through the Galaxy to Earth undergoing collisions with 
interstellar matter en route and creating secondary . As discussed, the ratio of the secondary e"*" 
to the primary e~ from the sources should decrease with energy, in contrast to the behaviour seen 
by PAMELA. We follow Ref. [243] in explaining this by invoking a new component of which is 
produced through cosmic ray interactions in the SNRs, and then shock accelerated, thus yielding 
a harder spectrum than that of their primaries. We discuss these components in turn below and 
calculate their relative contributions by normalising to the gamma-ray flux from the SNRs, which 
provides an independent measure of the hadronic interactions therein. 

2.4.1. Primary electrons 

The radio and X-ray emission observed from SNRs is interpreted as synchrotron radiation of elec- 
trons accelerated up to energies of 0(100) TeV [274]. The spectrum of this radiation then de- 
termines the spectrum of the underlying relativistic electrons. Moreover, the theory of diffusive 
shock acceleration [276, 86] predicts similar spectra for the accelerated protons and nuclei as for 
the electrons. If the gamma-ray emission observed by HESS from a number of identified SNRs is 
assumed to be of hadronic origin, we can use the measured spectra to constrain both the relativistic 
proton and electron population. 
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Figure 2.6.: Schematic description of contributions to the galactic cosmic rays observed at Earth in 
the present framework. 



Table 2.2 shows a compilation of gamma-ray sources observed by HESS that have been identified 
as SNRs. We have included all identified shell-type SNRs and strong SNR candidates in the HESS 
source catalogue [292] (as of September 2009), and also added the SNRs IC 443, Cassiopeia A and 
Monoceros. Actually it is not clear that the acceleration of secondaries does occur in all the SNRs 
considered, especially when the gamma-ray emission is associated with a neighbouring molecular 
cloud rather than coming from the vicinity of the shock wave. In fact the gamma-rays could 
equally well be due to inverse-Compton scattering by the relativistic electrons responsible for the 
observed synchrotron radio and X-ray emission. Therefore, we have considered three possibilities 
— including all sources implies a mean power-law spectral index for the protons of (7) = 2.5, while 
excluding steep spectrum sources with 7 > 2.8 gives (7) =2.3 and excluding sources with 7 > 2.6 
yields (7) = 2.4. In the following we adopt the central value, 7 = 2.4, for the electron population 
too, unless stated otherwise. This requires a compression factor of r < 3.3 in contrast to the value 
of r = 4 expected for a strong shock, so there is clearly some tension between the DSA theory 
and observations. This can possibly be resolved if we consider only a subset of the SNRs in Table 
1 to be hadronic accelerators, or if the gamma-ray spectrum is steepened, e.g. by the onset of 
an exponential cut-off in the electron spectrum. Our model assumptions are intimately connected 
to the production of neutrinos, the detection of which will therefore provide an independent test, 
see Sec. 2.6. In this work we adopt a cut-off of Ecut ~ 20 TeV which is consistent with DSA 
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theory [274]. The source spectrum of primary electrons is then: 

E,-=i?o.(^)"'e-^/^-. (2.60) 

The normahsation is determined by fitting the electron flux at Earth resulting from our Monte 
Carlo computation to the preliminary measurement by PAMELA at 10 GeV [293]; the secondary 
fluxes can be neglected for this normalisation. We find = 1.8 x 10^*^ GeV~^ for 7 = 2.4 which 
corresponds to a total injection energy of 

^20TeV 

/ dEE R^- {E)c^l X 10'^'^ erg. (2.61) 

JlGeV 

This compares well to the value of 9.2 x 10^'' erg said to be required to power the GCR electrons [274] . 

Solar modulation which is important below ~ 10 GeV, has been accounted for using the force 
field approach [60], with a charge-independent potential ol(j) = 600 MV. However, our simple model 
ignores convection and (re) acceleration in the interstellar medium which become important below 
~ 5 GeV, hence the electron fiux cannot be predicted at lower energies. The primary e~ fluxes as 
measured on Earth for 30 different source configurations are shown in the top panel of Fig. 2.7. 
With an injection power-law index 7 ~ 2.4 it 0.1 as required for consistency with the gamma-ray 
data, there clearly is a deficit at high energies compared to the e'^ + e~ flux measured by Fermi-LAT 
and HESS. 

2.4.2. Secondary electrons and positrons from propagation 

Positrons in GCR are generally assumed to be of purely secondary origin, arising through the decay 
of pions and kaons produced in the interactions of GCR protons (and nuclei) with the interstellar 
medium (ISM) [126]. The neutral pions decay into gamma-rays which then contribute to, if not 
dominantly constitute, the galactic gamma-ray background. The charged pions on the other hand 
decay into neutrinos and muons, the latter subsequently decaying into electrons and positrons. 
Assuming that spatial and temporal variations in the GCR proton flux Jp and the ISM gas density 
?^ISM are small, the source density of these secondary background e~ and e"*" is also homogeneous, 
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Figure 2.7.: Predicted spectra of electrons and positrons with data from Fermi-LAT [169] (red 
circles) and HESS [186, 192] (blue squares &: green triangles). The diagonal arrows show the 
energy scale uncertainty. Top: Primary electrons after propagation to Earth. Middle: Secondary 
electrons and positrons from cosmic ray interactions, created during propagation (dashed line) and 
created during acceleration in SNRs (full lines^^ Bottom: The sum of primary and secondary 
electrons and positrons. 



both in space and in time: 

„ISM _ „ ^ ATP'^^ T / F^^ ^'^PP">c^+^ /'O R9\ 

Je,^, /3c d^; 

where dcrpp_>e±_|_x/di^' is the partial differential cross-section for production and /3 ~ 1 is the 
velocity of the GCR. We can then integrate the Green's function for a single source over space and 
time to calculate 

where ~ 0.1 kpc is the height of the galactic disk. 

We calculate the flux of secondary background e~ and e"*" from the Solar-demodulated flux of 
GCR protons as derived from the BESS data [294] and model the cross-sections according to 
Ref. [112]. The contribution from kaon decay is subdominant and is therefore neglected. The 
presence of He both in GCRs and in the ISM is taken into account by multiplying the proton con- 
tribution by a factor of 1.2. Our results are in good agreement with Ref. [157], taking into account 
the different diffusion model parameters and keeping in mind that convection and reacceleration 
have been neglected here. These fluxes are shown (dashed line) in the middle panel of Fig. 2.7 and 
are clearly a subdominant component which cannot account for the deficit at high energies. 

Moreover, the positron flux is falling faster than the primary e~ flux at all energies whereas the 
PAMELA data [163] clearly show a rise above a few GeV. One way this can be resolved is if there is 
a dip in the electron spectrum between ~ 10 and 100 GeV. It has been suggested that Klein-Nishina 
corrections to the Thomson cross section for inverse Compton scattering [110] or inhomogeneities in 
the distribution of sources [271] can produce such a dip. However the former would require a rather 
enhanced interstellar background light (IBL) field [110], while the latter calculation [271] assumes 
an incomplete source distribution (see Sec. 2.3.2) and moreover adopts diffusion model parameters 
quite different from those derived from the measured nuclear secondary-to-primary ratios [275] and 
the measured galactic magnetic field and IBL [107]. 

The other, perhaps more straightforward possibility is to consider an additional component of 
GCR positrons with a harder source spectrum that results in a harder propagated spectrum and 
therefore leads to an increase in the positron fraction. 
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2.4.3. Secondary electrons and positrons from the sources 

Following Refs. [243, 295], the parameters are those typical of an old SNR: ui = 0.5 x lO^cms"^, 
"-gas,! = 2cm~'^, B = IfiG. Choosing r = 3.1 to recover — 2.4 the characteristic niomentci Pcross 
and Pbreak (sec Sec. 2.2) turn out to be, 

Across = 427 K^' (^^) GeV, (2.64) 
Pbreak = 7.77^3! TeV. (2.65) 

What is still missing is the normalization of the injection spectrum, R^, in the sources which is 
proportional to the normalisation of the GCR protons, A'^gcRi through Eq. 2.3. Usually a factor 
Kcp — 10^"^ — 10~^ is introduced to normalize the electron component with respect to the protons; 
this depends on how particles are injected from the thermal background into the acceleration process 
and is not reliably calculable from first principles. We can get around this by assuming that the 
gamma-rays detected from known SNRs by HESS are of hadronic origin, as is expected in this 
framework. Thus we can use the total luminosity of individual sources in gamma-rays, 

= 4ttL'^J^ , (2.66) 

to determine the normalization of the proton component and therefore also the secondary injection 
rate q^. if we know their distance L. 

The compilation of 7 ray data on SNRs from HESS, see Table 2.2, suggests an average value 
Q'^ ~ 5.7 X 10'^^ s~^ TeV^^. We find then for the total spectrum 

Rl = tsnrQ^ ^ rsNR|±Q° , (2.67) 

where E+ (and analogously E^) is defined by Eq. 2.9, or explicitly 

Rl = 7.4 X W^^( ^) ( AGey-\ (2.68) 

+ VlOVy V5.7 X 1033s-iTeV-V 

In the Monte Carlo code we have explicitly input the experimentally measured pp cross-section 
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Table 2.3.: Summary of parameters used in the Monte Carlo simulation, for an 
injection spectral index 7 ~ 2.4. 



Diffusion Model 



Do 


1028 cm2 s-1 


1 


5 


0.6 


> from GCR nuclear secondary-to-primary ratios 


■^max 


3 kpc 


) 


u 


1 n~i6 HpV"! sj^i 

J-U \JKt V 


t rK \ rK 1 , Q T1 ri r-t C^y~i t^VCfTT' (H tin Gl "l" 1 t^G 
V^IVJ-JJ , 1 1 J 1 J dllU. -U ciici ^ V CLCllol tiCo 


Source 


Distribution 




^max 


1 X 108 yr 


from -Emin — 3.3 GeV 


TSNR 


10"^ yr 


from observations 


N 


3 X 10^ 


from number of observed SNRs 


Source 


Model 






1.8 X 10^°GeV-i 


fit to e" flux at 10 GeV 


7 


2.4 


average gamma-ray spectral index 


-E'max 


20TeV 


typical gamma-ray maximum energy 


Ecut 


20TeV 


DSA theory 


Rl 


7.4 X 10^8 GeV-i 


cf. Sec. 2.4.3 


Kb 


15 


only free parameter (for fixed 7) 



which gives a similar normalisation as the estimate presented above assuming Feynman scaling. 
Also, we have taken into account the cut-off of the underlying primary protons. Their maximum 
energy is determined from the average maximum gamma-ray energy -Emax — 20TeV (see Table 2.2) 
through the inelasticity of the pp 7 X process as ~ 20TeV/0.15 ^ 100 TeV [281]. The 
normalisation for secondary electrons is computed similarly. 

The middle panel of Fig. 2.7 shows an example of the flux of secondary source and e+ for 30 
realisations of the SNR density in our Galaxy. Clearly this component can potentially match the 
high energy Fermi-LAT and HESS data. 

We note that in our model, the contribution from secondary electrons and positrons to the total 
flux is about twice as large as in Ref. [243] where the primary injection spectrum was assumed to 
be oc E~'^, motivated by DSA theory. However this is not consistent with gamma-ray observations 
of SNRs as seen from Table 2.2. 
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Figure 2.8.: The three best fits (out of 30 source reaUsations) to the total spectrum of electrons 
and positrons measured by Fermi-LAT [169] (red circles) and HESS [186, 192] (blue squares & 
green triangles), and the corresponding prediction for the positron fraction for different values of 
7 and K-q, for both charge-sign independent (full line) and charge-sign dependent (dashed line) 
Solar modulation (see text for details). The PAMELA data [163] is shown for comparison (open 
red circles). 



2.5. Results 



The parameters used in the Monte Carlo simulation are given in Table 2.3. For an assumed 
injection spectral index 7, the only free parameter is Pcross (cf- Eq. 2.12) or, equivalently, the factor 
(cf. Eq. 2.13) which is determined by fitting the total flux of electrons and positrons to the 
Fermi-LAT and HESS data (see Fig. 2.7). 

We have calculated the with respect to the combined Fermi-LAT and HESS data for each 
realisation m of source distances and times, {Li,ti}m, over all energy bins. The three best "fits" 
are shown in Fig. 2.8 for different values of and for 7 = 2.3, 2.4 and 2.5 (see Table 2.3). The 
corresponding predictions for the fraction are shown in the bottom panels. Adopting 7 = 2.4, 
we find good agreement for Kb — 15, which corresponds to a cross-over of the accelerated and 
the non- accelerated secondary components from the sources at pcross — 28 GeV and a spectral 
break at Pbreak — 510 GeV (cf. Eqs. 2.64 and 2.65). The predictions for the fraction agrees 
reasonably well with the data down to 6 GeV; we would not expect agreement at lower energies 
since we have neglected convection and reacceleration during interstellar propagation. In fact the 
PAMELA measurements of the e"*" fraction are systematically lower than previous measurements, 
e.g. AMS-01 [148] or HEAT [165], and it has been noted that this discrepancy can be resolved by 
considering charge-sign dependent Solar modulation with (/>+ = 438 MV for and = 2 MV for 
e~ [296] (rather than (p^ = (p^ = 600 MV). This however seems to be at odds with preliminary 
PAMELA data on the absolute electron flux [293] which does show substantial Solar modulation. 
Accordingly in Fig. 2.8 we have shown the predicted e~^ fraction for both cases; note that this does 
not affect our predictions for energies above 10 GeV. 

2.6. Gamma-Rays and Neutrinos 

Our fits to both the PAMELA (absolute e~) and the Fermi-LAT (total + e^) spectra, allow to 
predict the PAMELA positron fraction by including secondary e~*" accelerated in SNRs and thus 
provides a consistent picture of current data on cosmic ray e~ and e"*" between a few GeV and tens 
of TeV. Turning the argument around, since a large fraction of the e~ and e"*" observed in GGR 
above hundreds of GeV are required to be secondaries in this model, there must be a large number 
of hadronic cosmic ray accelerators in our Galaxy, some of which should be quite nearby. 
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An independent test of the model is provided by the usual 'messengers' of such hadronic accel- 
eration environments, namely gamma-rays and neutrinos. Taking the known distribution of SNRs 
in the Galaxy (see Sec. 2.3.2) we have calculated the column depth in SNRs in the galactic disk as 
seen from Earth, 



and show this in the top panel of Fig. 2.9. As expected, the column depth is largest towards the 
galactic centre. However, the quantity that is more important for observations is the brightness 
of sources. We have therefore weighted the integrand in Eq. 2.69 by and this flux weighted 

column depth is also shown in the top panel of Fig. 2.9. We note that although the maximum 
brightness is still expected around the galactic Centre, the brightness in other directions is smaller 
by only ~ 30% because the sources in the closest spiral arms are then dominant (if they are actually 
there, of course). 

This is illustrated in the bottom panel of Fig. 2.9 by an example distribution of SNRs from 
the Monte Carlo simulation, denoted by circles. The position of the circle denotes the galactic 
longitude and the radius is proportional to the brightness in units of the Crab Nebula, i.e. an 
integrated flux of (1.98 ± 0.08) x lO^^^cm^^g-i ^bove ITeV [297]. For a source of luminosity of 
Q° = 5.7 X lO^^ TeV"^ s^^ (see Table 2.2) at distance L, the integrated flux above 1 TeV is. 



i.e. about 40 % of the Crab Nebula flux at L = 2 kpc. It is seen that although most of the sources 
are clustered towards the galactic centre, there are several bright sources at large longitudes as 
well. We find typically ~ 3 sources brighter than the Crab (or ~ 7 brighter than 50 % Crab). 

The adopted distribution of SNRs (Sec. 2.3.2) and the average luminosity per source determined 
from a compilation of known sources (Table 2.2) thus leads to the prediction of several nearby 
SNRs with fluxes of the order of the Crab Nebula. Note, however, that close sources could be 
rather extended and thus have escaped detection by HESS in one of its surveys of the Milky 
Way [286, 298, 299]. For example, a diameter of ~ 50 pc which is a typical value for a very old 
SNR, corresponds to 1.5° at 2 kpc. 




(2.69) 




(2.70) 
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Figure 2.9.: Top: The column depth and flux weighted column depth of the SNR density in the 
galactic plane. Bottom: Example of a distribution of SNRs in gamma-rays/neutrinos from the 
Monte Carlo simulation. The position of a circle denotes the galactic longitude of the source and 
the radius is proportional to the brightness in units of the Crab nebula. One source whose circle 
exceeds the vertical scale is ~ 500 pc from Earth and has a total integrated flux above 1 TeV of 
~ 6 times the Crab Nebula. 

Extended gamma-ray luminous SNRs can however be detected by MILAGRO [300] with its larger 
field of view. A survey in galactic longitude I £ [30°, 220°] and latitude b G [—10°, 10°] has revealed 
6 new sources at a median energy of 20 TeV, several of which are spatially extended. The flux 
from a SNR of the above luminosity at L = 2kpc is QyiiircP) ~ 1.2 x 10"^^ TeV"-^ cm"^ s"^ 
at 1 TeV. Scaled with a spectral index of 2.4 to 20 TeV, this gives Q°/(47rd2) 20-2-4 ~ 9.0 x 
-|^Q-i5 XeV^^ cm^^ s^^ which is in the range of the unidentified MILAGRO sources [300]. We note 
that the MILAGRO source MGRO J1908+06 was recently confirmed by HESS [301], though with 
a smaller angular extent of ~ 0.7°. However, correlating unidentified MILAGRO sources with the 
FERMI Bright Source List [302, 303] seems to favour associations with pulsars, although several 
new unidentified extended sources have also been found. 
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Hadronic sources of cosmic rays should also be visible by their neutrino emission. On general 
grounds, the neutrino luminosity (from vr^ decay) can be directly related to the gamma-ray lumi- 
nosity (from TT^ decay) and should be of the same order of magnitude since p-p interactions produce 
vr"'", tt" and tt~ in roughly equal numbers. Each of the three neutrinos produced in the decay chains 
vr"^ —7- fi^i'^ —7- e^VeV^v^ and vr^ — t- fi^v^ — t- e~i>eV^v^ carries about half of the energy of each 
photon produced in the decay vr*^ — )• 77. Hence, the ratio of neutrinos to photons produced on 
average is ~ 3 : 1 and the total neutrino luminosity is 




Presently the largest cosmic neutrino detector is the IceCube observatory [246] under construction 
at the South Pole. IceCube observes high energy neutrinos via their interactions with nucleons in the 
vicinity of the detector and subsequent Cerenkov light emission of energetic charged particles in the 
transparent glacial ice. The most important signal for neutrino astronomy is the Cerenkov radiation 
by muons produced via charged current interactions of muon neutrinos. Since the muon inherits the 
large boost of the initial neutrino the point source resolution is ~ 1°. The large background signal 
of atmospheric muons is efficiently reduced for upward-going muons, i.e. neutrino sources which are 
somewhat below the horizon. Hence, IceCube is mainly sensitive to neutrino point sources in the 
northern sky, which excludes SNRs in the direction of the galactic centre. 

Neutrino emission associated with galactic TeV gamma-ray sources has been investigated by 
many authors [304, 305, 306, 307, 308, 309, 310, 311, 312] including also the HESS sources used in 
our analysis. In particular, Ref. [308] investigates the prospects of neutrino detection for the SNRs 
HESS J0852.0-463, J1713-381, J1804-216, J1834-087 (see Table 2.2) in the proposed KM3NeT 
detector in the Mediterranean which will see the galactic centre region. The muon neutrino rate is 
expected to be a few events per year for such sources. 

Due to flavour oscillations of neutrinos with large mixing angles, the initial flavour composition 
Qve '■ Qu^ ■ Qut ^1:2:0 from pion decay is expected to become ~ 1 : 1 : 1 at Earth. The TeV 
muon neutrino point flux from a hadronic gamma-ray source located at a distance L and with a 
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power-law index 7 ~ 2.4 is thus F^^{> 1 TeV) ~ 2^-^'F^{> 1 TeV), hence 

F^(> 1 TeV) ~ 3.2 x lO^^^ f -J^\ 011'^ s~^. (2.71) 

\ 2 kpc / 

This should be compared to the results of searches for neutrino point sources in the northern sky, 
in particular the close-by SNR Cassiopeia A (see Table 2.2), using data taken with AMANDA-II 
(the predecessor of IceCube) during 2000-2006 [313] and, more recently, with the first 22 strings of 
IceCube during 2007-08 [314]. The average 90% C.L. upper limit on the integrated flux in the 
energy range 3 TeV to 3 PeV is [314] 

Fu^ < 4.7 X 10-^2 cm-^g-i, (2.72) 

i.e. still above ~ 7 x 10^^'^ cm^^g-i expected from a SNR at 2 kpc, assuming 7 = 2.4. 

The full 80 string configuration of IceCube thus has excellent prospects to identify these SNRs. 
A point source in the northern sky with an E~''^ muon neutrino flux, 

F^^ ~ 7.2 X IQ-^^cm'^ s-\ (2.73) 

in the TeV-PeV range can be detected with a 5a significance after three years of observation. This 
does depend somewhat on the spectral index and energy cut-off, since the signal (after "level 2 
cuts") peaks at an energy of ~ 10 TeV [246]. As mentioned previously, our analysis predicts on 
average ~ 3 nearby gamma-ray sources stronger than Crab with corresponding muon neutrino 
fluxes larger than ~ 7 x 10^^^ cm^^ g^^. Note that although the galactic centre is not in the field of 
view of IceCube, SNRs following the spiral arm structure of the Galaxy are expected to be detected 
also in the galactic anti-centre direction, as seen in the example distribution shown in the bottom 
panel of Fig. 2.9. 

2.7. Comments on the Time-Dependent Picture 

It was recently questioned [315] whether the acceleration of secondaries also holds true in the 
time-dependent picture of DSA. Intuitively, one would expect that the energy-dependence of the 
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diffusion coefficient still leads to an increase of the size of the diffusion zone and therefore to a 
harder spectrum for charged secondaries. The authors of [315] however argue that both numerical 
simulations and a simple analytical argument show that the secondary fractions do not rise. In the 
following, we critically examine both these claims. 

2.7.1. An analytical argument 

The authors of [315] present a simple 'gedankenexperiment' which they claim shows that the ac- 
celeration of secondary mechanism does not lead to a rise in any secondary-to-primary ratio. The 
argument compares the flux of secondaries produced and accelerated in the SNR to the flux of 
primaries in the case where no secondaries get produced (which would be the case, for example 
for 

TXgas — 0)' For the particular case of the anti-proton fraction, assuming that in an inelastic pp 
collision, the anti-proton produced takes all the energy of the incoming proton, the interactions 
would in fact only convert part of the proton flux into an anti-proton flux. As both particles are 
affected in the same way by DSA, the total anti-proton spectrum would be proportional to the 
proton spectrum in a secondary-less case. 

This comparison is however not appropriate. What needs to be compared with the proton 
spectrum /* in a secondary-less case is not the antiproton flux fp but the sum of the proton and 
antiproton flux {fp + fp) for the case with secondaries. If the assumptions about the inelasticity 
being 100 % were in fact correct, the sum of antiproton and proton flux {fp + fp) should indeed be 
comparable to the proton flux /*. This, however, does not tell us anything about the individual 
spectra, fp and fp. In fact, the hardness in the antiproton spectrum, fp oc pf* would be balanced 
by a softer proton spectrum fp = {1 — k {p/pcT))fp where pcr is the momentum at which fp = fp and 
K is a normalisation constant. Therefore, the antiproton-to-proton ratio, fp/ fp would become even 
harder and in this sense our above argument is conservative in that it ignores the —pfp reduction in 
the proton flux. However, for the case of antiprotons, the latter correction is negligible as fp » fp, 

i.e. Per Pmax' 

2.7.2. The numerical model 

The authors of [315] further present the results of a Monte Carlo calculation of time-dependent 
DSA of electrons and protons including the production and subsequent acceleration of secondary 



97 




Figure 2.10.: Antiproton-to-proton ratio pro- 
duced in a time-dependent Monte Carlo calcu- 
lation of diffusive shock acceleration in the test 
particle approximation (from [315]). The red 
(A) and blue (B) lines denote the ratio from the 
secondaries which have been produced upstream 
and downstream, respectively, and the black line 
(A+B) is the sum of both. 



Figure 2.11.: Shock velocity of the SNR as a 
function of time, both in units of values charac- 
teristic for a particular SNR model, see [79]. The 
blue solid and green solid lines show the analyt- 
ical model of [79] for the ejecta-driven (ed) and 
the Sedov- Taylor (ST) phase, respectively. The 
effective value for the shock velocity adopted in 
our calculation is 500kms~^. 



antiprotons and positrons. Although the secondary-to-primary ratios of the total spectra N{p) oc 
J dVp'^f{p) clearly show the presence of a harder component compared to the primary one, the 
former one seems to be subdominant and not influence the total spectra too much. This is shown 
in their Fig. 2 which we here reproduce in Fig. 2.10. 

However, we believe that this discrepancy is merely due to a different set of parameters adopted. 
As shown in Eq. 2.8, the normalisation of the harder component depends on the diffusion coefficient 
D, the velocity of the shock front squared, uf and the gas density 

Tigas through Q-^ . Not only is the 

value of the diffusion coefficient in [315], -DBohmj smaller than the one adopted in our calculation, 
KsD-Qohm with Kb ~ 20, but also is the velocity larger: The dynamical SNR model adopted [79] 
predicts a time-dependent velocity, both for the free expansion and in the Sedov- Taylor phase of 
a SNR. Our time-independent analysis, however, needs to adopt an effective velocity. As we show 
in Fig. 2.11, the time-dependent velocity is higher than our effective value most of the time. As 
the velocity enters as the inverse square into the normalisation of the harder component, a small 
difference in velocity can change the relative importance rather drastically. 

We therefore conclude that the discrepancy between our analytical result and the outcome of 
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the Monte Carlo calculation performed in [315] is most likely due to different parameters for the 
adopted SNR model. We note that it has not been shown that the time-dependence introduces any 
additional physical effects to the time-independent picture discussed above. 

The true value of the parameters is, however, a separate issue. One possibility, in fact the one 
pursued in [315], is to adopt an (analytic) model of the SNR dynamics that is believed to reflect the 
physical behaviour of actual SNRs and feed it into a numerical calculation of the secondary spectra. 
However, even the parameters of the analytical model are affected by uncertainties. Therefore, one 
needs to consider different SNR models to cover all possibilities of effects. 

Another possibility is to adopt the simplest analytical picture possible, like the simple time- 
independent test particle approximation above, and to treat the parameters, e.g. diffusion coef- 
ficient, velocities, maximum energy, as effective ones. As explained above, these effective values 
only enter the normalisation in a certain combination. Choosing a value that allows to reproduce 
measured data, i.e. basically fitting these parameters to the data, allows to make independent, 
testable predictions for other observables. We have followed this approach already for the charged 
lepton channels in Sec. 2.4, and we will further use this approach for nuclear secondary-to-primary 
ratios in Chapter 3. 

2.8. Alternative Astrophysical Explanations 

We have already addressed the possibility of DM annihilation or decay accounting for the excesses 
in the positron fraction and total electron-positron flux. Furthermore, we have discussed some of 
the constraints, for example, from antiprotons, radio waves and gamma-rays, and we have argued 
that a large fraction of the parameter space necessary to explain the leptonic anomalies are already 
ruled out. 

To round off our discussion of the acceleration of secondaries in old SNRs, we present some 
alternative astrophysical explanations that have been claimed to explain the excesses. 

2.8.1. Pulsars 

The idea that pulsars might be responsible for the lepton excess [123, 124] has been around since 
the first indications of a rise in the positron fraction [316]. Pulsars are highly- magnetised and 
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quickly spinning neutron stars that transfer (part of) their rotational energy into a magnetised 
wind. Electrons are accelerated to very high energies by the electric fields present and synchrotron- 
radiate gamma-rays in the extreme ambient magnetic fields. These gamma-rays are energetic 
enough to pair-produce high energy electrons and positrons on the ambient magnetic fields or the 
thermal X-rays form the pulsar itself. The high energy gamma-ray emission is usually interpreted 
as curvature radiation or ICS from high-energy electrons and positrons. There is disagreement 
about where exactly the gamma-ray emission is produced: close to the surface of the neutron star 
(polar cap models) or further away in the light cylinder along the last open field lines (outer /slot 
gap models). 

The total output in electrons and positrons is usually connected to the spin-down power of the 
pulsar and with an efficiency of a few per cent, most of the pulsars in the ATNF catalogue [317] 
are predicted to provide outputs of lO^^---^'^ erg. In fact, probably only mature pulsars, i.e. with 
s-gss > 10^ yr, can inject electrons and positrons into the ISM as for younger ones all charged 
particles are trapped inside the pulsar wind nebula (PWN). The spectral index a of the electrons 
and positrons can be estimated from radio emission interpreted as synchrotron (7 = 1 . . . 1.6) or 
from the gamma-rays (7 = 1.4 . . . 2.2). The range usually considered is 7 = 1.4 . . . 2.2. 

As already mentioned, the idea was first developed in Ref. [123], stressing the importance of 
young and nearby pulsars for the electron-positron flux and the positron fraction at high energies 
and suggesting Geminga as one of these possible sources. The positron excess can be reproduced 
for a few mature, nearby pulsars, e.g. Geminga and Monogem, complemented by a distribution of 
pulsars further away [318]. This is in supported by [170] where it is concluded that the scenario 
with a single pulsar is disfavoured. Other authors focus on the one-source model [319] although 
this might require too high an efficiency for the conversion of spin-down power into electrons and 
positrons. It was also shown that the sum of all pulsars covering a broad range of parameters can 
reproduce the Fermi-LAT and PAMELA data [185]. 

The way to pinpoint to a (few) pulsar(s) as the origin of the lepton excesses suggested is to look 
for anisotropics in the arrival directions of electrons and positrons, e.g. with Fermi-LAT. A simple 
estimate of the anisotropy shows that most favoured scenarios could be detected within a few years 
time to 2 or 3 fj significance [318, 170]. Nearby pulsars as the source of the positrons are of course 
quite consistent with the absence of antiprotons. 
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2.8.2. GRBs 

Any astrophysical positron source with an injection distinctively harder than the primary electrons 
from SNRs can in principle explain the PAMELA positron excess. It was recently suggested [320] 
that a galactic gamma-ray burst (GRB) could be such a source. The GeV electrons and positrons 
would be pair-produced by scattering of TeV gamma-rays on the afterglow eV background. The 
energy contained in TeV gamma-rays is about 5 % of the total GRB energy, typically ~ lO^'' erg 
and the spectral index of the electron-positron pairs is 1.8 it 0.7 [320]. 

In fact, an injection of 0{W^^) erg a few times 10^ yr ago and with a spectral index between 1.6 
and 2.2 nicely fits the data. A harder spectral index can also reproduce a ATIC/PPB-BETS-like 
feature in the total electron-positron flux while a softer index is in agreement with the Fermi-LAT 
data. If in fact a GRB or another singular event of similar characteristics (pulsar, hard SNR, 
microquasar) was responsible for the lepton excesses one would again expect to see a rather large 
anisotropy in electron and positron fluxes which would be detectable by Fermi-LAT or AMS-02 [321] 
in a few years' time [320]. 

2.8.3. Very old supernova remnants 

Another idea [322] invokes dense gas clouds around very old SNRs as sites of a production of harder 
positrons. The authors consider several simultaneous SN explosions 10^' '^yr ago inside a dense gas 
cloud that need to have taken place only hundreds of parsecs from the solar system but could explain 
the local bubble or other similar structures like Loop I [323, 324]. Once the Sedov- Taylor phase is 
coming to an end, the compression ratio r (see Sec. 1.2.3) is not given by the Rankine-Hugoniot 
relation any more but by the ratio of shock velocity and upstream Alfven velocity, r ~ ^/2vs/va, 
and is much larger than 4. The spectral index 7 of the protons accelerated by DSA is therefore 
smaller than 2 and the interaction of the protons with the gas in the surrounding cloud produces 
positrons of a similarly hard spectrum. 

For gas clouds of tens of parsecs radius and ~ 100 cm~^ density, both the positron fraction and, 
depending on the exact value of the proton spectral index 7 considered, even the Fermi-LAT or 
ATIC/PPB-BETS total electron-positron spectrum can be fitted. However, the maximum proton 
energy assumed is 100 TeV - a somewhat unrealistic situation, given that such high energy protons 
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can only be accelerated in the Sedov- Taylor phase when the magnetic field amplification is strong 
and even the highest energy particles are confined to the shock region. 

The model predicts a similar rise in the antiproton-to-proton ratio, due to the contributions from 
antiprotons produced in a similar manner. Testing this model by the boron-to-carbon ratio might 
not be conclusive, as the metallicity can vary throughout the Galaxy. 

2.8.4. Inhomogeneous source distribution 

We have already alluded to the suggestion that the inhomogeneous distribution of SNRs in the 
galactic disk is responsible for the rise in the positron fraction [271]. The idea is that the SNR 
density in the spiral arms is amplified by a factor of 4 with respect to the disk, and the propagation 
cut-off from the nearest spiral arm at ~ 1 kpc leads to a softening of the primary electron spectrum. 
The secondary production is spatially more homogeneous because it follows the cosmic ray proton 
density and is therefore not affected by the details of the source distribution. The softening of the 
electron spectrum therefore leads to a rise in the positron fraction. 

While we agree with the call for a realistic source distribution, we argued above in detail that its 
effect is only to be seen at hundreds of GeV. In particular, we expect the propagation cut-off from 
sources at distances of ~ 1 kpc to occur at 0(1) TeV, i.e. more than a magnitude above the energies 
at which the positron fraction starts rising. This disagreement can be traced back to propagation 
parameters used in [271] that are rather unrealistic. In particular, the diffusion coefficient adopted, 
D{E) = 6 X 10^'' cm^ s~^, is too low and the energy loss rate at 1 GeV, bo = 1.8 x 10^^^ GeV~^ s~^ 
too high, which strongly affects the diffusion-loss length see Eq. 1.48. In addition, the predicted 
soft electron spectrum seems to be in disagreement with preliminary PAMELA data on the absolute 
electron flux [164]. 

2.9. Conclusion 

Supernova remnants have long been suspected to be the sources of galactic cosmic rays. We have 
discussed a recent proposal [243] that proton-proton interactions in the shocks of SNRs followed 
by the diffusive shock acceleration of the secondary positrons produced can flatten the spectrum 
of the secondaries relative to that of the primaries. These hard spectra may be the origin of the 
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recently observed cosmic ray "excesses" — both the e+ fraction observed by PAMELA [163] and 
the e" + e+ flux measured by Fermi-LAT [169] and HESS [186, 192]. 

We have investigated how gamma-ray emission of SNRs - assumed to be of the same hadronic 
origin as the positrons ~ together with cosmic ray data, constrain the acceleration of positrons. 
We have accounted for the spatial and temporal discreteness of SNRs via a Monte Carlo exercise, 
drawing samples from a realistic galactic distribution with the observed SN rate. For the diffusion 
parameters we have adopted standard values derived from cosmic ray nuclear-to-primary ratios, as 
well as the energy densities of galactic radiation and magnetic flelds. 

We have compiled a list of all gamma-ray emitting SNRs observed by HESS and determined the 
mean value of the flux, which fixes the hadronic interaction rate in the SNR. Low energy data from 
PAMELA on the absolute e~ fiux were used to normalize the primary flux of . The contribution 
from accelerated e"*" was then found by fitting the e~ -|- fiux to Fermi-LAT and HESS data, 
adjusting the (only) free parameter which determines the diffusion rate near SNR shocks. 

The spectra of e"*" and e~ thus derived agree well with the e"*" fraction observed by PAMELA 
in the range 5 — 100 GeV. The apparent discrepancy at lower energies can be attributed to the 
uncertainty in solar modulation (charge-sign dependent or independent). Furthermore, convection 
and diffusive reacceleration of primary electrons that become important at these energies were 
neglected in our analysis. The flux of and e~ becomes dominated by the accelerated secondary 
component at high energies; the corresponding fraction levels out at ~ 0.4, reflecting the relative 
multiplicity of and e~ produced by p-p interactions. 

To be consistent with our overall framework the gamma-rays observed from SNRs have been 
assumed to be of hadronic origin. The known spatial distribution of SNRs then implies (on average) 
several nearby sources with a gamma-ray flux comparable to the Crab. We have speculated that 
some unidentifled MILAGRO sources [300] might correspond to such old SNRs. Moreover, the same 
hadronic processes in SNRs will inevitably produce high energy neutrinos which can be detected 
in cubic-km telescopes such as IceCube [246]. The neutrino luminosity can be directly related to 
the gamma-rays and is not connected to the hypothetical acceleration of e"*" and e~ in the sources 
as in our present model. Nevertheless, similarly to the previous argument, we expect on average a 
few nearby sources, some of which may also lie within the fleld of view of IceCube and can thus be 
detected with high statistical signiflcance after three years of data taking. 
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While our calculational framework is based on first-order Fermi acceleration by SNR shock waves, 
we have noted that in detail the observations do not fit the theoretical expectations, e.g. the shock 
compression ratio inferred from the observed gamma-ray spectrum (~ is 3.1 rather than 4 

as is expected for a strong shock [276]. Going beyond the test particle approximation, the generic 
expectation in such a process is for particle spectra which are much flatter than those observed 
(~ E~^'^ and slightly concave), when the back reaction of the cosmic rays on the shock is taken 
into account [86]. By contrast, the observed radio spectrum of Cassiopeia A is slightly convex and 
this, as well as the morphology and time evolution of radio emission from such young SNRs, can be 
well explained in terms of second-order Fermi acceleration by plasma turbulence behind the shock 
wave [325]. Moreover the observed spatial correlation between the gamma-ray emission and the 
hard X-ray emission from some SNRs argues for a leptonic rather than hadronic origin and further 
observations are necessary to resolve this issue [326]. It has been argued that cosmic ray protons 
and nuclei may well have different sources (e.g. "superbubbles" formed by multiple supernovae) 
than the cosmic ray electrons [327]. The additional predictions made in this work concerning the 
visibility of hadronic accelerators in gamma-rays and neutrinos, tied to the expectations for the 
fluxes of the accelerated secondary positrons in cosmic rays, will hopefully enable further consistency 
tests of the SNR origin hypothesis for galactic cosmic rays. 
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3. Acceleration of Secondary Nuclei 



3.1. Introduction 

As we discussed in Chapter 2, there is a wealth of suggestions on how to explain the apparent 
excesses in the positron fraction and in the total electron-positron flux as measured by PAMELA 
and Fermi-LAT, respectively. In particular, we have devoted our analysis to the investigation 
of the acceleration of secondaries mechanism which can explain the excesses by considering the 
production and subsequent acceleration of secondary electrons and positrons in the cosmic ray 
sources, i.e. supernova remnants (SNRs). This process is guaranteed in the sense that secondary 
particles will necessarily be generated by the spallation of primary cosmic rays on ambient matter 
in the sources. As discussed above, the crucial point is the normalisation of this harder component 
which depends on many parameters that are not directly accessible to observations, for example, the 
ambient gas density. One parameter in particular, the diffusion coefficient in the cosmic ray source, 
plays an important role in determining the relative size of the effect of acceleration of secondaries. 

One possible way to resolve this problem is to determine the normalisation from other observables. 
If it is in fact acceleration of secondaries that leads to a rise in the positron fraction then both 
the ambient density and the diffusion coefficient are large enough such that other secondaries are 
produced in equally important abundances and also get accelerated. Therefore, we should see 
similar features of a rise in other secondary-to-primary ratios. We will use observations of such 
observables as the antiproton-to-proton ratio, titanium-to-iron (Ti/Fe) and boron-to-carbon (B/C) 
to determine the normalisation of the harder electron-positron component and thereby test the 
model of the acceleration of secondaries. 

In particular, nuclear secondary-to-primary ratios like B/C not only allow us to test this particular 
model but also to discriminate it against other explanations of the rise in the positron fraction 
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Table 3.1.: Truth table representing the logical possibilities 
to explain rises in different secondary-to-primary ratios. 
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/e^ , e.g. dark matter annihilation or decay (see Sec. 1.3.1) and pulsars (see Sec. 2.8). For 
instance, if a feature similar to the rise in the positron fraction was also observed in the antiproton- 
to-proton ratio p/p, this would rule out^ the pulsar explanation of the positron excess as pulsars 
are not expected to produce any antiprotons. Conversely, and as such a rise is not observed up 
to the maximum energies of current observations, ~ 130 GeV [149], this sets some constraints on 
models which in principle predict such a rise in p/p. There have, however, been efforts to build 
dark matter (DM) models that are leptophilic, that is DM predominantly annihilating or decaying 
into leptons, see Sec. 1.3.1. Similarly, a rise in nuclear secondary-to-primary ratios, like B/C would 
rule out the DM explanation of the positron excess as DM annihilation/decay should not produce 
nuclei. Again, the non-observation of such a rise up to energies higher than currently accessible, 
would in turn rule out the acceleration of secondaries model^. These relations are summarised in 
Table 3.1. 

3.2. Antiproton-to-Proton Ratio 

The calculation for the antiproton-to-proton ratio [295] is conceptually very similar to that for the 
positron fraction. Of course, both the production cross section and the inelasticity are different for 
antiprotons than for positrons. Furthermore, as both protons and antiprotons are stable and their 
energy losses are negligible above a few GeV, their fluxes can be calculated in the simple leaky box 



There is a caveat: If in fact neither a dark matter model nor the acceleration of secondaries could explain the full 
positron excess when fitted to the potential p/p feature there might still be room for a contribution from pulsars. 
^The acceleration of secondaries model could be saved by arguing that GCR protons and heavier nuclei are produced 
by different source populations and that the acceleration of secondaries only plays a role in the sources of protons. 
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model with a grammage parametrised similarly as in Sec. 1.2.7, 



E V 
10( 



^esc(^) = X, ( J , (3.1) 



with Xq = 6gcm~^ and 5 = 0.6. Above lOGeV, solar modulation is negligible, too. 

The antiproton-to-proton ratio now has two contributions: the usual one from the antiprotons 
produced in spallation of GCR nuclei on the ISM, 

Jp,ism{E) 2eXcsc{E) 



Jp{E) m. 

and one from the antiprotons produced inside the SNRs, 

Jp,SNRs{E) 



dEoEo'-^-'a,p{Eo, E) , (3.2) 



JpiE) 

where 



c^2nic[A{E)+B{E)], (3.3) 



j[^E) = a(l+A X [ dwL^'^-^^^^ / dEoEo^-''app{Eo,oj) (3.4) 
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is from antiprotons produced in the diffusion zone and 

B{E) = dEo Eo'-^ appiEo, E) . (3.5) 

is from the antiprotons produced downstream. Here, r = 3.8 is the compression ratio of the shock 
chosen such that the proton's source spectral index in phase space a = 3r/(r — 1) ~ 4.1. The 
inelasticity of the antiproton production ^ ~ 0.17, the SNR lifetime tsnr = 10^ yr, the upstream 
velocity ui = 0.5 x lO^cms^^ and the fudge factor Kb = 20 in the diffusion coefficient D = 
3.3 X 10^^i^'^(£'/GeV)'^ cm^ s~^. The maximum energy is fixed to -Bmax = lOTeV and for the 
production cross section app a parametrisation [112] is used. The factor e ~ 1.26 takes into account 
the contribution from CR nuclei and nuclei in the ISM on the production of antiprotons. We 
reproduce the antiproton-to-proton ratio calculated in Ref. [295] in Fig. 3.1. 
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Figure 3.1.: The antiproton-to-proton ratio (adapted from [295]). The A term (dotted line) and B 
term (dot-dashed Hne) are the secondaries produced in the diffusion zone and in the downstream 
region, respectively. The data points are from PAMELA measurements [328, 149]. The antiprotons 
produced by spallation of GCR nuclei on the ISM are shown by the dashed line. The sums of B term 
and ISM contribution do not show a qualitatively different behaviour than the ISM contribution 
alone. The A term finally leads to the upturn in the total antiproton-to-proton ratio. 

3.3. Nuclear Secondary-to-Primary Ratios 

3.3.1. Timescales of the problem 

The transport equation for any nuclear species i reads 



= Di—^ + Tifi + Qi, 

ox ox'^ 3 ax op 



(3.6) 



where fi is the phase space density and the different terms from left to right describe convection, 
spatial diffusion, adiabatic energy losses as well as losses and injection of particles from spallation 
or decay. We consider the acceleration of all species in the usual setup: in the frame of the 
shock front the plasma upstream {x < 0) and downstream (x > 0) is moving with velocity n_ 
and n+ respectively. We solve Eq. 3.6 analytically for relativistic energies Sk greater than a few 
GeV/nucleon such that p « i?, /3 ~ 1 and NidE ~ Airp'^fidp. At these energies ionization losses 
can be neglected and the spallation cross sections become energy independent. 
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There are three relevant timescales in the problem 
1. Acceleration time race (cf -[86]): 

3 rfD+ D;\dp' 



^'^'^ U- — u+ Jo V ^+ J P 

for Bohm diffusion and the parameter values mentioned later. 

2. Spallation and decay time Tj. 

^spall ^ -^/^spall ^ 2 X 10^ (^)~' (3-8) 

where an average ai of 0(100) mb has been assumed. The rest lifetime rf'^'^ of the isotopes 
considered ranges between 4 x 10"^ yr and 10^"^ yr. 

3. Age of the SNR under consideration (see Sec. 2.2) 

rsNR = a;max/^i+ ~ 2 X lO'^yr. (3.9) 

There are two essential requirements for a SNR to efficiently accelerate nuclei by the DSA mech- 
anism: 



(a) Tacc ^ T^^^"^"^, which is equivalent to 



20^^ < 1 ^ ek < 6.4 X lO^^^^cGeV. (3.10) 



u 



(b) TsNR ^ Ti which implies, 

— X < 1 . (3-11) 

n+ 1 i M+ 

The isotopes for which condition (b) is not satisfied at the lowest energy considered viz. ^^Ni, 
^''Co, ^^Fe, ^^Mn, ^^Cr, ^^V, ''^Ti and ''Be do not contribute significantly, so their decays in the 
source region are neglected. 
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3.3.2. Nuclear spectra at source 

We find that the general solution to Eq. 3.6 for x 7^ is 



ft = i: (^J^e^"^/^ + ^te^"^/^) + Gf, (3.12) 



With xf = ^{i-^Ji + ^Dfrf/ul) 



i 

where Gf is the asymptotic value and E^^ and -F^^ are determined by the recursive relations: 



DfAf -2nAf -4r± ' 



-4V F^*" r^*" 
DfKf-2uKf-4Tf 



We require that the phase space distribution function converges to the adopted primary composition 
Yi (at the injection energy pq) far upstream of the SNR shock: 

fi{x,p) ^z^z^ Yi,6{p - po) , dfi/dx (x,p) ^llz^ 0. (3.15) 

We also require the solution to remain finite far downstream. As the phase space density is con- 
tinuous at the shock front, we connect the solutions in both half planes to f-' = fi{x = 0,p) and 
find them to be: 



j<i 

f+ = jyt-/^ + Y.^^^ P^^"^' - ^''^'■^') + Gt (1 - e^^"/2) . (3.17) 

Using Eqs. 3.10-3.11, we can linearly expand and k,~ in Eq. 3.13 and the exponentials in 
Eqs. 3.16-3.17, 

^ U- ' 
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to obtain: 



X . 



(3.19) 



where qf^ denotes the downstream/upstream source term: qf = Ylj<i fj^f-^i- 



Finally we integrate the transport equation over an infinitesimal interval around the shock, 
assuming that q^ /q^ = Tf /T^ = n+ /ng^s = and that ~ D^: 



P^ = -af?-a{l + r')^^f^ + a 
op u_ 



(1 + ^2)C(^^0)D^^ _ 

u_ 



(3.20) 



where a = 3r/(r — 1). This is readily solved by 



^ dp' f p'^^ 



P \P 



-ail+r^){D-{p)-D-{p'))r-/ul 



X a 



[1 + r 



^ ^r(x = 0)Dr{p') 



+ Y,5{p'-po) 



(3.21) 



The Eqs. 3.19-3.21 should be compared to Eqs. 2.5-2.8 in Sec. 2.2 where the loss terms Fj/j for ion- 
isation was not taken into account but cooling losses of the electrons and positron were introduced 
in an ad hoc fashion by an exponential cut-off at E'cut — 20 TeV. This approach was justified for 
electrons and positrons as cooling losses have been shown to produce exactly such a functional be- 
haviour with Ecut in the right range [329]. The exponential in our Eq. 3.21 leads to a natural cut-off 
in both the primary and secondary spectra above the energy predicted by Eq. 3.10. However, due 
to the approximations we have made, the secondary-to-primary ratios cannot be predicted reliably 
for AFiDi/u^ ^0.1, i.e. much beyond ~ 1 TeV. 

Starting from the heaviest isotope, Eqs. 3.19 and 3.21 can be solved iteratively to obtain the 
injection spectrum after integrating over the SNR volume. 



Ni{E) = An 



dxp^ flip) Air dp' 



(3.22) 



3.3.3. Propagation of nuclei 

To account for the subsequent propagation of the nuclei through the ISM we solve the trans- 
port equation in the leaky box model, see Sec. 1.2.6, which reproduces the observed decrease of 
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secondary-to-primary ratios with energy in the range ~ 1 — 100 GeV by assuming an energy- 
dependent residence time. The steady state cosmic ray densities A/i observed at Earth are then 
given by recursion, starting from the heaviest isotope i = 1, 

= ^ — —( , (3.23) 

J-/ ' esc,j ~r J- j 

where TZsN ~ 0.03 yr^^ is the Galactic supernova rate. 



3.4. Parameters 

We calculate the source densities Ni and ambient densities A/i, taking into account all stable and 
metastable isotopes from ^^Ni down to ^^Cr/'^^Ca for the Ti/Fe ratio, and from ^^O down to ^'^Be 
for the B/C ratio. Short lived isotopes that decay immediately into (meta)stable elements are 
accounted for in the cross-sections. The primary source abundances are taken from Ref. [142] and 
we have adopted an injection energy of 1 GeV independent of the species. The partial spallation 
cross-sections are from semi-analytical tabulations and the total inelastic cross-sections is obtained 
from an empirical formula [330, 331, 332]. The escape time is modelled according to the usual 
relation: 

resc,^ = (pc)-' Xesc,^ = (p c)-' X^.^E/Z,)-' (3.24) 

where Xesc,i is the column density traversed in the ISM and p = 0.02 atomcm"^ is the typical mass 
density of hydrogen in the ISM. We have neglected spallation on helium at this level of precision 
as its inclusion will have an effect < 10%. The fit parameters are sensitive to the adopted partial 
spallation cross-sections, for example 6 ~ 0.7 for the Ti/Fe ratio but ~ 0.6 for the B/C ratio. 

The parameters are chosen as in Sec. 2.4.3: r = 4, u_ = 0.5 x lO^cms^^, = 2cm^^ and 
B = 1 /zG. The diffusion coefficient in the SNR is 

Di{E) = 3.3 X lO^^i^B B^^EQ^yZr^ cm^s"^ (3.25) 

where the fudge factor Kb is the ratio of the diffusion coefficient to the Bohm value and is deter- 
mined by fitting to the measured Ti/Fe ratio. 
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3.5. Results 



The calculated Ti/Fe ratio together with relevant experimental data is shown in Fig. 3.2. The 
dashed line corresponds to the leaky box model with production of secondaries during propagation 
only and is a good fit to the (reanalysed) HEAO-3-C3 data [334]. The solid line includes production 
and acceleration of secondaries inside the source regions which results in an increasing ratio for 
energies above ~ 50GeV/n and reproduces well the ATIC-2 data [333] taking Kb — 40. This 
is somewhat higher than the value of Kb = 10 ... 20 determined from a fit of the total electron- 
positron flux to the Fermi-LAT and HESS data (see Sec. 2.5). However, we note that the error 
bars in the ATIC-2 data are rather large and therefore the uncertainty in Kb is also. In fact, the 
la and 2a confidence intervals are Kb G [27,66] and Kb G [2.7,90], respectively. Therefore, the 
explanations for the rise in the positron fraction and the excess in the total electron-positron flux 
is consistent with the explanation for the rise in Ti/Fe. The value of Kb — 20 is also consistent 
with the antiproton-to-proton ratio, see Sec. 3.2. 

Fig. 3.3 shows the corresponding expectation for the B/C ratio with the diffusion coefficient scaled 
proportional to rigidity according to Eq. 3.25. The CREAM data [336] do show a downward trend 
as has been emphasized recently [339], but the uncertainties are still large so we await more precise 
measurements by PAMELA which has been directly calibrated in a test beam [340]. Preliminary 
data show indeed a rise above ~ 100 GeV see Fig. 3.3. Agreement with our prediction would confirm 
the astrophysical origin of the positron excess as proposed in Ref. [243] and Chapter 2 and thus 
establish the existence of an accelerator of hadronic cosmic rays within a few kpc. 

3.6. Conclusion 

We have presented a nice and interesting test of the proposed acceleration of secondaries that could 
in principle explain the GCR lepton excesses, by means of nuclear secondary-to-primary ratios. If 
we saw a rise in such ratios, this would clearly point at the importance of this effect, allow us to 
determine the normalisation and perform cross checks with the positron analysis. Furthermore, we 
would be able to extend our knowledge about the conditions prevailing in old supernova remnants, 
for example the level of magnetic turbulence, through this additional handle. 

Unfortunately, at the moment, the experimental situation is not clear. Although ATIC-2 has 
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Figure 3.2.: The Ti/Fe ratio in cosmic rays along with model predictions — the 'leaky box' model 
with production of secondaries during propagation only (dashed line), and including production 
and acceleration of secondaries in a nearby source (solid line - dotted beyond the validity of our 
calculation) for Kb = 40 as determined by fitting to the ATIC-2 data; the 1 a and 2 a uncertainty 
bands are shown by the shaded dark grey and light grey areas, respectively. The data points are 
from ATIC-2 (triangles) [333] and HEAO-3-C3 (circles) [334]. 
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Figure 3.3.: The B/C ratio in cosmic rays along with model predictions — the 'leaky box' model 
with production of secondaries during propagation only (dot-dashed line) , and including production 
and acceleration of secondaries in a nearby source (dashed line), la (dark grey) and 2a (light 
grey) uncertainty bands for Kb = 40 as determined from Ti/Fe; solid lines for Kb = 10, 15 and 
20 as suggested by the electron-positron fit (see Sec. 2.5). The data points are from IIEAO-3-C2 
(circles) [142], ATIC-2 (triangles) [335], CREAM (squares) [336] and PAMELA (diamonds) [337, 
338]. 
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clearly observed a rise in Ti/Fe, the older HEAO-3-C3 data are in strong disagreement. Therefore, 
the ultimate test will be B/C for which the PAMELA experiment is continuously taking data and the 
official result of their analysis is therefore eagerly anticipated. Finally, the AMS-02 experiment [321] 
with its superior charge identification, broad energy range and improved statistics (due to the 
extension of the mission to 10 or even 18 years [341]) will clarify this issue once and for all. 
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4. Systematic effects in the extraction of the 
'WIVIAP haze' 



4.1. Introduction 

Local cosmic ray measurements are not the only possible probe of the galactic lepton population. 
Electrons and positrons of GeV and TeV energies produce synchrotron radiation and gamma-rays 
by scattering off the galactic magnetic fields and interstellar radiation fields (ISRFs), respectively. 
Therefore every contribution in addition to the standard electrons and positrons from supernova 
remnants (SNRs) must also reflect in galactic diffuse radio and gamma-ray backgrounds. 

We have already mentioned the possibility to constrain dark matter (DM) models invented to 
explain the positron and electron-positron excesses by gamma-ray and radio measurements. The 
bounds on the annihilation cross section from diffuse measurements are at the moment still at 
least an order of magnitude above the expectation for a thermally produced weakly interacting 
massive particle (WIMP). These studies are mostly restricted by modelling of the astrophysical 
backgrounds, like diffuse emission. It is therefore necessary, to improve our knowledge of these 
backgrounds to be able to subtract them from the data and uncover possible, exotic contributions. 

Some studies try to abbreviate this process by invoking proxies for different physical processes 
that contribute to the backgrounds. Following the assumed correlation, the data can be fitted 
for and the background subtracted. One particular type of these studies has found an excess in 
microwave sky maps - the 'WMAP haze' (see Sec. 1.3.4). A crucial ingredient of both studies 
[232, 234] that identify a haze is the extrapolation of the morphology of the synchrotron radiation 
template from 408 MHz to the WMAP bands at 23 (K), 33 (Ka), 41 (Q), 61 (V) and 94 (W) GHz, 
i.e. over two orders of magnitude in frequency. In fact the spatial distribution of the radiating 
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CR electrons is likely to differ significantly given their energy dependent diffusive transport in the 
Galaxy. Instead of attempting such a bold extrapolation, other studies, including the analysis by 
the WMAP collaboration [237], employ the K-Ka difference map as a tracer of synchrotron emission 
(despite some contamination by free-free emission and an anomalous component which has been 
interpreted (see, e.g. [342]) as spinning dust [228]). However although both maps are dominated 
by synchrotron radiation, such a template could also contain any unidentified radiation, such as a 
possible haze, and therefore cannot exclude it. 

CR transport in the Galaxy is dominated by diffusion through interstellar magnetic fields with an 
energy- dependent diffusion coefficient D(E) = DqE^ where 6 = 0.3 . . . 0.7 [275]. Taking the energy 
loss rate b{E) = dE/dt = h^E"^ as is appropriate for synchrotron and ICS, the diffusion length A is 



for the standard values Dq = lO^^cm^s"^ and 6o = 10^^^ s~^ [275]. Therefore, the distance that 
GeV energy electrons can diffuse is comparable to the kpc scale on which the source distribution 
varies; moreover it changes by a factor of 2.4 (1.5) for 5 = 0.3 (0.7) in the energy range ~ 4 — 50 
GeV (corresponding to peak synchrotron frequencies between 408 MHz and 50 GHz for a magnetic 
field of 6/iG). As a consequence the ~ 50 GeV electrons will trace the source distribution much 
better than the ~ 4 GeV electrons which diffuse further away from the sources and wash out their 
distribution. The synchrotron map at 408 MHz cannot therefore be a good tracer of synchrotron 
radiation at much higher, in particular WMAP, frequencies. Relying on such a crude extrapolation 
of the morphology of synchrotron emission can thus potentially introduce unphysical residuals. We 
estimate these by simulating synchrotron sky maps at 408 MHz and the WMAP frequencies and 
feeding these into the template subtraction process [233]. We show that this leads to residuals of the 
same order as the claimed haze, which can in fact be matched for a particular source distribution in 
the galactic disk. We conclude therefore that the WMAP haze might be an artefact of inappropriate 
template subtraction rather than evidence of dark matter annihilation. 
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4.2. Template Subtraction 

The subtraction method is based on a multihnear regression of the CMB subtracted WMAP data 
using foreground templates for free- free (f), dust correlated (d) and synchrotron emission (s). Tech- 
nically this can be achieved by assembling the maps represented each by a vector of all pixels, that 
is f, d and s, into one 'template matrix': P = (f , d, s). (The template for the haze, h, is appended 
later.) The pseudo-inverse allows the determination of the coefficients a = P"^w that minimise 
the = llw — Pa|p/cj^ for the different templates at the WMAP frequencies; a is the mean 
measurement noise in each frequency band. For details see Ref. [233]. 

Since we are interested only in the effect of the electron diffusion on the subtraction of the 
synchrotron foreground we do not use the free-free and dust templates or radio sky maps that are 
strongly affected by local structures such as Loop I [343] . Instead we simulate both the synchrotron 
sky map at 408 MHz and the sky maps in the WMAP frequency range with the GALPROP code [126]. 
We adopt the same mask as in Ref. [233] which excises pixels along the galactic plane, around radio 
sources and in directions of excessive absorption. 

To allow comparison with the results of Ref. [233] we apply the same fitting procedure over the 
whole sky. In order to determine the magnitude of the 'haze' we append a template h = (1/6 — 1 /9q) 
to P where = \fW^r^ is in galactic coordinates and = 45°. This corresponds to the "FS8" fit 
performed in [233] and adding the haze back to the residual maps gives the "FS8 -|- haze" maps. 
We determine the latitudinal profile of the residual for £ = 0° south of the galactic centre direction. 
As our simulated maps do not contain any localised structures, we do not need to divide the sky 
into several regions and fit them independently, as was done with the "RG8" fit [233] . We have 
checked explicitly that doing so does not change the profiles of the residual intensity or the spectral 
indices. 

We have checked that our procedure gives a residual 'haze' in agreement with Ref. [233] when 
we subtract the 408 MHz survey sky map from the WMAP sky maps. Although with the CMB 
estimator "CMB5" we find a residual intensity of the same magnitude at 23 GHz, its spectral index 
of about —0.7 is somewhat softer than in Ref. [233]. 
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4.3. Calculation 



The transport of CR electrons is governed by a diffusion-convection equation (cf. Eq. 1.39), 

-=V-{D^^Vn-.n) + -pD,,-^n--[pn--{V..)n)+q, 

where ndp is the number density of electrons with p G [p,p + dp], D^x = Dqxx{p/ ^GeVY is the 
spatial diffusion coefficient, v is the convection velocity, Dpp is the momentum diffusion coefficient 
and q is the source power density. This equation is numerically solved with the GALPROP code 
v50 . Ip in two dimensions, that is assuming azimuthal symmetry around the galactic centre and 
enforcing the boundary condition n = on a cylinder of radius R — 20 kpc and half-height -Smax 
(see below). 

The source power density q factorises into a source energy spectrum q^E"^ and a spatial variation 
^^j,^g-2/zscaic -with Zgcaie = 0.2 kpc. For the radial part we consider two possibilities. The distribution 
of SNRs is expected to be correlated with that of pulsars which is inferred by Lorimer to be [344] 

\ 2.35 

^ > -r/1.528kpc 



CTLonmer(r) = 64.6 — J ^-r^-^^^^VC (4 

However, the determination of pulsar distances from their rotation measures relies on knowledge 
of the thermal electron density throughout the Galaxy and different distributions lead to different 
functional forms for the inferred radial variation of the pulsar density [345]. Therefore we also 
consider an exponential source distribution 

fTexp(r) = aoe-'-/2kpc, (4.2) 

following Refs.[346] and [347]. 

The normalisation Dqxx^ the scale height ^max of the CR halo and the spectral index 5 of the 
diffusion coefficient are usually determined from measurements of CR nuclei and nuclear secondary 
to primary ratios (see Sec. 1.2.7). The measurement of CR 'chronometers' like ^^Be/^Be is still not 
precise enough to break the degeneracy between Dqxx and Zmax (see Sec. 1.2.7), so we vary Zmax 
between 4 kpc and 8 kpc and vary Dqxx only a little, checking that we have rough agreement with 
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the measured fluxes of nuclei and nuclear secondary-to-primary ratios. On theoretical grounds [103] 
one expects a spectral break in the diffusion coefficient at ~ 1 GeV. We fix the break energy to 1 
GeV and vary 6i and 62 below and above the break (keeping 5i > 62), again trying to satisfy all 
local CR measurements. 

The source electron spectrum is usually assumed to have a break around 4 GeV so we fix the 
electron source normalisation qoao and the spectral indices ai and 02 below and above the break 
by fitting the propagated flux to the electron spectrum as measured at Earth [348, 169]. We apply 
Solar modulation in the spherical approximation [60] with a median potential of (p = 550 MV. 
Reacceleration and convection play a role at energies below 10 GeV and are therefore important 
for the 408 MHz map. For the Alfven velocity v\ which determines the strength of reacceleration 
via Dpp oc we consider the range 0-50kms~^ GALPROP assumes the convection velocity to 
vary linearly with distance from the galactic plane and we vary the slope dvcom/dz between and 
20kms-ikpc"^ 

Since the random component of the galactic magnetic field is known to dominate over the regular 
component [349], we neglect the latter. For the radial dependence we adopt the usual exponential 
fall-off where the radial scale p and the (perpendicular component of the) field strength Bq at 
the galactic centre are chosen to reproduce the 408 MHz sky map [350]. Although it was initially 
believed [198] that an exponential dependence on z could give a satisfactory fit to the 408 MHz lati- 
tude profile, the galactic field model was later refined [351] by considering different, non-exponential 
behaviours which in fact give better fits. We therefore apply the method described in Ref. [352] of 
determining the emissivity dependence on r for galactic longitude £ = ±180° (towards the galactic 
anti-centre). With an estimate for the electron density this translates into a z-dependence of the 
form a + 6exp [(— this is iterated to convergence where we find a/b = 0.27, ^ = 0.51 
and K = 0.68. 

4.4. Results 

4.4.1. Lorimer source distribution 

The parameters of the diffusion model, the magnetic field and the electron source spectrum have 
been adjusted as described above and the values are shown for the Lorimer source distribution. 
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Table 4.1.: Parameters of source and diffusion models. 



Lorimer exponential 



distribution c.f. Eq. 4.1 c.f. Eq. 4.2 

ai,a2 1.2,2.2 1.2,2.2 

Doxx 5.75 X 10^8 cm^ s~^ 5.75 x lO^^ cm^ s" 

Zmax 4 kpc 8 kpc 

Si,S2 0.34,0.34 0.1,0.4 

VA 50kms^^ 36kms~^ 

d^conv/d-z; 10kms~^kpc~^ 15kms~^kpc^^ 

Bo 6.3 6.8 /iG 

p 5 kpc 50 kpc 



Eq. 4.1, in Table 4.1. 

The electron flux measured locally and at the positions {(r,z)} = {(1,0), (4,0), (0,4)} (in kpc) 
are shown in Fig. 4.1. We note that close to the galactic centre the electron flux responsible for 
synchrotron radiation at 408 MHz is not only much softer but also suppressed by over an order 
of magnitude with respect to its locally measured value. Fig. 4.2 shows the latitudinal proflles of 
the synchrotron radiation at 408 MHz; in general, the fit is good for b < 50° but underestimates 
the emission at larger latitudes. It has been shown [352] that this can potentially be overcome by 
increasing the scale height of the synchrotron emissivity at larger galactic radii. The remaining 
discrepancies between the simulated and measured profiles are probably due to the assumption of 
rotational symmetry. This leads to an underestimation of the synchrotron radiation along tangents 
of the spiral arms and an overestimation between them. For example, the Carina arm is tangent at 
75° and the Sagittarius arm at —40°, so both i = +60° and £ = —60° are between spiral arms and 
thus slightly overestimated, in particular in the galactic plane. It is also clear that point sources 
(that have not been subtracted from the 408 MHz data) are not accounted for in our calculation 
(e.g., Fornax at £ ~ 120°, b ~ -57°). 

The sky map of the residual r{i, b) (Fig. 4.3) shows a deficit for \i\ < 40° and \b\ < 20°. Further 
away from the galactic centre direction there is a slight excess. The residual specific intensity 
(Fig. 4.4) is of opposite sign but its absolute value is of the same order of magnitude as the 'haze' 
at 23 and 33 GHz. 
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Figure 4.1.: Left: The electron (plus positron) flux measured locally by AMS-01, CAPRICE, 
HEAT [348] and Fermi-LAT [169], compared with the expectation for the Lorimer source dis- 
tribution, Eq. 4.1; the dotted line is the calculated interstellar flux while the solid line is its Solar 
modulated value (with = 550 MV). The dashed vertical lines show the energy corresponding to 
peak synchrotron frequencies of 408 MHz and 23 GHz for the local magnetic field. Right: The 
calculated electron (plus positron) flux at the positions {(r, z)} = {(1,0), (4,0)} (in kpc). 
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Figure 4.2.: The calculated latitudinal profile of galactic synchrotron radiation at 408 MHz (black 
dashed line) for galactic longitudes |£| = 0°, 30°, 60°, 90°, 120° and 150°. The red (blue) solid line 
is the observed profile [350] for positive (negative) i. 



4.4.2. Exponential source distribution 

For the exponential source distribution, Eq. 4.2, the electron fluxes are shown in Fig. 4.5. Close to 
the galactic centre, it is larger by about an order of magnitude than measured locally and slightly 
harder. The latitudinal profiles of the synchrotron radiation at 408 MHz are shown in Fig. 4.6. 
The residual sky map contains a roughly spherical excess around the centre of the map, although 
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Figure 4.3.: Residual sky map in galactic coordinates for the Lorimer source distribution, Eq. 4.1. 
The grey scale is linear from —4 to -|-4kJysr~^ (corresponding to —0.25 to -1-0.25 mK at 22.8 GHz). 
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Figure 4.4.: Latitudinal profile of the K band residual outside (solid curve) and under (dot-dashed 
curve) the mask at 23 GHz (left panel) and 33 GHz (right panel). The square (circle) data points 
are the 'haze' as extracted in Ref. [233] ([239]). The dotted line shows the extrapolated emission 
at 23 (33) GHz from scaling the simulated 408 MHz emission and the dashed line shows the actual 
simulated 23 (33) GHz emission. 



somewhat more extended in longitude than in latitude (see Fig. 4.7). 

We note that the systematic uncertainty of the residual intensity (as determined from real sky 
maps) induced by chance correlations between the 'haze' template and the CMB has been estimated 
in Ref. [233] and can be read off their Fig. 8 as ±11.8hkJysr-i (±23.7hkJy sr^^) in the 23 GHz 
(33 GHz) band. We therefore allow for an offset of our calculated residual relative to the 'haze' 
template in this range when fitting the residuals from real sky maps. The residual intensity (Fig. 4.8) 
matches the claimed WMAP haze. 

To compare our results to those of Ref. [233], we also determine the average spectral index (for 
details see Appendix B) in a region south of the galactic centre, b G [—45°, 0°], £ € [—25°, 25°]. The 
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Figure 4.6.: Same as in Fig. 4.2, but for the exponential source distribution, Eq. 4.2. 



cofour maps of spectral indices scaled by intensity are shown in Fig. 4.9, both for the synchrotron + 
residual and for the residual alone. Not only is the synchrotron emission much more disk-like than 
the residual, but the spectral index of the residual is also considerably harder than the synchrotron 
spectral index. This is to be compared with Fig. 7 of Ref. [233] which exhibits the same qualitative 
behaviour. 

Furthermore, we show the spectral index for the unmasked pixels in the region south of the galac- 
tic centre (as defined above) as a function of latitude in Fig. 4.10, again both for the synchrotron 



+ residual and for the residual alone. With average indices of (/3^ 



23,33/ 



-0.44 for the residual and 



of (/Sfs 33) = —1.03 for the residual + synchrotron in this region, we find that the residual index 
is harder than the synchrotron index by 0.6, which is in excellent agreement with the findings of 
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Figure 4.7.: Residual sky map in galactic coordinates for the exponential source distribution, 
Eq. 4.2. The grey scale is linear from —4 to -|-4kJysr~^ (corresponding to —0.25 to +0.25 mK at 
22.8 GHz). The contour lines are logarithmically spaced in intensity between 0.02 and 2kJysr~^. 
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Figure 4.8.: Same as in Fig. 4.4, but for the exponential source distribution, Eq. 4.2. We have added 
an offset to the calculated residual of +11.8 hkJy sr~^ (+23.7hkJy sr~^) in the 23 GHz (33 GHz) 
band reflecting the systematic uncertainty from chance correlations between the 'haze' template 
and the CMB. 

Ref. [233]. The values of our different model parameters are shown in Table 4.1. 

4.5. Discussion 

To qualitatively understand these results, consider the longitudinal profile of the synchrotron inten- 
sity I{i,b); for simplicity let us constrain ourselves to the galactic plane, i.e. 6 = 0. The intensity 
in any direction i is given by the integral of the synchrotron emissivity over the line of sight and 
this samples the radial distribution of the relativistic electron density in the range r £ [dsini, R], 
where d is the distance of the Sun from the galactic centre. Since the fitting procedure minimises 
the square of the difference in the maps, the sign and size of the residual is determined not by the 
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Figure 4.9.: Colour maps of spectral indices between 23 and 33 GHz defined in Eq. B.l scaled by 
the 23 GHz intensity for synchrotron + residual (top panel) and residual only (bottom panel). 
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Figure 4.10.: Spectral indices of the unmasked pixels in the region south of the galactic centre 
{b G [— 45°,0°], £ € [—25°, 25°]) as a function of latitude for residual -|- synchrotron (large beige 
circles) and residual alone (small blue squares). The average spectral indices, /3s and /3h) are shown 
in the upper right corner. 
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absolute difference but by the difference in tlie radial slope of the emissivity e(r) at 408 MHz and 
the WMAP frequencies. The difference in the slopes reflects the energy dependence of the electron 
diffusion — higher energy electrons lose their energy more quickly, hence their emissivity traces the 
source distribution more closely than does the emissivity of low energy electrons. 

For the pulsar source distribution, the low energy electrons peak at the galactic centre whereas 
the high-energy electrons peak further away (see left panel of Fig. 4.11). This leads to a deficit for 
small radii (translating to small longitudes) and a slight excess further away from the galactic centre 
(see also Fig. 4.3). For the exponential source distribution the radial distribution of synchrotron 
emissivity is steeper at higher energies. The template subtraction therefore yields a residual with 
an excess around the centre direction and a deficit further away along the galactic plane (see right 
panel of Fig. 4.11). 

We note that the size and morphology of the residual is thus sensitive not only to the source 
distribution but also to the parameters of the diffusion model. For instance, decreasing the Alfven 
speed below the value given above reduces the importance of reacceleration, and therefore effec- 
tively limits the number of GeV electrons around the galactic centre where otherwise energy losses 
dominate. 

4.6. A Related Study 

A related study [353], published at the same time as the present work [3], also tries to determine 
whether the 'WMAP haze' could have a purely astrophysical origin and how to distinguish between 
this and a DM explanation. As there are some important differences, we will briefiy comment on 
this work. 

At first sight, the astrophysical explanation invoked in Ref. [353], a large number of supernova 
remnants at the galactic centre, might look not too different from our model with the exponential 
source distribution. However, while we also agree on the importance of diffusion-loss steepened 
electron spectra for producing the haze there is a major difference between our approaches — 
while the authors of Ref. [353] consider the haze to be physical, we argue that it might in fact 
be an artefact of the foreground subtraction. Our models are also more constrained insofar as 
we reproduce the observed radio emission at 408 MHz and match the direct measurements of the 
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electron spectrum at the Solar position. Furthermore, we allow for spatial dependence of the B 
field, and convection and reacceleration of cosmic ray electrons, which are all essential in order 
explain all these datasets simultaneously. 

4.7. Conclusion 

We have investigated systematic effects in WMAP foreground subtraction stemming from the naive 
extrapolation of the 408 MHz map. To this end we have considered two illustrative cosmic ray dif- 
fusion models assuming different source distributions, the first one based on a pulsar survey, and 
the second one exponential in galactocentric radius. Both models are able to reproduce the syn- 
chrotron radiation at 408 MHz, the locally measured electron flux and are furthermore consistent 
with nuclear cosmic ray fluxes and secondary-to-primary ratios. When our 'foreground' 408 MHz 
map is subtracted from the 23 GHz map, we find a residual whose size and morphology depends on 
the source and diffusion model adopted. Thus the energy-dependent diffusion of relativistic elec- 
trons makes the 408 MHz sky map a bad tracer of synchrotron radiation at microwave frequencies, 
as had been suspected earlier [354]. Such a template subtraction produces a residual, which for 
certain values of the source and propagation model can reproduce the 'WMAP haze' in intensity 
and spectrum. 

The residual obtained from the exponential source distribution does not perfectly reproduce 
the morphology found in Ref. [233] (although it is not disk-like but rather clustered around the 
galactic centre). However, a quantitative assessment of the discrepancy is not straightforward, 
mainly because Ref. [233] does not provide any objective measure, e.g. the ellipticity of equal 
intensity contours. On the other hand, even the numerical GALPROP model we employed for our 
analysis is very likely too simple to fully capture the complexity of synchrotron emission in the 
Galaxy. For instances, not only the source density but also the galactic magnetic field is supposed 
to be correlated with the galactic spiral arms, which will break the symmetry in r (and hence 
in i) and can therefore considerably modify the morphology. Furthermore, much of the 'diffuse' 
synchrotron emission from the disk may originate in the shells of old supernova remnants which 
have grown very large in their radiative phase [355]. Exactly the same argument concerning the 
energy-dependent diffusion length that we applied to the cosmic ray source distribution can be 



128 



applied to such localised structures too. Therefore the 408 MHz survey sky map is not expected to 
trace the emission from the latter at higher frequencies either. One can easily imagine that such 

localised structures (of which Loop I may be a nearby example) can at least in part modify the 
morphology of the residual and bring the simulated map into agreement with the one determined 
from the subtraction of real templates. 
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Figure 4.11.: The (scaled) synchrotron emissivity at 408 MHz and 23 GHz, and their difference, for 
the Lorimer (left panel) and the exponential (right panel) source distribution. 
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Epilogue 



Nearly two years after the excesses in the positron fraction and the total electron-positron flux 
were reported by the PAMELA and Fermi satellites respectively, the dark matter explanations 
proposed for these 'anomalies' seem to be almost ruled out. One the one hand, the necessarily 
powerful injection of additional electrons and positrons from annihilating/decaying dark matter in 
an extended halo is in conflict with radio/microwave and 7-ray observations. One the other hand, 
a broad range of astrophysical effects have been identified which can account for these excesses 
and would naturally dominate over plausible dark matter signals. Furthermore, signatures that 
have been believed to be typical of dark matter have in fact been found to be rather generic: 
features in the total electron-positron flux at tens and hundreds of GeV can easily be induced by 
the discreteness of astrophysical sources. Even sharp shoulders in the energy spectrum can be due 
to a nearby source with a hard spectrum, like a pulsar. Conversely, it appears that an arbitrary 
spectral signal, even if it has nothing to do with dark matter, can be fitted in a dark matter scenario 
by harnessing the multitude of free or uncertain model parameters. 

Prospects for dark matter indirect detection in charged lepton cosmic ray channels looks therefore 
rather bleak at the moment. In the long run with the advent of advanced satellite detectors such 
as AMS-02, we might hope for a better understanding of the astrophysical backgrounds down to 
perhaps even the per cent level. This would allow us to subtract them and uncover possible dark 
matter contributions, which in WIMP scenarios would naturally be expected to contribute at this 
level. Other high energy channels might, however, provide better prospects. Photons, for example, 
do not only sample a much larger volume than locally measured charged leptons which, as we 
have seen, originate within a few kiloparsecs of the Solar system. Even more importantly, because 
of their rectilinear propagation, it is possible to focus on specific targets that have little or no 
astrophysical backgrounds, e.g. the Galactic halo at high latitude or prominent substructures such 
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as dwarf spheroidal galaxies. The Fermi satellite as well as the forthcoming Cerenkov Telescope 
Array will provide important data in this connection. 

Of course, the attention given to the lepton excess because of the claimed connection with particle 
dark matter has played an important sociological role in helping the astroparticle community to 
appreciate the complexities of galactic cosmic ray physics. The first hints of an anomalous positron 
fraction were reported over 30 years ago, but only a few prescient works have taken this signal 
seriously and tried to explain it, e.g. by the contribution from nearby pulsars. We hope very much 
therefore that the momentum which cosmic ray physics has gained of late will foster improvements in 
their astrophsyical modelling and help us to be open-minded and alert to other possible anomalies 
that might show up in the data. An example could be the recent observation by the CREAM 
detector of different spectral indices for protons and helium nuclei. Possible interpretations range 
from different source populations for protons and heavier nuclei, to secondary production and 
acceleration of helium in supernova remnants, in much the same way as has been suggested in this 
thesis for secondary electrons and positrons. 

Although dark matter thus seems to be ruled out as an explanation for the lepton excesses 
this puzzle is far from being solved. Not only is it unclear which sources are responsible for the 
additional positrons, one might also wonder whether these excesses are local (i.e. particular to the 
environment of the Solar system) or global (i.e. similarly present in other regions of the Galaxy). 
In this context, we would like to stress the two-fold role of the nuclear component of galactic cosmic 
rays: Firstly, if a rising secondary-to- primary ratio (e.g. B/C) is observed, this will clearly show 
that the acceleration of secondaries is the likely cause of the rise in the positron fraction. Secondly, 
as the local fluxes of stable nuclei obtain contributions from larger distances in the Galaxy than do 
charged leptons, such a rise would imply that the secondaries are also produced and accelerated in 
cosmic ray sources elsewhere in the Galaxy. 

It has also been suggested that anisotropics in the arrival directions of charged leptons in cosmic 
rays can help identifying the origin of the additional positrons. We believe however that the 
uncertainties here are too large and that anisotropics are therefore unlikely to settle the issue. As 
mentioned above, models that explain the increasing positron fraction and excessive total electron- 
positron flux by a (few) nearby pulsar (s), usually find anisotropics of at most a few per mil and 
rely on assumptions concerning the homogeneity of energy losses and diffusion. It is therefore 
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conceivable that small variations of the diffusion model can drastically alter the predictions. A 
systematic orientation of the local magnetic fields, for example, can shift the predominant arrival 
directions in such a way that no correlation with known objects can be made; or it can lead to a 
spurious anisotropy even for an isotropic source distribution. 

The only possible way to investigate the high energy electrons and positrons elsewhere in the 
Galaxy is to look for their secondary radiation, i.e. radio/microwaves synchrotron radiation on 
Galactic magnetic fields and 7-rays from inverse- Compton scattering on interstellar radiation fields. 
As wc have seen, a prime application of this is again the study of dark matter annihilation or decays 
into leptons. Of course, this raises the important question of how to improve the modelling of diffuse 
galactic backgrounds which are of general interest also for other fields of physics, e.g. studies of 
the polarisation of the cosmic microwave background (CMB). At the moment, large efforts are 
being made by the Fermi collaboration in explaining the astrophysical contribution to the 7-ray 
sky. Their approach accounts not only for the ICS contribution using the GALPROP model but 
also uses a very flexible fitting of ISM densities in galactocentric rings to the 7r° decay generated 
7-rays. First results show that the residuals thus obtained are very small but the obvious question 
is whether this is because of too many model degrees of freedom. Theoretical models might not 
have fewer parameters but can at least provide some insight into their physical significance. At the 
moment, however, such models are not able to incorporate, let alone predict, local structure such as 
the radio 'loops' clearly visible in sky maps over a wide range of frequencies. Furthermore, a high 
degree of degeneracy is to be expected, for example between the GCR source distribution and the 
structure of the magnetic fields or ISRFs. Attempts at solving the inverse problem, i.e. determining 
the three-dimensional distribution of emissivity would help in extracting the rclativistic electron- 
positron density throughout the Galaxy as a function of the assumed radiation backgrounds and 
would bring us closer to identifying individual sources of GCRs. 

On a more general level and irrespective of whether the explanations we have suggested for 
the apparent signals will be confirmed by future data, the moral is that we need both better 
understanding of the astrophysical sources and of the propagation of galactic cosmic rays. 

Many of these efforts will probably be data-driven. For example, Fermi-LAT is presently closing 
a crucial energy gap in GeV 7-rays that lACTs could not address. This will allow looking for 
a specific bump in the 7-ray spectra of SNRs that would implicate pion production and hence 
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hadronic acceleration processes. In general, many arbitrary extrapolations from TeV energies as 
well as the underlying physics models will be tested. In addition, the improved statistics with 
respect to, e.g. the EGRET satellite will allow sampling of the spectra of even comparatively 
faint and distant sources. Furthermore, the superior spatial resolution will allow and/or necessitate 
going beyond zero- or one-dimensional source models and considering the internal structure of the 
sources. It might turn out that some common assumptions about the sources are in fact prejudices 
and cannot be sustained any longer. 

With respect to improving the cosmic ray propagation models, what is probably needed is a 
completely new approach. Many important steps towards realistic models have been taken in the 
last couple of years. For example, the inclusion of Klein-Nishina corrections to the interactions 
of electrons and positrons with radiation fields has been recognised to lead to spectral effects at 
the highest energies. Two other improvements aim at a more realistic modelling of the ISM and 
its effect on cosmic propagation: Firstly, the structure of the large-scale ordered magnetic field is 
being given due importance for synchrotron radiation studies. Secondly, the spatial dependence of 
energy loss rate and diffusion coefficients is beginning to be accounted for in numerical propagation 
codes. However, what these models lack is a clear physical intuition for how these quantities should 
vary across the Galaxy. Here, the hitherto followed approach of calibrating diffusion models by 
measurements of local cosmic ray nuclei necessarily reaches its limits: Stable nuclear cosmic rays 
average over a large volume of the cosmic ray halo (typically several kiloparsecs at GeV energies) and 
can therefore not resolve structures at smaller scales. At the highest energies, however, local sources 
become more important and therefore the large bulk of the galactic volume remains inaccessible. 
Generally speaking, it is difficult to imagine how an enlarged set of parameters that necessarily 
accompanies every more complex propagation model could be constrained by the same number of 
observations. 

One possible resort could be to devise new sets of observables which are particularly sensitive 
to some of the parameters. We have already encountered one such example in our analysis of the 
systematic effects of template subtraction. The degeneracy encountered here was that between the 
source distribution and the magnetic field configuration. If one can neglect exotic contributions, 
a new observable like the relative differences of maps at different frequencies could be used to pin 
down one particular combination of source distribution and magnetic fields. 
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The more radical approach would be to search for a new propagation model altogether. This 
would need to provide a self-consistent explanation for many different mechanisms that today's 
heuristic models invoke in an ad hoc fashion. The combination and description of cosmic rays, 
magnetic fields and the interstellar medium in one framework will require deeper theoretical efforts 
as well as immense computational power. One can hardly overestimate the difficulties such a 
program would face, nevertheless it would be a fitting challenge to undertake on the eve of the 
centenary of the discovery of cosmic radiation by Victor Hess in 1912. 
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A. Moments of Functions of Random Variables 



Suppose we have two random variables, x and y with the probabiUty densities fxix) and fy{y). A 
function (f){x, y) of x and y is also a random variable, 



z = (f){x,y), — = (l)'{x,y). 



(A.1) 



Let's further assume that <p has an inverse function with respect to x, 
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B. Determination of Spectral Indices From 



Synchrotron Sky Maps 



In general, a spectral index /3(x) between two different frequencies, z^i and 1^2, can be defined for 
each given pixel x by assuming a power law behaviour of the specific intensity, I(z^, x). 



However, it turns out that the template method applied to the WMAP data and the 408 MHz sky 
map leads to a residual with negative intensities for some pixels (see e.g.. Fig. 6 of Ref. [233]), 
partly due to over-subtraction and partly because the sky maps are mean-subtracted. We also find 
negative intensities for some pixels when applying the template subtraction to our mock microwave 
data and radio template. This does not necessarily imply that the residual is not physical but that 
a global offset A/(z^) exists between the residual intensity, /', as determined from the template 
subtraction and the intensity of the actual, possibly physical residual, /: 



This makes the determination of the spectral index non-trivial. 

At first sight, the analysis presented in Ref. [233] seems to avoid this difficulty by determining 
the average spectral index in the region south of the galactic centre from the average ratio r' of the 
intensities at two different frequencies vi and 1^2, e.g. I'l = 23 GHz and 1^2 = 33 GHz. This ratio 
can be determined from a scatter plot of the pairs of residual intensities {Il-^^jl'^.^} (as determined 
from the template subtraction), to which a straight line, = r'l'j^^ + I^^; is fitted, allowing 

for the ordinate offset I'^^ because of the unknown global offset AI(z/). The average spectral index 




(B.l) 



A/(zy) = /(i/,x) - /'(i/,x) . 



(B.2) 
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{Pvi,v2) defined by this procedure is then simply log (r')/ log (1^2/1^1). Alternatively, if the spectral 
index is determined from a scatter plot of the actual residual intensities {lu^^In^}, then there is 
no ordinate offset, so the straight line is Ii,^^^^) = rly^ and the actual average spectral index 
{1^1^1,1^2) = log (r)/ log {V2/U1). 

In general, these two descriptions cannot be expected to give a similar spectral index. Even 
assuming that with an appropriate 'haze' template h the amount of over-subtraction is much 
smaller than the offset due to the use of mean-subtracted maps, the answer is in general different. 
In this case, the offset is simply the mean over the n pixels, A7(z^) = (/(z^, x)). The coordinate 
system {11^,11^} is therefore centred at the centre of gravity of the data {/(i/i, x), /(f2, x)}, and 
the ordinate offset I'^^ is zero. As usual, the slope of the linear regression = r'l^^ is 

, ^ E»4(^04(^») - n(4(x,))(4(x,)) 

Unless the covariance of I^^ and I^^ is much larger than the product of their mean values, which 
is for example the case if the spectral index is constant in the region of interest, this is in general 
difi^erent from the slope r of the straight line I^,^ {In^ ) = rln^ , 

However, since we cannot determine the offset A/(z/) from data, we need to define an offset A7(i/). 
We choose it to be: 

= min [I'(i^,x)] , (B.5) 

such that the intensity is always positive, allowing us to define the spectral index in each pixel. 
(The exact value chosen for A/(z/) is actually (1 -|- 10~^) min [/'(z^, x)] to prevent the spectral index 
from diverging in the pixel where 7(23 GHz) attains its minimum.) 
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